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I.  INTRODOCTIOH 


A.  OVERVIEW 

Briefly  stated,  direct  time-domain  techniques  for 
transient  radiation  and  scattering  are  broadly  concerned 
with  the  response  of  soae  object  of  interest,  either 
equipment  or  a  target,  to  non-monochromatic  electromangetic 
sources.  The  analysis  is  done  in  the  tine-domain,  with  time 
treated  as  an  explicit  indepandant  variable  in_  most  methods. 

In  the  past,  frequency  domain  analysis  has  been  the 
primary  means  of  analyzing  non-monochromatic  EH  phenomena. 
Time-domain  solutions  in  closed-form  were  almost  impossible 
to  obtain,  except  in  certain  very  simple  cases.  The  amount 
of  computational  work  required  was  correspondingly 
significant.  Ho  practical  aides  were  available  to  assist  in 
performing  the  operations  prior  to  the  development  of  fast 
computer  systems. 

Another  limitation  was  the  requirement  that  the 
time-domain  response  be  as  nearly  an  impulse  as  possible  and 
thus  possess  a  broad  band  of  frequencies.  Hardware  was 
required  to  produce  very  short  pulse  waveforms. 
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The  situation  changed  radically  with  the  widescale 
introduction  of  computers  capable  of  performing  rapid 
computations  involving  large  numbers  of  iterations  in 
relatively  short  periods.  This,  coupled  with  the  rapid 
maturing  of  short-pulse  hardware  technology,  soon  made 
possible  the  application  of  time-domain  techniques  to  a  wide 
variety  of  eletromagnetic  problems. 

Siller  and  Landt  [Ref.  1]  list  some  advantages  of 
direct  time-domain  solutions  over  frequency  domain 
treatments  of  transient  problems.  They  are: 

1.  Greater  solution  efficiency  for  many  types  of  EH 
problems. 

2.  Sore  convenient  handling  of  non-linearities. 

3.  Improved  physical  insight  and  interpretablility  of  EH 
phenomenon. 

4.  Availability  of  wide-band  information  from  a  single 
calculation. 

5.  Opportunity  to  isolate  interactions,  using  time  range 
gating  (e.g.,  pulse  reflections  from  distant  objects 
not  of  interest,  etc.,). 

The  major  disadvantage  is  the  increased  complexity  of 
the  computer  implementation  of  a  time-domain  program  as 
opposed  to  one  involving  frequency  domain  techniques.  There 
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is  also  a  corresponding  difficulty  in  the  software 
development  and  its  practical  use.  Finally,  the 
computational  time  may  be  significantly  greater  than  that 
required  of  other  techniques. 

The  applications  for  time-domain  techniques  are  diverse. 
Transient  analysis  can  be  used  to  directly  evaluate  the 
transient  characteristics  of  new  antenna  designs, 
determining  such  parameters  as  driving  point  currents, 
radiated  far  fields  and  near  fields,  and  assisting  in 
determining  antenna  gain.  These  methods  may  also  be  used  in 
determining  the  degree  of  coupling  due  to  electromagnetic 
pulses  (EBP)  generated  because  of  lightening  during  a 
nuclear  detonation.  Pinally,  time-domain  techniques  may  be 
directly  applied  in  the  solution  of  the  inverse  scattering 
problem  involving  the  receiving  and  scattering  properties  of 
a  large  category  of  target  types. 

ks  tine- domain  transient  techniques  are  further  refined, 
it  can  be  expected  that  a  graater  number  of  uses  will  be 
identified. 

3.  2ETIBB  OP  TBiHSXEIT  BB  1N1L7SIS  TECHNIQUES 

There  are  three  broad  categories  of  techniques 
identifiable  in  the  literature  involving  transient  analysis. 


They  are  direct  (tine- domain)  techniques,  transform 
(frequency  domain)  techniques,  and  hybrid 

techniques.  [Ref.  2]  This  report  is  only  interested  in 

direct  techniques  and  will  not  review  the  others. 

1 .  Physical  Optics  Inverse  Scattering 

The  earliest  techniques  developed  involved  physical 
optics  inverse  scattering.  Kennaugh  and  Cosgriff  [Ref.  3] 
in  1958  proposed  that  the  impulse  response  of  a  scattering, 
conducting  body  is  simply  the  second  derivative  of  the 

projected  area  function  of  the  body  if  physical  optics 
currents  are  postulated  on  the  surface.  If  a  target  is 
illuminated  with  an  incident  ramp  plane  wave,  then  the 
projected  area  function  of  the  target  is  directly 

proportional  to  the  far-scattered  ramp  response.  This  is 
only  an  approximate  relationship,  however,  since  it  is  based 
on  the  assumption  of  physical  optics  currents. 

A  more  direct  and  precise  ramp  response  technique 
was  suggested  by  Kennaugh  and  Moffat  [Ref.  4]  in  1965. 
They  predicted  satisfactory  approximate  ram?  response 
signatures  could  be  obtained  with  narrower  bandwidth 

interrogating  signals  than  that  required  for  the  impulse 
response. 


19 


The  above  techniques  were  generally  applicable  for 
targets  which  were  rotationally  syaaetric  and  for  whoa  the 
direction  of  the  electroaagnetic  field  incidence  was  axial. 
In  1976,  Toung  (  Bef .  5]  deaonstrated  that  a  voluae  estiaate 
can  be  obtained  froa  the  raap  response  wavefora.  Re  further 
extended  the  raap  response  technique  to  arbitrarily  shaped 
targets  with  the  incident  field  being  other  than  axial. 

2.  "Exact1*  Inverse  Scattering  Method 

Another  aethod  which  iaproves  upon  the  physical 
optics  inverse  scattering  theory  is  the  "exact"  inverse 
scattering  aethod  first  suggested  by  Bennett  et.  al., 
[Bef.  6]  in  1974,  and  further  refined  in  1977.  [Bef.  7] 
The  inverse  scattering  problea  is  solved  using  an  inversion 
of  a  space-tiae  integral.  The  aethod  takes  into  account  the 
backscatter  response  of  the  target,  the  direct  effect  of  the 
incident  electroaagnetic  field,  and  the  correction  currents 
flowing  on  the  target  surface.  Theoretically,  b7  coapletely 
describing  the  interactions  of  the  two  components  of  the 
surface  currents,  a  close  approxiaation  of  the  target  shape 
is  possible.  The  actual  iapleaentation  is  liaited  to 
"perfectly"  conducting  scattaters.  It  is  this  aethod  which 
has  been  selected  for  use  in  the  Naval  Postgraduate  School 
Tiae  Doaain  Badar  Laboratory  (TDHL) . 
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Other  aethods  that  have  been  developed,  although  to 
a  less  satisfactory  degree  of  accuracy  than  the  inverse 
scattering  aethods,  include  radar  target  identification 
aethods.  Only  a  finite  number  of  solutions  (targets  shapes) 
are  allowed  as  opposed  to  the  nearly  infinite  number 
possible  using  inverse  scattering.  The  models  of  targets 
are  stored  for  parameter  comparison  with  actual  returns. 
Certain  parametric  aspects  of  a  wideband  illuminating 
signal,  such  as  natural  resonance  frequencies  and 
polarization  of  the  return,  ar§  used.  Other  similar  methods 
use  amplitude  and  phase  differences  or  special  properties  of 
Rayleigh  region  scattering.  Ill  these  methods  are  limited 
in  practical  application. 

Still  another  important  and  powerful  method 
considers  multifrequency  target  illumination  using  existing 
radars  to  simultaneously  illuminate  targets  with  VHP,  UHP,  L 
and  S  band  radar  frequencies  with  circular  polarization. 
[Bef.  8j  Circular  polarization  eliminates  any  polarization 
angle  dependance  the  signal  might  have.  This  approach 
allows  for  the  inclusion  of  all  radar  system  parameters 
(such  as  antenna  gain,  noise  level,  range,  etc.)  in  the 
Identification  algorithm.  The  technique  further 
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demonstrates  that  the  target  shape  information  may  be 
obtained  with  a  high  degree  of  accuracy  from  radar  systems 
correctly  tailored  to  the  problem  of  general  target 
identification. 

C.  REVIEW  or  TDBL  PROGRESS  AT  UPS 

Three  requirements  were  initially  identified.  First,  a 
method  had  to  be  chosen  from  the  techniques  listed  above. 
Secondly,  the  construction  of  a  physical  laboratory  to 
perform  the  required  measurements  in  support  of  the  chosen 
technique  had  to  be  designed  and  built.  Thirdly,  an 
alogorithm  had  to  be  developed  employing  the  selected 
technique  for  laboratory  use. 

As  stated  above,  the  "e-ract"  inverse  scattering  method 
suggested  by  Bennett  was  the  scheme  chosen.  To  support 
this,  work  began  in  January  1980  where  evaluation  of  various 
designs  to  be  used  in  constructing  a  Time  Domain  Radar 
Laboratory  at  the  Naval  Postgraduate  School  was  carried  out 
by  Capt.  C.  ».  Hammond,  OSSC,  as  a  thesis  project.  rforlc  was 
essentially  completed  with  the  erection  of  the  laboratory 
and  interfacing  of  hardware  in  September,  1930.  [Ref.  9] 

In  support  of  the  TDRL,  LCDR  Heir  Horag,  Isreali  Navy, 
began,  as  &  thesis  project,  the  developmsnt  of  an  "exact” 
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1980 


The 


inverse  scattering  algorithm  in  September, 
algorithm  was  completed  in  .larch  1981  [Bef.  10]  However, 
the  program  was  developed  in  FORTH  AH ,  a  language  not  used  by 
the  TD8L  computer  system.  A  translation  to  3ASIC  and 
adaptation  of  the  code  was  necessary  to  reduce  memory 
storage  requirements  and  to  ensure  computational  times  were 
of  reasonable  lengths. 

The  final  step  in  the  establishment  of  the  TDRL  was  the 
validation  of  the  laboratory  in  confirming  the  work  of 
Hammond  and  Morag.  This  is  one  of  the  principal  aims  of 
this  project. 

D.  THESIS  OBJECTIVES 

The  objectives  of  this  thesis  are  to: 

1.  Develop  a  working  software  program  which  will  acquire 
and  prepare  transient  response  data  from  various 
targets  for  input  to  the  time-domain  integral 
equation. 

2.  To  show  that  through  signal  processing,  the  effects 
of  the  receiving  antenna  can  be  neglected,  and  the 
response  to  arbitrary  specified  exciting  waveforms 
can  be  constructed,  or  "synthesized”,  using  PFT 
techniques. 
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3.  Validate  the  TDRL. 

The  main  effort  has  been  directed  toward  the  development 
of  a  user  oriented  computer  program  which  will  allow  the 
determination  of  the  second  objective  and  prepare  the 
transient  domain  data  for  input  into  the  numerical 
time-domain  integral  equation.  The  shape  of  the  object 
being  observed  is  displayed  both  before  and  after  correction 
factors  are  applied.  The  shape  of  the  object  was  developed 
from  the  ramp  response  of  the  impulse  function  applicable  to 
the  target.  The  program  was  also  made  compatible  with  the 
input  of  data  to  an  alternate  natural  resonance  frequency 
target  identification  technique,  using  Prony's  method,  which 
is  being  developed  by  LT  Demetrius  Papaspiridakos,  Hellinic 
Navy,  as  a  thesis  project. 

This  report  first  describes  the  physical  structure  of 
the  TDHL.  It  then  derives  the  relevant  time-domain  integral 
equations  for  both  the  radiation  and  inverse  scattering 
solution  for  simple  axisymmetrical,  perfectly  conducting, 
closed  surfaces  in  Chapter  III.  Next,  it  discusses  and 
develops  the  numerical  solutions  to  the  integral  equations. 
In  Chapter  V,  the  implementation  of  the  complete  computer 
code  is  presented.  And  finally,  measurements  that  are 
relevant  to  the  validation  and  confirmation  of  previously 

developed  concepts  are  reported. 

2h 


This  thesis  is  the  third  to  be  produced  for  this  major 
project  at  the  Naval  Postgraduate  School.  It  marks  a  point 
of  departure  in  that  essentially  all  the  preliminary  work 
which  will  allow  for  the  effective  use  of  the  TDRL  has  been 
performed  with  validation.  The  laboratory  is,  for  the  most 
part,  ready  for  applied  research  involving  transient 


time-domain  analysis. 
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INTflODOCTION 

is  stated  in  Chapter  I.,  the  primary  thrust  of  this  work 
is  to  proride  a  practical  working  alogorithm  to  serre  as  the 
principal  link  between  the  acquisition  of  transient  data  for 
targets  situated  on  the  TDRL  scattering  range,  and  the  use 
of  that  data  in  the  inrerse  scattering  algorithm. 

The  work  performed  on  this  project  is  a  direct  extension 
of  the  preliminary  work  performed  by  Hammond  in  the 
construction  of  the  Naval  Postgraduate  School  Time  Domain 
Sadar  Laboratory  (TDRL) ,  and  by  Morag  in  the  development  of 
the  numerical  integral  equation  algorithm  for  the  solution 
of  the  inrerse  scattering  problem  involving  axisymmetric 
targets. 

The  results  obtained  by  Hammond  and  Morag  are 
essentially  complementary.  However,  the  independant 
approaches  each  took  in  their  work,  particularly  in  the 
selection  of  unrelated  processing  .  machines  for  data 
processing,  made  necessary  certain  practical  modifications 
to  the  TDBL.  It  is  germane  to  review  the  basic  work  of 
Hammond  and  discuss  those  modifications  made.  Such  a  review 
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will  also  allow  for  the  setting  of  a  frame  of  referance  to 
aid  in  the  explanation  of  the  linking  algorithm  developed 
and  results  obtained  in  this  project. 

3.  PHYSICAL  DESCRIPTION  OP  TDRL 

The  physical  set-up  of  the  TDRL  is  shown  in  Figure  2.  1. 
As  can  be  seen,  the  TDRL  can  be  broken  functionally  into 
three  groups: 

1.  The  Impulse  Generator  Source  Group 

2.  The  Imaging  Plane  Group 

3.  The  Sampling  Receiver  and  Signal  Processing  Group 
Basically,  the  TDRL  is  an  open  system  designed  to  determine 
the  transfer  function  of  a  target  from  its  backscattered 
radiation.  Scattering  measurements  are  bi-static.  A 
description  of  each  functional  group  follows. 

saiiLsa.  aiflue 

This  group  consists  of  a  single  item,  a  Video 
Communications  UHF  Impulse  Generator,  Model  1000.  This 
generator  is  the  source  of  the  transmitted  pulse  and  also 
provides  the  triggering  signal  which  synchronizes  the  Signal 
Processing  Group  to  the  transmitted  signal. 

The  Impulse  Generator  is  specified  to  provide  a 
fixed  amplitude  45-  volt  pulse  with  a  maximum  rise  time  of 
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Pigara  2.1.  Time  Doaain  Radar  Laboratory 
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Figure  2.3.  Source  Pulse  froe  Generator 

370  picoseconds.  Figure  2.3  is  a  typical  wavefora  as 
aeasured  on  a  Tektronix  S-6  Saapler  Head  by  the 
aanufacturer . 

The  output  of  the  Iapulse  Generator  as  aeasured  at 
the  TDRL  as  aeasured  by  the  Operating  Systea  is  given  in 
Figure  2.4,  and  Figure  2.5  as  photographed  froa  the  DPO 
directly.  To  set  full  wavefora  presentation  on  the  aeasuring 
oscilloscope,  30dB  attenuation  was  applied  to  the  generator 
output.  Taking  the  attenuation  into  consideration,  the 
aeasured  peak  value  of  voltage  (1.4604  volts)  gives  the 
actual  output  of  the  generator  to  be  46.2  volts.  This  is 
about  3%  better  than  rated.  The  rise  tiae  of  the  output  is 


30 


TARGET:  RUN  DATE  DIST 
1  6-26-81  N/A 

NANI NUN  PEAK  UALUE  (UOLTS) 

NINIHUN  PEAK  UALUE  (UOLTS) 

RMS  UALUE  (UOLTS) 

MEAN  UALUE  (UOLTS) 

NUMBER  OF  HAUEFORMS  AUERAGED  ■ 


ANT  TGT  REMARKS 

H«  ???  GENERATOR  UAUEFORH 

1.1284 

-8.3328 

8.4214 

8.8808 

1  OPTIMIZATION  UALUE  •  8 


Figure  2.8.  source  Waveform  Measured  at  TDRL 

about  250  picoseconds/  also  better  than  specified  by  the 
manufacturer,  by  about  48*. 

Soue  notable  features  of  the  generator  used  in  TDRL 
■easureuents  are  a  variable  trigger  repetition  rate  frou  10 
Hz  to  1  9 Hz  and  an  adjustable  delay  from  50  to  150 

nanoseconds,  nil  outputs  are  matched  into  a  50-oha  load. 
The  generator  allows  rapid  acquisition  of  highly  coherent/ 
repetitive  waveforms.  This  is  an  important  feature  when 
attempting  to  aaxiaize  signal-to-noise  ratio  through  signal 
averaging. 
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Figure  2.5.  DPO  Waveform — DBF  Impulse  Generator 


The  Video  Communications  Impulse  Generator  replaces 
a  Tektronix  Type  109  impulse  generator  used  by  Hammond.  The 
completely  solid-state  electronics  of  the  Mark  1000  has 
proven  more  reliable  in  producing  a  relatively  noiseless, 
pulse-to-pulse  coherent  signal  than  did  the  reed-switched 
Type  109.  The  Type  109  had  a  fixed  pulse  repetition  rate  at 
about  720  pulses  per  second.  This  is  a  very  limiting  factor 
when  attempting  to  noise  average  a  signal. 

The  impulse  source  is  connected  to  the  transmitting 
antenna  via  a  3. 2  meter  RG-8A/U  coaxial  cable.  This  cable 
serves  both  as  a  conductor  of  the  source  pulse  and  a  short 
delay  element  to  the  transmitted  signal. 
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2.  laagiag  Plane  Group 


The  Imaging  Plane  Group  consists  of  four  elements 
the  transaitting  antenna,  the  receiving  antenna,  the 
target,  and  the  iaaging  plane. 

a.  Transaitting  Antenna 

The  transaitting  antenna  is  a  6.4-meter  endfed 
aonopole.  It  extends  vertically  froa  the  ground  plane  and 
produces  a  vertically  polarized  EH  field. 

b.  Receiving  Antenna  —  Transverse  Electroaagnetic 

Horn  (TEH) 

Haaaond  explored  a  variety  of  receiving  antenna 
types  to  be  used  in  the  TDRL.  The  aost  useful  is  the  TEH 
horn  antenna  in  that  it  has  the  following  advantages  over 
the  other  antennas  considered: 

1.  Better  directivity. 

2.  Less  Attenuation. 

3.  Hobility  on  the  iaage  plane. 

4.  A  readily  variable  characteristic  inpedance. 

5. "  Relatively  frequency  independant  over  a  wideband. 

The  TEH  Horn  used  in  the  TDRL  is  physically  an 
open  sided  pyrimidal  structure  when  deployed.  It  is  a 
single  piece  of  flat  copper,  about  Ian  thick  and  triangular 
in  shape.  Figure  2.6  gives  the  relevant  dimensions.  The 
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Figure  2.6.  TEB  Born  Dimensions 

hola  at  the  apex  of  the  antenna  allows  for  connecting  to  the 
receiving  cable. 

The  receiving  horn  antenna  is  mounted  in  a  three 
dimensional  trapezoidal  support  mount  with  quarter-inch 
plexiglass  walls  on  the  top,  bottom,  and  sides.  The  front 
and  rear  are  open  to  allow  easy  insertion  of  the  copper 
plate.  The  dimensions  and  shape  of  the  mount  are  given  in 
Figure  2.7  ill  dimensions  are  outside  wall  measurements. 

The  flat  copper  antenna  is  slipped  into  the 
front  of  the  plexiglass  mount  and  protrudes  through  the 
back.  &  coaxial  cable  is  connected  to  the  back  with  an  rf 
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433  an- 


C.  FRONT 


Pigure  2.7.  TEH  Horn  Sapport  Bount  Dimensions 


end-connector  and  hand-tightaned  screw.  See  Figure  2.13  for 
a  diagra.  of  the  asseably  of  the  connection  point.  The 
antenna  is  physically  supported  within  the  aount  by  grooves 


cut  into  and  running  along  the  length  of  each  side.  The 
grooves  serve  to  fix  the  flare  angle  of  the  antenna  at 
approximately  24.1®. 

The  dimensions  of  the  aount  were  chosen  for 
reasons  of  practicality,  being  the  aost  appropriate  for  the 
aaterial  available  and  providing  a  reasonable  flare  angle 
and  antenna  characteristic  impedance  for  signal  reception. 
The  aount  also  provides  stabilty  to  the  Horn  under  various 
weather  conditions.  Thus,  by  placing  the  horn  in  the 
plexiglass  support  aount,  physical  and  electrical 
characteristics  are  fixed,  assisting  in  correlating  run  to 
run  aeasuraaents .  The  PSM  Horn  as  it  would  be  mounted  on 
the  imaging  plane  is  drawn  in  Figure  2.8.  The  support  mount 
is  not  shown  for  clarity. 

The  addition  of  the  plexiglass  mount  should  have 
negligible  effect  on  the  received  signal.  Figure  2.9  shows 
the  response  of  the  supported  horn  antenna  to  a  pulse 
generated  by  the  source.  This  compares  favorably  with  the 

■v 

pulses  obtained  by  Hammond  without  the  plexiglass 
aount.  [Bef.  11] 

ks  noted  in  Figure  2.8,  the  electric  field  at 
the  aperture  of  the  horn,  assuming  a  uniform  TFM  wave  being 
excited  between  the  horn  and  the  ground  plane,  is  given  by: 
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(2.1) 

=  amplitude  of  the  incident  voltage  at 
tha  feed  point. 

=  height  of  aperture 

*  wave  number  =  2TT/\ 

»  distance  of  observer  from  the 
referance  point. 

It  can  than  be  easily  shown  that  the  receiving 
antenna  impulse  response  in  the  far  field  is:  [Bef.  12] 

Sr(t)  =  -  c  £q  WZL  -5  (t  -  rQ/c)  (2.2) 

where:  &  =  impulse 

*  load  impedance 

and  the  impulse  response  of  the  transmitting  antenna  is: 

t  60e  W  . 

S  (t)  =  -  — p2-  3  (t-rQ/c) 

where:  r  =  distance  from  the  aperture  to  an 
observer  on  the  image  plane  in  the 
boresight  direction. 

Thus  the  transmitting  impulse  response  is  simply 
the  time  derivative  of  the  receiving  impulse  response.  The 
impulse  response  is  a  delayed  impulse,  delayed  in  time  by 
r#/c,  and  multiplied  by  a  negative  constant  which  is 
dependant  on  the  width  of  the  antenna  aperture,  w,  and  the 
load  impedance  ZL. 


where:  V0 

h 

!c 
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TARGET:  RUN  DATE  DIET  ANT  TGT  REMARKS 

1  6-26-81  N/A  H8  ???  GENERATOR  MAUEFQRM 


MAXIMUM  PEAK  VALUE 

(VOLTS) 

8.2328 

MINIMUM  PEAK  UALUE 

(VOLTS) 

-8.0932 

RMS  VALUE  < VOLTS) 

9.8796 

MEAN  VALUE  (VOLTS) 

9.8998 

NUMBER  OF  WAVEFORMS  AUERAGEO  -  1  OPTIMIZATION  UALUE  >  9 


9.59  1.99  1.38  2.99  2.39  3.99  3.39  4.89  4.39  3.99 
Figure  2.9.  typical  TEH  Horn  Receded  Wavefora 


Generally, 

the 

following  limitations 

are 

applicable 

to  the 

r  eh 

Horn  during 

most 

TDRL 

aeasureaent 

s:  [Ref.  13  ] 

1.  The 

antenna  is 

not 

in  the  far 

field  of 

the 

transmitting  antenna. 

2.  The  flare  angle  of  the  TEM  Horn  is  not  wide  enough  to 
receive  energy  from  all  positions  along  the 
transaitting  antenna. 

3.  Wavelengths  at  high  frequencies  aay  not  propagate  in 
a  true  TEH  aode,  resulting  in  attenuation  and 
distortion. 
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The  net  result  is  a  deviation 


from  the 


ideal 


results  of  the  impulse  responses  as  stated  above.  When 
Figure  2.9  is  compared  with  Figure  2.4,  distortion  of  the 
pulse  as  received  by  the  antenna  is  noted.  However,  for 
TDRL  purposes,  the  reproduced  waveform  is  still  reasonably 
faithful  to  the  original  source  pulse. 

c.  Receiving  Antenna  --  aonopoles 

Monopole  antennas  are  also  used  on  the  TDRL 
scattering  range.  One  of  eleven  aonopoles,  ranging  in 
length  froa  85  millimeters  to  285  millimeters,  in  20 
ailliaeter  lengths,  can  be  chosen.  Each  monopole  is 
constructed  of  brass  rod  and  is  threaded  on  one  end  to  allow 
for  ease  of  connection  to  the  receiving  elements. 

The  monopolas  may  be  considered  center-fed 
dipoles  with  their  iaages  in  the  ground  plane.  In  effect, 
they  are  oscillating  electrical  dipoles  in  free  space. 
Oscillation  at  frequencies  related  to  the 
guarter-wavelengths  of  the  aonopoles  is  to  be  expected. 

Figures  2.10,  2.11  and  2.12  confirm  the  above. 
The  effects  of  radiation  damping,  dispersion  and  reflection 
due  to  the  impedance  discontinuity  of  the  antenna  are 
clearly  demonstrated.  From  Figure  2.4  the  pulse  width  of 
the  source  pulse  is  seen  to  be  about  1.2  nanoseconds.  This 

corresponds  to  a  wavelength  of  about  360  millimeters. 
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TARGET:  RUN  DATE  OIST  ANT  TGT  REMARKS 

4  6-29-81  1.27M  2  ???  ACQUIRE  83HH  ANTENNA  RESPONSE 

MAXIMUM  PEAK  UALUE  (UOLTS)  8.0577 

MINIMUM  PEAK  UALUE  (UOLTS)  -8.0782 

RMS  UALUE  < UOLTS)  0.8237 

MEAN  UALUE  < UOLTS)  0.0800 

NUMBER  OF  MAUEFORMS  AUERAGED  -  1  OPTIMIZATION  UALUE  •  0 


0.30  1.00  1.30  2.00  2.30  3.00  3.30  4.00  4.30  3.00 

Figure  2.10.  85m  Honopole  Pulse  Response 


Figure  2.10  is  the  short-pulse  response  of  the 
85mm  monopole.  This  antenna  is  almost  exactly  one 
quarter-wavelength  in  height  in  relation  to  the  source 
pulsewidth.  It  correspondingly  produces  a  nearly  pure 
sinusoidal  variation,  lamped  in  time,  as  predicted  by 
classical  antenna  theory. 

Figures  2.11  and  2.12  are  the  observed  responses 
for  the  185  millimeter  and  285  millimeter  monopoles.  As  can 
be  seen,  the  antennas  tend  to  resonate  at  their 
quarter-wavelength  periods  of  2.5  and  3.8  nanoseconds 
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TARGET:  RUN  DATE  DIST 
3  6-29-91  1.27M 

MAX  I  HUN  PEAK  UALUE  (UOLTS) 

H1NINUH  PEAK  UALUE  (UOLTS) 

RHS  UALUE  (UOLTS) 

NEAN  UALUE  (UOLTS) 

NUMBER  OF  WAUEFORMS  AUERAGED  ■  1 


ANT  TGT  REMARKS 

???  ACQUIRE  109MN  ANTENNA  RESPONSE 
9.8839 
-9.8907 
0.8328 
8.0888 

OPTIMIZATION  UALUE  •  8 
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respectively.  Although  the  amplitudes  of  these  two 
aonopoles  do  not  differ  greatly  (less  than  .08%  between 
their  maximum  peak  values),  the  longer  antenna  tends  to 
distort  the  incident  wave  to  a  greater  extent.  This  is  most 
readily  observed  when  comparing  the  root-mean-square 
voltages  of  the  two  waveforms.  The  185  millimeter  antenna 
provides  more  effective  power  to  the  receiving  circuits  than 
does  either  of  the  other  two  aonopoles.  The  285  millimeter 
antenna  dissipates  a  considerable  amount  of  energy  in  the 
long  period  oscillations;  the  35  millimeter  antenna  is  too 
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saall  to  acquire  the  full  power  of  the  transaitting  antenna. 
There  may  then  be  a  "best”  nonopole  length  for  a  particular 
source  pulse.  Shat  this  length  sight  be  and  whether  the 
aonopoles  are  selective  of  targets  according  to  their 
dimensions  will  be  discussed  more  fully  in  Chapter  7. 

TARGET:  RUN  DATE  D1ST  ANT  TGT  REMARKS 

2  6-29-01  1.27ft  12  fVT  ACQUIRE  283NH  ANTENNA  RESPONSE 

MAX I HUN  PEAK  UALUE  <U0LT3>  0.0061 

MINIMUM  PEAK  UALUE  (UOLTS)  -0.0934 

RHS  UALUE  (UOLTS)  0.0209 

MEAN  UALUE  (UOLTS)  0.0000 

NUMBER  OF  HAUEFORNS  AUERAGED  ■  1  OPTIMIZATION  UALUE  -  0 


Figure  2.12.  285aa  Honopole  Pulse  Response 


411  receiving  antennas  are  coupled  to  the  Signal 
Processing  Group  through  a  coaxial  end-coupling  fitted 
through  the  iaaging  plane.  The  connection  poir*  is 
diagraaaed  in  cut-away  in  Figure  2.13.  The  antenna  is 
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Figure  2.13.  Receiving  Antenna  Connection  Point 


coupled  as  close  to  the  ground  plane  as  possible.  This  is 
particularly  important  for  the  TEM  Horn  to  prevent 
undesirable  oscillations.  As  drawn,  the  threaded  portion  of 
the  lonopole  (or  hand  tightened  screw  of  the  TEH  Horn)  is 
joined  to  a  connector  pin  by  butting  the  surfaces  of  the 
antenna  and  pin  together.  A  coaxial  cable  is  then  connected 
to  the  end  connector  and  the  rf  signal  is  transferred  to  the 
receiving  circuits. 

This  arrangeaent  has  worked  well.  No 

reflections  or  excessive  attenuation  of  the  incosing  signal 
due  to  the  joint  were  observed.  Oscillations  on  the  TEH 
Horn  were  sinisal. 

Table  I  lists  the  antennas  currently  available 
to  the  TDRt,  user.  The  significant  dimensions  and  pertinent 
characteristics  are  reviewed, 
d.  Targets 

The  targets  used  on  the  TDRL  scattering  range 
are  of  various  dimensions.  Currently,  eleven  are  used  in 
imaging  experiments.  They  are  all  axisyametric,  closed 
surfaces  with  edges  aad  vertices.  The  targets  say  be  broken 
into  three  catergories  according  to  general  shape: 

1.  Half-cylinders 

2.  Half- spheres 


TABLE  I 


TDBL  Antennas 


ANTENNA 

NUMBER 

DESCRIPTION 

DIMENSIONS 

COMMENTS 

1 

TEN  Horn 

Pigure  2.6 

Broadband 

2 

Dipole 

85mm 

Oscillating  Dipole 

3 

Dipole 

1 05mm 

4 

Dipole 

125mm 

5 

Dipole 

145mm 

6 

Dipole 

165mm 

7 

Dipole 

185mm 

8 

Dipole 

205mm 

9 

Dipole 

225mm 

10 

Dipole 

245mm 

11 

Dipole 

265mm 

12 

Dipole 

285mm 

3.  Half- cones 

Table  IX  lists  the  targets  and  provides  general  information 
about  each. 

There  are  six  half-cylinders  used  as  targets. 
They  are  all  machined  from  aluminum  stock.  The  largest  is 
12-inches  long  with  a  3-inch  radius.  The  smallest  is 
3.06-inches  long  and  1.5-incnes  in  radius. 

A  single  half-cone  is  also  used  as  a  target. 
The  cone  is  machined  from  the  same  stock  as  the 


half-cylinders 


Four  spheres,  varying  in  radius  from  5-inches  to 
1.5-inches  are  the  four  remaining  targets.  The  spheres  are 
basically  styrofoam  balls  cut  in  half.  3ach  half  is  covered 
with  heavy-duty  aluainua  foil  to  provide  the  proper 
reflection  of  transient  SB. 

2ach  of  the  targets,  when  placed  on  the  ground 
plane,  appears  to  the  incident  pulse  to  be  airrored. 
Therefore,  the  receiving  antenna  sees  not  the  half-shapes, 
but  the  full  cylinder,  cone  cr  sphere. 

No  particular  problems  were  noted  when  using  the 
targets  except  for  the  need  to  anchor  the  spheres  with 
weights  during  inclement  weather.  Otherwise,  they  have  a 
tendency  to  move  on  the  plane  due  to  the  wind, 
e.  Image  Plane 

The  final  element  of  the  Image  Plane  sroup  is 
the  image  plane  itself.  It  is  a  sguare  structure  of  welded 
Reynolds  Type  IN-190  aluainua  sheet.  It  measures  about  11 
meters  on  a  side  for  an  area  of  about  121  sguare  aeters.  It 
is  four  millimeters  thick. 

Based  on  the  above,  the  range  clear  time,  which 
establishes  the  low  frequency  cut-off  of  the  measurement,  is 
greater  than  36  nanoseconds.  The  clear  time  establishes 
the  first  arrival  at  the  measurement  point  of  a  reflection 

froa  objects  surrounding  the  plane  and  the  end  of  the  plane. 
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TABLE  II 


TDRL  Targets 


TARGET 

TYPE 

RADIOS 

LENGTH 

HATERIAL 

1 

Cylinder 

3" 

12" 

Aluminum  Stock 

2 

Cylinder 

1.5" 

12" 

Aluminum  Stock 

3' 

Cylinder 

3" 

6" 

Aluminum  Stock 

4 

Cylinder 

1.5" 

6.06" 

Aluminum  Stock 

5 

Cylinder 

3" 

3.06" 

Aluminum  Stock 

6 

Cylinder 

1.5" 

3.06" 

Aluminum  Stock 

7 

Sphere 

5" 

Al  and  Styrofoam 

8 

Sphere 

4" 

Al  and  Styrofoam 

9 

Sphere 

3" 

Al  and  Styrofoam 

10 

Sphere 

1.5" 

Al  and  Styrofoam 

11 

Cone 

3.18” 

5.94" 

Aluminum  Stock 

The  main  purpose  of  the  image  plane  is  to 
isolate  the  targets  and  the  antennas  from  the  underlying 
structure  of  the  TDRL.  Instrumentation  cables  are 

introduced  from  under  the  ground  plane  and  are  thus  kept 
fro*  interfering  with  the  aeasurenent.  The  plane  also 
serves  as  the  primary  support  structure  for  both  the 
transmitting  and  receiving  antenna. 

Ho  significant  problems  were  encountered  with 
the  image  plane.  However,  soma  disadvantages  of  this  plane 
arrangement  are  that  only  objects  having  a  symmetry  plane 

c 


can  be  measured,  and  polarization  and  incident  angles  are 
limited.  3ut,  for  the  work  intended,  the  image  plane  serves 
its  purpose  well. 

3.  Sampling  geceive£  gnd  Signal  grocessing.  Ggoug 

The  Signal  Processing  Group  consists  of  three 
elements: 

1.  Digital  processing  Oscilloscope  (DPO) 

2.  Mini-Computer 

3.  Hard  Copier 

The  Signal  Processing  circuits  are  connected  from 
the  output  of  the  receiving  antenna  to  the  input  of  the  DPO 
via  2.58  meters  of  HG-8A/0  coaxial  cable.  The  cable 

provides  some  attenuation  and  delay  of  the  source  signal 
provided  by  the  impulse  generator.  Additionally,  a  trigger 
signal  is  sent  to  the  DPO  from  the  pulse  generator  via  .54 
meters  of  RG-58C/U  cable. 

The  lengths  of  the  signal  and  trigger  cable  are 
chosen  for  convenience  and  for  minimum  attenuation.  They 
may  be  varied  for  special  purpose  applications  such  as  to 

accommodate  a  different  trigger  source  by  providing  a 

necessary  delay.  Care  must  be  taken  not  to  introduce  too 
much  delay  or  the  received  signal  will  not  be  windowed  on 

the  CHT  display  of  the  DPO.  The  maximum  practical  length  of 
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cabling,  including  that  to  the  transiitting  antenna  from  the 
pulse  generator,  should  be  lass  than  23  meters  for  the 
current  equipment  configuration  of  the  TDBL. 

a.  Digital  Processing  Oscilloscope 

The  Tektronix  Digital  Pror-  ,.,ing  Oscilloscope  is 
a  computer  compatible,  analog  to  digital  system.  It  is 
comprised  of  two  elemental  units;  a  Tektronix  7704A  General 
Purpose  Oscilloscope  System,  and  a  P7001  Processor. 
Modification  of  the  7704A  Oscilloscope  by  the  addition  of 
various  plug-in  nodules  allows  for  tailoring  of  the  whole 
DPO  system  to  the  specific  needs  of  the  TDRL.  A  cursory 
description  of  the  various  elements  of  the  DPO,  as 
applicable  to  the  TDRL,  follows. 

The  7704 A  Oscilloscope  system  is  composed  of  two 
parts;  the  D7704  Display  Unit  and  the  A7704  acquisition 
unit.  The  Display  unit  provides  results  of  signal 
acquisition,  both  real  time  and  as  stored  in  the  P7001.  It 
does  this  by  a  visual  interface  with  the  DPO  in  the  fora  of 
an  Oscilloscope.  Additionally,  the  CRT  readout  is  carried 
over  into  the  DPO  to  provide  appropriate  scaling  on  the 
displays. 

The  A7704  acquisition  unit,  through  its  modular 
plug-ins,  provides  the  means  to  acquire  the  analog  signal. 
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The  plug-in  units  used  in  the  TDHL  are  the  S-6  Sampling  Head 
and  the  S-53  Trigger  Recognizer,  inserted  into  a  7S12 
TDR/Saapler. 

The  7S12  provides  for  measurement  of  recurring 
fast-rise  tiae  signals;  i.e.,  the  source  pulse  and  the 
target  back-scatter.  It  features  the  aeans  of  determining 
vertical  deflection  factors  in  millivolts;  horizontal  sweep 
factors  in  seconds;  a  HIGH  RESOLUTION  switch  which  reduces 
the  waveform  noise  and  jitter  by  signal  averaging  (10 
averages  per  sample  point);  a  TIME-DISTANCE  scale  which 
allows  one-way  distance  measurements  in  air  dielectrics;  and 
a  locate  switch  which  increases  the  time/division  and 
intensifies  a  portion  of  the  display  to  locate  the  tiae 
window  relative  to  the  total  waveform.  Finally,  the  7S12 
provides  the  means  for  the  inserting  and  interfacing  of  the 
plug-in  heads  into  the  DPO. 

The  S-6  Saapling  Head  is  a  50-oha  loop  through 
input  saapling  unit.  It  provides  the  means  of  inputting  the 
signal  of  interest  to  the  DP3 .  This  unit  has  a  specified 
rise  tiae  of  30  picoseconds  or  less  and  a  bandwidth 
equivalent  to  dc  to  11.5  GHz  at  3  d3  down.  It  has  a  signal 
range  of  «■  Ivolt  to  -Ivolt  (dc  to  ac  peak)  or  Ivolt 
peak-to-peak .  The  S-6  establishes  the  baseline  measuring 
capabilities  of  the  TDRL. 
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The  S-53  Trigger  Recognizer  permits  use  of  the 
7s 1 2  as  a  general  purpose  sampler.  It  produces  a  stable 
trigger  signal  from  input  signals  from  DC  to  1  GHz.  It  is 
capable  of  recognizing  signals  having  amplitudes  ranging 
from  10  millivolts  to  2  volts  peak-to-peak  into  50  ohms.  A 
trigger  pulse  from  the  pulse  generator  is  fed  into  this  unit 
at  the  same  rate  as  the  source  signal.  Trigger  to  signal 
delay  is  specified  as  15  nanoseconds  or  less. 

The  second  element  of  the  DPO,  the  P7001 
Processor,  provides  the  link  from  the  DPO  to  the 
mini-computer.  This  unit  receives  the  analog  signal 
acquired  from  the  Acquisition  unit  and  digitizes  and  stores 
or  discards  the  information,  along  with  appropriate  scale 
factors,  as  selected  by  the  operator.  The  method  used  by 
the  digitizer  is  pseudo- random  sampling  with  one  sample 
taken  every  6.5  microseconds.  This  is  the  maximum  rate  of 
sampling.  Any  transient  longer  than  5  milliseconds  or 
COHERENTLY  REPETITIVE  signal  up  to  the  frequency  limit 
specified  by  the  Sampling  Head  that  can  be  displayed  on  the 
D7704  Display  Unit  can  be  stored  in  internal  memory  along 
with  its  scale  factors.  The  information  stored  can  then  be 
redisplayed  on  the  CRT  *  nd/or  sent  to  the  minicomputer  for 
processing.  There  are  a  maximum  of  512  samples  taken 
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horzontally  (9-bit  word)  per  waveform.  Although  either  X- 
or  T-  (TIME)  may  be  specified  for  this  axis,  the  TDRL 
primarily  employs  the  time  axis.  Vertical  resolution  of  the 
Processor  is  1024  levels  (10-bit  word) .  This  axis  is  used 
to  measure  the  voltage  of  the  back-scattered  signal. 

b.  Tektronix  4052  Mini-Computer 

The  minicomputer  used  in  the  TDRL  is  a  Tektronix 
4052.  It  is  a  high  performance,  integrated  graphics  system 
employing  LSI,  bi-polar,  16-bit  technology.  The  unit  has  a 
64k-byte  memory  with  a  300k-byte  magnetic  cartridge  built 
in.  Extended  BASIC  is  the  high  level  programming  language 
employed  to  perform  processing  functions.  Although  this 
unit  has  many  features  worthy  of  note,  the  above  are  the 
salient  ones  for  the  TDRL. 

c.  Hard  Copier 

The  Hard  Copier  used  is  the  Tektronix  4631  Hard 
Copy  Unit.  It  provides  permanent,  dry  copies  of  the  graphic 
and  alphanumeric  information  displayed  on  the  CRT  storage 
screen  of  the  4052  mini-computer. 

All  units  described  for  the  Signal  Processing 
Group  were  especially  selected  to  perform  as  an  integrated 
system  requiring  minimum  interfacing.  The  GPIS  bus  is  the 
means  of  transferring  information  from  one  unit  to  the  other 
in  the  TDRL. 
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III.  THSOBT 


This  chapter  is  devoted  to  the  general  derivation  of  the 
equations  which  are  applicable  to  the  tiae-domain  analysis 
of  the  impulse  and  tamp  responses  of  axisymmetric,  metallic 
bodies  which  are  illuminated  by  an  axially  directed  incident 
electromagnetic  field.  Where  sources  providing  detailed 
derivations  are  numerous  and/or  adequate,  the  derivations  in 
this  chapter  will  be  general,  providing  only  the  major 
results  of  interest  with  specific  constraints  noted. 

The  time-domain  integral  equations  which  provide 
solutions  for  the  transient  electromagnetic  problem  will  be 
considered  first.  They  will  be  followed  by  a  discussion  of 
the  measurement  equations  and  procedures  employed  in  the 
TDEL.  Pinally,  various  methods  used  to  reduce  the  noise 
error  in  acquired  signals  will  be  discussed. 

A.  TIHE-DOBAIN  INTEGRAL  EQUATIONS 

One  primary  problem  in  transient  analysis  of  target 
impulse  responses  is  to  determine  the  shape  of  an  unknown 
object  when  the  incident  field  and  the  response  of  the 
3catterer  to  that  field  are  a  priori  information.  Numerous 
tiae-domain  techniques  have  been  developed  which  perform 
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this  function.  As  noted  in  Chapter  I,  the  earliest  method 
used  involved  the  physical  optics  response  of  a  body  to  an 
incident  electromagnetic  field.  The  surface  currents 
produced  an  approximate  impulse  response  that  was 
essentially  the  second  derivative  of  the  projected  area 
function  of  the  scatterer.  [Bef.  5]  However,  this  approach 
was  limited  in  its  accuracy  in  that  the  relationship  between 
the  impulse  response  and  the  derivatives  is  exact  only  at 
the  leading  edge  of  the  scattered  field  response,  a  single 
point  in  time.  Interactions  between  currants  on  the  target 
body  which  continue  to  radiate  for  a  significant  time  after 
the  incident  field  leading  edge  has  propagated  further  in 
space  alters  the  bacfcscattered  response  of  the  target. 
Therefore,  for  an  accurate  solution  to  the  problem,  these 
additional  currents  must  be  of  consideration. 

Bennett  (Bef.  IB]  first  proposed  a  direct  method  which 
"corrected"  the  response  by  computing  the  currents  flowing 
on  the  scatterer  surface.  ?rom  these  currents,  the 
scattered  field  could  be  calculated  "exactly". 

In  general,  the  remainder  of  this  section  follows  the 
format  of  the  method  used  by  Bennett  in  the  derivation  of 
integral  equations  for  exact  solutions. 
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1.  *5221  lalegtjL.  Solution 

The  expression  for  the  magnetic  field  H  at  an 
arbitrary  point  in  space  (not  on  the  scatterer  surface)  is 
given  by: 


H(r,t)  -  Hl (r ,t)  *  ^  -fs  *  jj  J  C’1  . x  aR 

K 


(3.1) 


where:  H(r,t)  *  total  magnetic  field  at  (r,t) 

(?,  t)  =  incident  magnetic  field  at  (r,t) 
J(r',t)  =  surface  current  at  (£•,■£) 

r  =  position  vector  to  the 
observation  point 

a  position  vector  to  the 
integration  point 
a  =  I r-r ’  | 

A  T  -  f‘ 

=  ft 

t  »  tiae  in  light-meters  (one 

light-meter  is  the  tiae  it  takes 
light  to  travel  one  meter) 

T  =  t-R  =  retarded  tiae 

By  specializing  the  arbitrary  space  point  to  "r,  a 
point  on  the  surface  of  the  scatterer,  and  then  applying  the 
boundary  conditions  (i.e.,  by  causuality,  incident  field  is 
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zero  prior  to  arrival  of  the  incident  pulse  at  a  sample 
point) ,  an  integral  equation  for  the  current  density  on  the 
surface  of  the  scatterer  can  be  found: 


2an  x  H1(r,t)  ♦ 


4  an  x  {[^7  *  i.  3T1  ^(r'T) 


x  4r)  ds’  (3.2) 

where:  "r  =  position  vector  to  a  point  located  on 

the  surface  of  the  scatterer 
*  unit  vector  normal  to  the  surface 
The  first  term  on  the  right-hand  side  of  equation 
(3.2)  is  the  source  tera.  It  represents  the  direct 
influence  of  the  incident  field  on  the  current  at  the 
observation  point  (i^,t) .  It  is  the  physical  optics 
approximation  for  the  surface  current  when  applied  to  the 
illuminated  side  of  the  scatterer.  The  integral  on  the 
right  side  of  equation  (3.2)  represents  the  influence  of 
currents  at  other  surface  points  on  the  current  at  (?,t). 
Note  that  the  influence  of  other  currents  on  the  current  at 
(*r,t)  is  delayed  by  R.  This  allows  the  surface  current 
density  equation  to  be  solved  by  an  iterative  numerical 
procedure  in  the  time-domain,  rather  than  the  familiar 
frequency-domain  matrix  inversion  process. 

Once  the  surface  currents  have  been  found,  the 
far-scattered  field  can  be  calculated  by  using: 
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(r' ,t)  x  ar> 


(5.5) 


wherei  re  =  distance  to  the  far-field  observer 

£r  *  unit  vector  from  the  integration  point 
to  the  far-field  observer 

With  the  substitution  of  the  surface  current 
expression  into  the  above  equation,  and  the  assumption  that 
the  incident  field  is  a  rasp,  the  result  for  the  backscatter 
direction  is: 


ToHR(?’t)  *  27  S(V*H  *  57  Tc  fs 


t(r',T)  x  af)  ds 


(3.4) 


where:  roH^  -  backscatter  raap  response  of  the 
target 

s(ts>  *  physical  optics  silhouette  area  of 
the  target 

r  a  distance  of  the  far-field  observer 

O 

fron  the  origin 


t  +  r 
$  o 


Jft ,  correction  currents  resulting 


from  the  incident  raap  waveform. 
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This  result  provides  the  exact  relationship  between 
the  target  response  and  target  geometry.  A  particular 
unknown  sample  value  at  a  given  point  in  space  and  time  is 
determined  by  the  exciting  field  at  that  same  space-time 
point  and  by  the  scattered  fields  from  earlier,  more  distant 
locations.  note  that  the  interactions  between  unknown 
samples  are  displaced  in  time  by  an  amount  equal  to  the  time 
required  for  a  field  to  propagate  between  the  samples  at  the 
speed  of  light.  The  unknown  samples  can  be  solved  at  any 
time  step,  provided  all  sample  values  at  earlier  times  are 
already  known.  Thus,  by  determining  the  unknown  samples, 
i.e.,  the  correction  currents  ,  and  adding  the 

contribution  of  the  target  area  function,  s(t^),  the  target 
ramp  response  can  be  utilized  to  recover  information 
olving  the  unknown  dimensions. 

2.  I&2££S£  Ssatleiiai  Solution 

The  targets  employed  in  the  TDRL  are  all 

rotationally  symmetric  scatterers,  similar  to  the  one 

diagramed  in  Pigure  3.1.  The  general  constraints  applicable 
to  the  problem  are  that  the  scatterer  is  symmetric  about  the 
z-axis,  the  incident  field  is  axially  incident,  and  the 
far-field  is  computed  in  the  backscatter  direction. 
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The  contour  of  a  rotationally  symmetric  target  is 
completely  described  by  the  radius  vector  p  (z)  which 
varies  as  a  function  of  2.  The  projected  area  function  in 
the  plane  orthogonal  to  the  direction  of  the  propagation  of 
the  field  (z-axis)  can  be  expressed  simply  as: 

S  *  ffp2(2) 

Substituting  the  projected  area  function  into  the 
equation  for  the  backscatter  ramp  response  of  the  target  and 
solving  for  (z)  ,  an  inversion  equation  for  the 
rotationallly  symmetric  case  is  obtained: 


p(z) 


t)  -  27  3 V  £  iJCR(T',T)  x  ar}  t3,S^ 


1  3 


1/2 


Thus,  p  (z)  is  given  in  terms  of  the  measured  ramp  response 
(a  field  measurement) ,  and  in  terms  of  correction  currents 
at  earlier  times,  which  have  been  previously  computed  or  are 
known  to  be  zero. 

In  the  TDRL,  an  iterative  approach  is  employed  using 
estimates  of  the  entire  target  geometry  to  solve  for  the 
contour  function.  The  procedure  involves  five  steps: 

1.  Neglecting  the  surface  integral  providing  the 
correction  currents,  obtain  the  physical  optics 
estimate  of  0  (z). 
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{  2.  Determine  the  correction  currents  using  the  numerical 

solution  of  aquation  (3.2)  for  surface  current 
density. 

3.  Apply  the  correction  currents  to  equation  (3.5)  to 
obtain  a  next  estimate  of  (2)  • 

pl(2)  “  f2r0  H*|  (r,t)]1/2 

4.  Compare  the  nee  value  ^  (z)  with  (z)  to  see  if 
the  change  is  ] ess  than  some  small  number,  an  error 
factor.  The  normalized  mean  square  error  is  defined 
as: 

,2 .  <"»<«>  • 

(Hjjct))2 

where:  »  the  value  to  be  minimized 

in  order  to  obtain  the 
most  accurate  shape 
description 

I  S 

;  H  (t)  *  the  kth  estimation  of  the 

H, 

U-l 

backscattered  field. 

Thus,  the  difference  between  the  estimated  and  actual 
backscattered  fialds  are  directly  linked  to  the 

c 
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difference  between  successive  contour  functions,  and, 
therefore,  successive  shapes  of  the  body. 

5.  Hepeat  the  iteration  process  until  a  minimum  error  is 
obtained. 

The  above  procedure  was  used  successfully  by  florag 
[Ref.  15]  in  the  development  of  an  algorithm  involving 
theoretical  analysis  of  targets.  This  algorithm  has  been 
integrated  with  the  overall  TDRL  algorithm  for  inverse 
scattering. 

B.  BACKSC ATTERED  IATEPORB  HEASOEEBBBTS 

The  TDRL  algorithm  is  designed  to  minimize  the  effects 
of  the  impulse  responses  of  the  transmitting  and  receiving 
antennas.  Figure  3.2  diagrams  the  parameters  involved. 

As  can  be  seen,  there  are  two  measurement  paths  of 
interest.  The  first  is  the  Direct  Path  between  the 
transmitting  and  receiving  antenna.  The  only  factors 
affecting  the  transmitte  waveform  are  the  intervening 
distance  and  the  impulse  responses  of  the  receiving  and 
transmitting  antennas  to  the  source  pulse.  The  second  path 
is  the  Augmented  Path,  from  the  transmitting  antenna  to  the 
target  and  back  to  receiving  antenna.  Factors  effecting  the 
transmitted  pulse  in  this  path  include  the  antenna  impulse 
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responses,  the  impulse  response  of  the  target  and  the 

corresponding  intervening  distances. 

Three  waveforms  are  measured:  the  Direct  Waveform,  the 

Incident  Waveform  and  the  Augmented  Waveform.  During  the 
|  2 

measurements,  h  (t)  and  h  (t)  ,  responses  due  to  the  delays 
of  the  intervening  distances,  are  windowed  out  by  using  the 
time-range  gating  feature  of  the  DPO. 

The  Direct  Waveform  is  measured  first.  The  transmitting 
and  receiving  antennas  are  boresighted.  No  target  is  on  the 
imaging  plane.  The  measured  wave  is  then  the  convolved 
responses  of  the  source,  and  transmitting  and  receiving 
antenna  impulse  responses.  By  using  the  properties  of 

convolution,  the  frequency  domain  representation  converts 
the  convolution  process  into  a  product.  The  frequency 

domain  representation  cf  the  Direct  Waveform  is: 

A(ejw)  -  S(ejw)Ht(ejh,)Hr(ej‘‘':) 

where:  A(eiU)  *  Direct  Waveform 

S  (e'w  )  *  source  pulse 

H^e*w)  =  transmitting  antenna  response 
H*(ei‘'>)  »  receiving  antenna  response 

The  effects  of  the  intervening  distance  have  been  time 
gaxed  out.  Figure  3. 3  is  a  representative  Direct  Waveform 
measured  in  the  TDRL. 
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target:  run  date  oist  ant 

1  8-11-91  1.27M  1 

MAXIMUM  PEAK  UALUE . 

MINIMUM  PEAK  UALUE . . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OF  WAUEFORMS  AUERAGED  -  21 


TGT  remarks 

1  CYLINDER  ACQUISITION  AUERACE 
♦289.21  nU 
-96. S3  nU 
♦99.58  nU 
0.99  nU 

OPTIMIZATION  UALUE  «  9.5 


Figure  3.3.  Typical  Direct  JTatrefora — Cylinder 
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The  Incident  Waveform  is  measured  next.  The  imaging 
plana  is  as  set-up  for  the  Direct  Waveform,  but  the 
time-gating  is  further  delayed  by  JAt  to  eliminate  the 
effects  of  the  distance  intervening  between  the  receiving 
antenna  and  the  target  site.  The  resulting  measured 

waveform  is  the  Direct  Waveform  delayed  by  dUlt: 

B(e^“)  =  A(e^)  e(jwdAt) 

where:  B(e3‘',)  =  Incident  Waveform 
Figure  3.4  provides  a  representative  time-domain  Incident 
Waveform  measured  in  the  TDRi. 

The  final  waveform  measured  is  the  Augmented  Waveform. 
It  is  a  convolution  of  the  Direct  Waveform  with  the  impulse 
response  of  the  target  and  summed  with  the  incident 

Waveform.  For  this  measurement,  the  transmitting  and 
receiving  antenna  remain  boresighted.  A  target  is  placed  on 
the  image  plane.  The  measured  frequency  domain 

representation  is: 

C£eJ“)  -  [B(eJ“)  ♦  A(e^)  HS (e^)  ] eJ“dAt  C3,6) 

where:  C(ew)  =  Augmented  Waveform 
H  (eiw)  =  Target  Response 

Figure  3.5  provides  a  representative  time-domain  Augmented 
Waveform  measured  in  the  TD3L  for  a  axisymmetric 
half- cylinder. 
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TARGET:  RUN  DATE  D1ST  ANT  TGT  REMARKS 

1  8-11-81  1.27M  1  1  CYLINDER  ACQUISITION  AVERAGE 


MAX  I  HUM  PEAK  UALUE .  *4.48  ftU 

MINIMUM  PEAK  VALUE .  -3.83  nO 

RMS  UALUE .  +1.92  nU 

MEAN  UALUE .  8.80  nU 


OPTIMIZATION  UALUE  •  8.5 


NUMBER  OF  HAUEFORNS  AVERAGED  -  21 


HUMBER  OF  WAUEFORMS  AUER AGED  »  21 


As  can  be  seen  from  aquation  (3.6),  the  target  response 
is  not  independant  of  the  effects  of  the  impulse  responses 
of  the  antennas.  To  eliminate  this  dependance,  the 
backscattered  waveform  is  determined  and  then  the  remaining 
antenna  terms  are  theoretically  eliminated  by  a  division 
process. 

The  Backscattered  Waveform  is  the  difference  of  the 
Augmented  and  the  Incident  Waveforms: 

D(oj(U)  «  C(ej“)  - 

-  A(eju)  HS  (ejw) 

Figure  3.6  is  the  time-domain  Backscattered  waveform 
resulting  from  the  subtraction  of  the  waveforms  in  Pigures 
3.5  and  3.4  respectively. 

By  next  dividing  the  Backscattered  Wave  by  the  Direct 
Wave,  all  antenna  influences  are  eliminated: 

HSfejw)  -  Ellil 
A(eJ“) 

.  5 

Tae  target  impulse  response,  h  (t)  ,  may  now  be  determined  by 
deconvolution  of  the  target  impulse  frequency  domain  values. 
Figure  3.7  is  a  representative  time-domain  target  impulse 
response  for  a  cylinder.  Noise  errors  have  been  minimized 
by  an  adaptive  filtering  method  to  be  described  shortly. 
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target:  run  pate  dist  ant  tgt  remarks 

1  8-11-31  1.27M  1  1  CYLINDER  ACQUISITION  AUERAGE 


MAXIMUM  PEAK  UALUE .  ♦4.63  nU 

MINIMUM  PEAK  UALUE .  -4.44  nU 

RHS  UALUE .  +1.71  nU 

MEAN  UALUE . 8.80  nU 


OPTIMIZATION  UALUE  -  0.3 


NUMBER  OF  UAUEFORMS  AUERAGEO  «  21 


/ 


target:  run  date 

i  a-u-ai 

MAXIMUM  PEAK  UALUE. 
MINIMUM  PEAK  UALUE. 
RMS  UALUE . 


PtST  ANT  T GT 

1.27M  1  I 


TYPICAL  CYLINDER  RESPONSES 


0.0084  nU 


-0.0004  nU 


0.0002  nU 


MEAN  UALUE . . . 

NUMBER  OF  UAUEFORMS  AUERACED  •  1 


0.0000  nU 

OPTIMIZATION  UALUE  •  10 


nU 

0.0004 

0. 0002 

0.0000 
-0.0002 

0.S0  1.00  1. SO  2.00  2.30  3.00  3.30  4.00  4.30  3.00  nS 


IMPULSE  RESPONSE 


Figure  3.7.  farget^lipaise  ae*P°nS8  Soise 
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C.  NOISE  CO VSIDEBATON S 


As  noted  in  the  previous  section,  deconvolution  of  the 
target  bacicscattered  impulse  response,  less  the  impulse 
responses  of  the  antennas,  is  a  central  operation  to  the 
effective  use  of  the  r DRL.  Ihe  reliability  of  the  acquired 
signal  to  be  a  true  response  of  a  target  to  an  incident 
field  bears  directly  on  the  accuracy  of  the  final  results  as 
determined  by  solution  of  the  inverse  integral  equations  for 
time-domain  analysis.  Any  errors  in  the  data  representing 
the  Incident  Waveform  and  Augmented  Waveform  degrade  the 
cancellation  of  the  division  process  and  introduce  errors  in 
the  bacicscattered  target  response.  Also  present  at  the 
input  to  the  DPO  is  noise  due  to  thermal  electron  excitation 
in  the  transmission  lines  and  antenna  surfaces  as  well  a 
that  due  to  external  atmospheric,  cosmic  and  man-made  noise 
sources  being  received  by  the  horn  antenna.  When  the 
bacicscattered  impulse  response  is  transformed  to  the 
time-domain,  high  frequency  noise  in  most  cases  forms  the 
dominant  feature  of  the  waveform.  An  example  of  a  typical 
impulse  response  of  a  cylinder  is  given  in  Figure  3.8. 
Noise  errors  have  not  been  removed.  Figure  3.9  is  the  ramp 
response  for  the  same  target. 
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target:  run  date  dist  ant  tgt  remarks 

1  3-11-81  1.27M  1  1  TYPICAL  CYLINDER  RESPONSES 

MAXIMUM  PEAK  VALUE .  8.0324  nU 

MINIMUM  PEAK  VALUE .  -8.0314  nO 

RMS  VALUE .  8.8128  nU 

MEAN  VALUE .  8.0888  nV 

NUMBER  OF  WAVEFORMS  AVERAGED  ■  1  OPTIMIZATION  UALUE  >  8 


IMPULSE  PESPONSE 


0.8288 


0.8888 


-0.8288 


irfftfiTiiiil  1111 


pm  ti> 


8.38  1.00  1.38  2.88  2.38  3.88  3.38  4.88  4.38  3.88  nS 


Figure  3.8.  lapulse  Response— >  Noisy  Signal 
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TARGET:  RUN  DATE  OIST  ANT 
1  3-11-31  1.27N  1 

MAXIMUM  PEAK  VALUE . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OF  NAVEFORHS  AVERAGED  «  1 


TGT  REMARKS 

1  TYPICAL  CYLINDER  RESPONSES 
♦53. EE  nU-N 
-733.00  wU-M 
♦373.03  nU-H 
-235.42  mU-M 

OPTIMIZATION  UALUE  -  0 


RAMP  RESPONSE  —  TIME 


Figure  3.9.  Raap  Response— loisy  Signal 
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*0+  'WMMK. 


In  the  TDBL,  typical  sources  of  errors  in  the  Incident 
and  Augmented  Jfaveforas  are  sampling  errors  and  noise  in 
signal  acguisition,  leakage  and  aliasing  errors  in  the 
Fourier  Transf omations ,  and  rounding  errors  in  the  computer. 
The  error  of  primary  concern  is  that  involving  noise  in 
signal  acguisition  and  in  particular  the  enhanced 
high-frequency  noise  resulting  from  deconvolution  of  the 
antenna  responses  to  obtain  the  impulse  and  ramp  responses 
of  the  target.  This  noise  has  been  directly  traced  to  the 
structure  of  the  data  within  the  input  array,  as  will  be 
noted  in  Chapter  IV. 

Two  methods  are  employed  to  reduce  noise  to  acceptable 
levels:  acquisition  averaging  and  an  optimal  compensation 
technique.  [  Hef .  16] 

i •  Acquisition  hzsL aaiaa 

Acquisition  averaging  is  aimed  at  reducing  signal 
acquisition  errors  and  noise  by  averaging  several 
acquisitions  taken  within  a  short  period  of  time.  During 
the  process,  a  more  serious  error  may  be  generated,  that  due 
to  signal  delft  with  time.  Signal  averaging  has  been  most 
useful  in  ensuring  a  good  target  ramp  response.  It  has  had 
noticeable  effect  on  making  the  information  content  of  the 
impulse  response  more  available.  Figure  3.10  demonstrates 
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an  iapulse  waveform  with  210  acquisitions  for  the  saae 
target  as  for  Figure  3.8.  rheor et ically,  the  RMS  noise 
level  should  be  reduced  b7  the  square  root  of  the  number  of 
acquisitions,  N .  M  *  10  in  Figure  3.8,  and  N  =  210  in 
Figure  3.10.  The  noise  reduction  should  be  about  4.5  times. 
Mote  that  the  noise  level  in  Figure  3.10  is  indeed  about  one 
fourth  that  in  3.8. 

Figure  3.11  is  the  corresponding  ramp  response. 
Signal  acquisition  only  is  employed. 

2.  OatiJll  £2®2§asati2£  £a s^nigue 

The  optimal  compensation  technique  involves  the 
design  of  a  compensator  (deconvolution)  function  operating 
on  the  convolution  output  of  the  Incident  and  Backscattered 
Waveforms,  to  produce  a  relatively  noise-free  estimated 
Iapulse  Hesponse  for  the  target.  The  compensator  is  applied 
to  the  convolution  product  of  the  target  iapulse  response 
and  the  Direct  Waveform  to  produce  the  required 
deconvolution  result.  The  design  as  used  in  the  TDRL 
involves  an  iteration  process  on  a  single  variable  with  a 
"man-in-the-loop"  to  determine  optimal  values  of  the 
variable. 


Pigure  3.12  diagrams  the  design  of  a  frequency 
domain  optimal  compensator.  The  compensation  principal  in 
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TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

l  8-11-81  1.27M  1  1  CYLINDER  ACQUISITION  AVERAGE 

MAXIMUM  PEAK  VALUE .  0.0094  nU 

MINIMUM  PEAK  VALUE .  >0.0091  nV 

RMS  VALUE.... .  0.0033  nU 

MEAN  VALUE .  0.0000  nU 

NUMBER  OF  WAVEFORMS  AVERAGED  »  21  OPTIMIZATION  VALUE  ■  0 


nU 


0.0050 


0.0000 


-0.0050 


0.50  1.00  1.50  2.00  2.50  3.00  3.50  4.00  4.50  3.00  nfl 


Figure  3.10.  Impulse  Response:  Ho  Optimization— Icquistion 
Iveraged 
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TARGET:  RUN  DATE  DIST  ANT 
4  8-11-81  1.27M  1 

MAXIMUM  PEAK  UALUE . . . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OF  UAUEFORMS  AUER AGED  ■  21 


TGT  REMARKS 

1  CYLINDER-ACQUISITION  AUG 
♦32.97  aU-M 
-441.33  aU-H 
♦233.42  aU-M 
-178.4 6  aU-M 

OPTIMIZATION  UALUE  >  B 


B. 50  1.30  1.50  2.00  2.50  3.00  3.50  4.00  4.30  3.00  nS 


Figure  3.11.  Rasp  Response:  tfo  Optisizatioa--Acquistion  Av¬ 
eraged 
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Figure  3.12.  Frequency  Doa&in  Optiaal  Compensator  Design 

deconvolution  is  to  design  a  transfer  function  Q  (e-*"*)  to  be 
applied  to  the  backscatt  er  response  D(eJh*)  to  yield  an 
estimate  of  the  target  impulse  response  That  is: 

Hg(eiai)  -  D(e^“)Q(e^“)  (3.7a) 

—  -  A(ejw).HS(eju)Q(ej“)  (3.7b) 

The  main  design  criteria  is  to  minimize  the  energy 
in  the  noise  terms  that  introduce  error  in  the  estimate.  The 
error  energy,  B  is  given  by: 

Ee  “  4  I  He  (e^“)  |2  d<*»  (3.8] 

where:  /J,  *  frequency  band  of  interest. 

When  the  error  energy  is  minimized,  H^(ei,J)  will 

equal  H*(e>u)  and  the  impulse  responses  will  be  identical, 

except  for  the  addition  of  highly  enhanced  random  noise.  In 

this  case,  the  compensator  transfer  function  is: 
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H®  (e3“)  *  HS  (e3<J) 


(3.9) 


Q(ej“) 


A(e^) 


3u t ,  this  will  result  in  an  infinite,  unbounded  value  for 
Q(e3,,,)  as  A  (e‘,‘J)  goes  to  zero,  yielding  a  noisy  fori#  for 
HS(e'1'*).  This  noise  can  be  limited  by  constraining  Q(e  ) 
to  be  a  bounded  function.  Noting  that  the  impulse  response 
of  the  target  is  bounded  by  physical  constraints,  it  is 
possible  to  keep  the  energy  due  to  the  convolution  of  the 
impulse  response  and  the  designed  transfer  function  bounded, 
i.  e . : 

Ec  »  £  |  H(ej“)  Q(ei“)  |2  du»  (3-10) 

The  problem  now  is  to  minimize  2  while  keeping  E 
finite.  The  total  energy  in  the  backscattered  wave  can  be 
defined  as: 


E  *  Ec  +  XEC  ,  X  >  0 


(3.11a) 


?  [|  H?  (e>)  -  HS  (e3<u)  |2  ♦  X  iH(ej“)Q(ejtJ)  1 2]du>  (3.11b) 
0  © 


where:  X*  optimization  parameter. 
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Thus,  by  minimizing  the  total  error  energy  through 


the  selection  of  the  optimal  value  of  the  error  energy 
in  the  estimate  is  also  minimized.  The  result  yields  an 
optimum  compensator  {eJ‘*>|  for  the  unknown  target  impulse 
response  h  (t) .  The  result  is  that  the  whole  system  output 
is  low-pass  filtered  with  minimal  effect  on  waveform  shape. 

Adhering  to  the  design  constraints  listed  above,  it 
can  be  shown  that  the  form  of  the  compensating  transfer 
function  Q  (e1"* )  for  a  given  Direct  8aveform  A(eau>)  is: 
CHef.  17] 


Q(eju) 


[|A(eJSr*X) 


13.12) 


where:  A*(e1‘’*)  =  complex  conjugate 
Proa  equation  3.7,  the  compensation  deconvolution 
process  yields: 


S  „  D(ej“;A»(ejM)  C3.13) 

e  [ACeJ“)|-*\] 


Pigure  3.13  shows  the  effect  that  optimization  has 
on  the  impulse  response  with  no  sgnal  acquisition  averaging. 
Optimization  is  0.5.  The  target  impulse  response,  while 
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TARGET:  RUN  DATE  D1ST  ANT  TGT  REMARKS 

1  3-11-31  1.27M  1  1  TYPICAL  CYLINDER  RESPONSES 

MAX IMUN  PEAK  UALUE .  0.0011  ptU 

MINIMUM  PEAK  UALUE .  -0.0010  nU 

RMS  UALUE .  0.0004  «U 

MEAN  UALUE .  0.0000  mU 

NUMBER  OF  WAUEFORNS  AUERACED  -  1  OPTIMIZATION  UALUE  ■  0.5 


nU 

0.0010 

0.0005 

0.0000 

-0.0005 


IMPULSE  RESPONSE 


0.50  1.00  1.50  2.00  2.50  3.00  3.50  4.00  4.50  5.00  nS 


Figaro  3.13. 


Impulse,  Resp 
No  Acquisiti 


oase:  .5 
on  Arerag 


Optimization- 
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still  noisy,  is  clearly  discernable.  Figure  3.13  should  be 
compared  with  Figure  3.11.  Figure  3.14  shows  the  effect  of 
optimization  increased  to  10.  The  noise  is  almost  totally 
eliminated.  The  trade-off  is  that  the  optimization  has 
reduced  the  amplitude  of  the  response  and  spread  it  in 
frequency.  This  is  a  general  observation.  However,  it 
appears  that  the  degree  of  effect  on  the  waveform  is  a 
function  of  the  ideal  optimization  factor.  Chapter  V  will 
explore  this  concept  more  fully.  However,  it  should  be 
noted  that  if  the  optimization  factor  is  large  enough  to 
reduce  the  large  scale  features  of  the  target  then  it  is  too 
large  and  should  be  reduced— this  is  the  '•man  in  the  loop's" 
job. 

The  combination  of  waveform  averaging  and 
optimization  on  the  impulse  response  is  demonstrated  in 
Figures  3.15  and  3.16.  With  a  relatively  small  amount  of 
optimization,  the  target  impulse  response  is  readily 
apparent.  See  Figure  3.16.  Signal  acquisition  averaging 
also  yields  an  important  advantage.  Since  noise  is  reduced 
directly  by  the  square  root  of  the  number  of  acquisitions, 
and  because  no  filtering  of  any  sort  is  employed  for 
averaging  alone,  the  detailed  features  of  the  actual  signal 
are  not  eliminated  as  with  filtering.  The  smaller  amount  of 
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optimization  needed  results  in  minimum  effect  by  the 
compensator  on  the  shape  of  the  output  waveform.  The 
trade-off  is  the  increased  liklihood  of  error  introduced  due 
to  signal  time  drift  during  the  acquistion  averaging  period. 

Figures  3.17  to  3.21  show  the  effects  of  signal 
averaging  and  optimization  on  the  Ramp  Responses.  Each 
figure  is  presented  in  the  same  order  as  were  the  Impulse 
Responses.  Signal  averaging  is  very  beneficial  in  obtaining 
an  accurate  ramp  response.  Optimization  has  a  greater 
tendency  to  effect  the  wave  shape  of  the  ramp  response  than 
the  impulse  response  of  the  target.  But  less  optimization 
is  required  to  reduce  the  error  noise.  It  is  generally 
possible  to  use  the  ramp  response  without  optimization.  The 
main  effect  of  optimization  is  to  spread  the  frequency 
content  of  the  ramp  response.  amplitude  effects  are  minor. 
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target:  run  date  oist  ant  tgt  remarks 

1  8-11-01  1.27N  1  1  TYPICAL  cylinder  responses 

0.0004  nU 
-0.0004  nU 
0.0002  nU 
0.0000  nU 

OPTIMIZATION  UALUE  ■  10 


nU 


0.0004 


0.0002 

0.0000 


-0.0002 


0.SO  1.00  l.SO  2.00  2.50  3.00  3.30  4.00  4.30  3.00  a 


MAXIMUM  PEAK  UALUE . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OF  HAUEFORMS  AUERAGED  -  1 


Figure  3.10.  Impulse  Response:  10  Optimization — Mo  Acquisi¬ 
tion  Averaging 
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TARGET!  RUN  DATE  DIST  ANT  TGT  REMARKS 

1  8-11-91  1.27N  1  1  CYLINDER  ACQUISITION  AVERAGE 

MAX  I  NUN  PEAK  VALUE .  8.8994  nU 

MINIMUM  PEAK  VALUE .  -0.0891  nV 

RMS  VALUE .  0.0033  nU 

MEAN  VALUE .  0.0000  nU 

NUNBER  OF  WAVEFORMS  AVERAGED  «  21  OPTIMIZATION  VALUE  «  0 


nV 


IMPULSE  RESPONSE 


0.0030 

0.0000 

0.0030 


0.30  1.00  1.30  2.00  2.30  3.00  3.30  4.00  4.30  3.00  *S 


Fiaora  3.15.  Iapalse  Basponsa:  No  optimization — 210  icqui 
*  sizions  lFetagad 
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TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

1  8-11-81  1.27M  1  1  CYLINDER  ACQUISITION  AVERAGE 

MAXIMUM  PEAK  UALUE .  8.8818  nU 

MINIMUM  PEAK  UALUE .  -8.8889  nU 

RMS  UALUE .  8.8883  *U 

MEAN  UALUE . 8.8888  nU 

NUMBER  OF  HAUEFORHS  AVERAGED  «  21  OPTIMIZATION  UALUE  -  8.5 


hU 

8.8818 

8.8885 

8.8888 

-8.8885 


IMPULSE  RESPONSE 


Pigure  3.16.  Ispalse  Response:  .5  Optimization — 210  icgais- 
tions  Averaged 
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TARGET!  RUM  DATE  DIST  ANT  TGT  REMARKS 

1  0-11-81  1.27N  1  1  TYPICAL  CYLINDER  RESPONSES 

MAXIMUM  PEAK  VALUE .  +33.  66  nU-M 

MINIMUM  PEAK  UALUE .  -700.00  mU-N 

RMS  UALUE .  ♦373.03  nU-M 

MEAN  UALUE .  -283.42  mV- M 


NUMBER  OF  WAVEFORMS  AVERAGED  ■  1  OPTIMIZATION  VALUE  -_0 

RAMP  RESPONSE  —  TIME 


Figurs  3.17.  Bas£  Response:  So  Optiaisation--Ho  Requisition 


TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

1  a- 1 1-81  1.27M  1  1  TYPICAL  CYLINDER  RESPONSES 

MAX  I  NUN  PEAK  UALUE .  *4?.  77  nU-H 

MINIMUM  PEAK  UALUE .  -647.48  nU-H 

RMS  UALUE .  +367.11  nU-M 

MEAN  UALUE . . .  -279.88  nU-N 


NUM8ER  OF  UAUEFORMS  AUERAGED  »  1  OPTIMIZATION  UALUE  •  8.5 

RAMP  RESPONSE  —  TIME 


Figure  3.18.  Baap  Response:  .5  Optiaization— So  Acquisition 
Averaging 
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target:  run  date  dist  ant  tgt  remarks 

1  8-11-81  1.27N  1  1  CYLINDER  ACQUISITION  AUERAGE 

MAX  I  HUH  PEAK  UALUE. .  +69.81  «U-M 

MINIMUM  PEAK  UALUE .  -481.43  nU-M 

RMS  UALUE .  +231.11  nU-N 

MEAN  UALUE . -170.38  nU-M 

NUMBER  OP  UAUEFORMS  AUERAGED  -  1  OPTIMIZATION  UALUE  >  18 


RAMP  RESPONSE  --  TIME 


Rigors  3.19.  Rasp  Response:  10  Optisization— Ro  Acqoisitioa 
Averaging 


TARGET:  RUN  DATE  OIST  ANT  TGT  REMARKS 

t  8-11-31  1.27M  1  1  CYLINDER  ACQUISITION  AVERAGE 

MAXIMUM  PEAK  VALUE .  *78.58  nU-M 

MINIMUM  PEAK  VALUE .  -4S7.93  nV-M 

RMS  VALUE .  *289.5?  nV-N 

MEAN  VALUE .  -223.73  nV-M 


NUMBER  OF  WAVEFORMS  AVERAGED  *  21  OPTIMIZATION  VALUE  »  8 

RAMP  RESPONSE  --  TIME 


Figure  3.20.  gy/gB"!  ‘e’ttiSi- 
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TARGET:  RUN  PATE  01  ST  ANT 
l  9-11-81  1.27N  1 

MAXIMUM  PEAK  VALUE . . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OP  WAVEFORMS  AUER AGED  -  21 


TGT  REMARKS 

1  CYLINDER  ACQUISITION  AVERAGE 
♦72.23  hU-M 
-465.94  nU-M 
♦288.81  nU-N 
-222.20  nU-M 

OPTIMIZATION  UALUE  •  8.3 


*rtl-M 

8 

-200 

-400 

0.38  1.00  1.30  2.00  2.50  3.00  3.30  4.00  4.30  3.00  nS 


Figar*  3.21.  *5  °PtiBizatioa — 210  Acsaisi- 
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IV.  Ogsg AJING  SISISfl  C05E0IS1  IHPLSHBHTATION 


A.  I5TH0D0CTI0M 

A  set  of  software  programs  called  the  "Operating  System" 
controls  the  input,  processing  and  output  of  information 
obtained  in  the  TDRL.  It  is  the  quality  and  throughness  of 
the  Operating  System  that  directly  determines  the  quality  of 
the  results  obtained.  In  this  chapter,  there  will  be  a 
discussion  in  detail  of  the  implementation  of  the  Operating 
System.  It  will  begin  by  discussing  the  physical  structure 
of  the  the  host  processor.  A  brief  overview  of  the 
Operating  System  will  follow,  defining  the  basic 
sub-programs  that  make-up  the  total  algorithm.  Finally,  a 
detailed  description  of  each  program  will  be  provided, 
emphasizing  the  method  of  implementation  of  the  numerical 
solutions,  where  applicable,  and  giving  details  of  a 
procedure  to  run  tha  individual  programs. 

B.  THE  TBKTBONIX  4052  GRAPHIC  COHPOTING  SYSTEM 

The  host  processor  used  in  the  TDRL  is  the  Tektronics 
4052  Graphic  Computing  system.  It  is  a  versatile,  stand 
alone  microprocessor.  As  such,  its  capabilities  and 
limitations  are  important  driving  functions  that  shape  the 

structure  of  the  Operating  System. 
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The  4052  has  three  memories.  They  are: 

1.  Head  Only  Memory  (ROM) — a  permanent  memory  that 
contains  the  System's  intelligence. 

2.  Random  Access  Meaory  (RAM) — working  memory  of  the 
processor's  CPU . 

3.  Line  Buffer — small  temporary  memory  that  allows  data 
to  be  written  and  edited  on  a  display  screen  prior  to 
releasing  it  to  the  RAM.  An  important  editing 
feature. 

Built-in  peripherals  for  the  Graphic  System  (GS)  are  a 
smart  keyboard  (primary  input)  ,  a  CRT  visual  display  screen 
(primary  output)  and  a  magnetic  tape  unit  (mass  storage). 
Three  extenal  peripheral  devices  extend  the  versatility  of 
the  system.  They  are  a  Tektronix  4631  Hard  Copy  Onit  (makes 
paper  copies  of  display  information),  a  4924  Digital 
Cartridge  Tape  Drive  (additional  mass  storage)  and  a  Digital 
Processing  Oscilloscope  (data  acquisition).  A  block  diagram 
for  the  system  is  shown  in  Figure  4.1.  [Ref.  18] 

The  Processor  is  the  main  computing  device.  It 
corresponds  to  the  Central  Processing  Unit  (CPU)  found  in 
larger  systems.  It  maintains  the  "firmware"  that  allows  the 
GS  to  direct  system  operations,  decode  instructions  and 
perform  arithmetic  and  logic  operations. 
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Figure  4.1.  Block  Diagraa  of  TDRL  Graphic  Systea 

The  General  Purpose  Interface  provides  the  aeans  for  the 
systea  to  coaaunicate  inforaation  to  or  receive  data  froa 
any  external  peripheral  that  is  conpatible  with  I2EE 
Standard  488-1975.  This  interface  is  a  bit-parallel, 
byte-serial  systea  capable  of  transfering  ASCII  code  or 
aachine  dependant  binary  code  up  to  8-bits  vide  at  a  aaxiaua 


A  magnetic  tape  data  cartridge  is  used  in  the  internal 
magnetic  and  peripheral  tape  drive  units  to  provide  up  to 
300-k  bytes  of  additional  mass  storage  per  tape.  Programs 
and  data  may  be  stored  on  tape  in  specifically  defined 
"files".  Piles  nay  range  in  size  from  763-bytes  up  to  the 
limit  of  the  magnetic  tape  storage  capability.  Files  may  be 
either  ASCII  or  binary  in  format.  All  files  in  the 
Operating  System  are  stored  in  Binary  format  to  take 
advantage  of  a  peripheral  "firmware"  package  called  the 
binary  loader.  Binary  data  transfers  between  the  processor 
compiler  and  magnetic  tapes  are  faster  than  when  stored  in 
ASCII.  Thi3  is  because  the  processor  uses  binary  format  for 
internal  processing  and  the  conversion  from  ASCII  to  binary 
format  is  eliminated  when  data  is  stored  binary. 
Additionally,  a  stored  binary  program  can  be  transfered  to 
the  processor  memory  from  a  specified  peripheral  device 
without  disturbing  variables  and  associated  values 
previously  defined.  This  feature  is  not  available  to  ASCII 
programs. 

The  GS  permanent  ROB  has  been  programmed  to  respond  to 
BASIC  (Begianer*3  All-Purpose  Symbolic  instruction  Code),  a 
high-level  programming  language  first  developed  at  Dartmouth 
College.  The  Operating  System  has  b9en  correspondingly 
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written  in  this  language.  However,  BASIC  as  used  by  the  GS 
has  been  extended  to  fit  the  language  to  the  specialized 
capabilities  of  Graphic  Systaa.  These  extensions  are  in  the 
areas  of  graphics,  file  systea  access,  unified  handling  of 
input/output  operations,  matrices,  character  string 
manipulation,  high-level  language  interrupt  handling,  and 
operating  systea  facilities.  A  key  difference  between 
extended  GS  BASIC  and  most  other  BASIC  languages  is  that 
most  keywords  and  their  paraaeters  can  be  evaluated 
independently  of  program  control.  The  result  is  a  rich, 
versatile,  yet  simple  to  use  programming  tool. 

Three  peripheral  BOB  firmware  packages  have  been  added 
to  the  TDRL  system  to  tailor  it  to  specific  needs  of  the 
laboratory.  They  are  the  Signal  Processing  Ron  Pack  Mo.  1, 
Signal  Processing  ROM  Pack  Mo.  2  (FFT)  ,  and  the  EDITOR.  A 
description  of  each  is  provided  in  Appendix  A. 

C.  OY1RYIER  OP  THE  TDRL  OPRRATIHG  SYSTEM 

The  TDRL  Operating  System  has  been  specifically  designed 
to  be  higly  operator-processor  interactive.  This  provides 
the  greatest  degree  of  flexibility  in  the  acquisition, 
processing  and  evaluation  of  electromagnetic  time  domain 
transients  for  target  imaging. 
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The  total  Operating  System  occupies  33818-bytes  of  aass 
storage  on  TD8L  Library  Tape  1  in  three  files.  It  requires 
an  additional  201628-bytes  of  storage  in  62-files  on  TDRL 
Library  Tape  2  (Data  Storage) .  The  Operating  Systea  is 
logically  divided  into  three-individual  programs: 

1.  INPUT— for  the  acquisition  and  averaging  of  signals 
of  interest. 

2.  MATH--perforas  the  significant  aarhaatical  processing 
of  the  acquired  signals  for  inverse  scattering. 

3.  GRAPH— provides  a  visual  display  to  the  operator  of 
the  results  of  the  processing. 

TABLE  III 

Hass  storage  TDRL  Prograas 


PROGRAM 

FILE  * 

BITES 

% 

INPUT 

2 

24552 

12 

HATH 

3 

14256 

7 

GRAPH 

4 

42624 

22 

Hass  storage  requirements  of  the  Operating  Systea 
prograas  are  given  in  TABLE  III.  The  storage  requireaents 
listed  in  TABLE  III  are  for  the  prograas  only  and  do  not 
include  requireaents  for  the  necessary  definition  of  arrays* 
strings  or  nuaeric  constants  and  variables. 
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ill  programs  are  designed  to  be  as  nearly  independant  of 
each  other  as  possible,  after  an  initial  run,  to  allow  for 
maximum  flexibility  and  utility  in  data  usage  and 
evaluation.  The  aost  dependant  prograa  is  HATH  which  cannot 
be  used  without  inputs  froa  either  INPOT  or  GRAPH.  The  aost 
independant  prograa  is  SR&PH  which  is  capable  of  being  run 
with  inputs  froa  aass  storage  only. 

&  general  flow  diagraa  of  the  operating  systea  is  given 
in  figure  4.2.  Specific  naaed  subroutines  have  not  been 
shown.  The  next  section  will  describe  in  detail  the 
individual  prograas  of  the  Operating  Systea  and  their 
specific  reguireaents  and  processing  methods. 

0.  OPERATING  SISTER  PROGRAM  DESCRIPTION 

This  section  will  describe  in  detail  the  individual 
prograas  and  their  subroutines  that  aaJce-ap  the  Operating 
Systea  for  the  TDRL.  &  general  description  of  the  purpose 
of  each  prograa  will  be  provided  first.  Sources  of  inputs 
and  outputs  will  be  noted.  i  general  flow  diagram  for  the 
specific  prograa  will  be  provided  showing  how  all 
subroutines  interact.  This  general  overview  will  then  be 
followed  by  a  more  detailed  description  of  each  subroutine 
in  the  program.  The  purpose  of  the  subroutine  will  be 
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General  Flow  Diagram  of  TDRL  Operating  Systea 


noted.  Sources  of  input  and  output  will  be  defined.  & 
discussion  of  applicable  algorithms  and  the  flow  of 
information  will  be  also  provided.  Examples  of  visual 
displays  applicable,  potential  problems  and  correction 
methods,  will  be  discussed  when  warranted. 

As  noted,  the  following  commentary  will  make  frequent 
referance  to  program  generated  messages.  Examples  of  these 
messages  will  not  be  given  in  the  body  of  this  text,  but  can 
be  found  in  Appendix  F,  a  general  run  of  the  Operating 
System. 

1.  Program  INPUT 

INPUT  provides  the  means  to  initialize  and  drive  the 
total  Operating  System.  This  sub-program  is  stored  on 
Library  Tape  1,  File  2.  Memory  space  required  for  its 
implementation  in  the  processor  CPU  is  24552-bytes.  With 
all  arrays  and  constants  defined,  a  total  of  42469-bytes  of 
memory  space  may  be  required.  A  listing  of  INPUT  can  be 
reviewed  in  Appendix  C. 

INPUT  i3  the  front-end  program  for  the  Operating 
System  and  as  such,  will  generally  be  the  first  program  to 
be  initiated  for  TDRL  measurements.  Program  start-up  is 
begun  by  placing  Tape  1  into  the  4052  built-in  magnetic  tape 
drive,  and  placing  Data  Storage  Tape  2  in  the  4924  fiemote 
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Tape  Drive.  INPOT  is  then  normally  placed  into  the  CPU's 
working  memory  under  operator  control  by  the  issue  of  the 
following  commands  from  the  keyboard: 

FIND  2 
CALL  "BOLD" 

Program  execution  is  next  begun  by  typing  the  command  "RON" 
from  the  keyboard. 

If  this  is  the  first  attempt  to  run  the  Operating 
System  in  conjunction  with  the  Digital  Processing 
Oscilloscope  (it  is  assumed  the  DPO  is  "ON"  whenever  INPOT 
is  employed) ,  the  following  message  may  appear  on  the  GS  CRT 
following  the  issuance  of  "FIND  2": 

"Ho  SRQ  on  unit  in  immediate  line  —  message  number  43." 

If  the  program  is  being  run  from  any  line  prior  to  line  140 
(the  usual  case),  this  message  may  be  ignored  and  the 
program  will  automatically  correct  the  error  message  and 
continue  processing.  If  the  program  is  being  run  subsequent 
to  line  140  (the  case  when  returning  to  INPUT  from  some 
other  internal  file) ,  the  following  command  must  be  entered 
from  the  keyboard: 

POLL  H,fl;1 

Following  this,  the  regular  sequence  of  commands  may  then  be 
issued  to  start  the  system. 
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As  noted,  INPUT  may  alternately  be  loaded  and 
initiated  under  prograa  control  froa  GRAPH  (Pile  4) . 
Execution  will  begin  following  line  180  in  this  case.  This 
return  to  INPUT  occurs  in  normal  Operating  System  (OS) 
operation  after  an  initial  run  has  been  completed  and 
information  is  then  desired  for  an  entirely  new  target. 

The  program  specifically  provides  for  the  means  to 
input  information  concerning  the  transient  electromagnetic 
response  of  a  target.  It  acquires,  averages  and  stores  all 
pertinent  waveforms.  INPUT  also  allows  for  the  easy  storage 
of  specific  antenna  and  target  dimension  parameters  for 
comparison  when  the  target  is  known  exactly.  This  data  is 
permanently  maintained  in  a  library  on  Tape  2  in  Piles  5-49. 
INPUT  further  initializes  most  of  the  parameters  used  later 
in  Prograa  GRAPH  for  the  visual  display  of  the  processed 
results. 

Data  inputs  are  froa  the  keyboard  or  internal 
memory.  Outputs  are  to  internal  memory,  mass  storage  and 
Pile  3  (HATH) .  INPUT  is  highly  operator-processor 
interactive.  Numerous  directives  and  visual  displays  guide 
the  user. 

The  structure  of  the  prograa  consists  of 
1 1 -subroutines.  One  (DRIVER)  is  of  general  application  to 
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the  program.  Five  (ACAVG,  INWAV,  DIHAV,  AOGWAV,  SCATtf  AV) 
are  specifically  related  to  input,  acquisition,  averaging 
and  storage  of  relevant  information  concerning  the 
time-domain  transient  response  of  the  target  in  question. 
The  remaining  five  subroutines  (ANTLOC,  TGTLOC,  NEWPAR, 
ANTIS,  TGTIN)  provide  special  services  for  the  operating 
system.  A  general  flow  diagram  for  input  is  provided  in 
Pigure  4.3  A  description  of  each  subroutine  follows, 
a.  DRIVES 

This  is  the  main  driving  routine  for  INPOT.  It 
consists  of  129  steps  (40*  of  the  program)  from  line  100 
through  line  1260.  It  is  initiated  directly  under  operator 
control  at  line  100,  or  under  program  control  from  GRAPH  in 
line  180.  Data  is  input  from  the  keyboard  or  from  the  mass 
storage  files  located  on  Tape  2.  OBITER  will  pass  program 
control  directly  to  HATH  when  directed. 

Lines  100-300  initialize  various  variables  and 
dimension  various  arrays.  Line  150  clears  the  General 
Purpose  Interface  Bus  (GPIB>  of  the  DPO  service  request. 
Flag  ?7  (indicates  whether  the  acquired  data  is  destined  for 
use  in  Time  Domain  or  Prony  Net hods)  is  set  to  0  (Time 
Domain)  in  line  160.  ?8,  the  run  number,  is  initialized  to 

zero  in  line  170  and  will  be  incremented  by  one  in  line  210 
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Iron  file 
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Figure  4.3.  General  Floe  Diagraa  of  IHPUT 


for  every  run  (including  the  first).  Various  string 
constants  are  initialized  in  line  240-250  and  280-290.  Plag 
P9  (indicates  the  need  for  initializing  input  parameters)  is 
set  to  1  (input  paraieters  required) .  Lines  320-530  provide 
a  menu  of  subroutines  and  functions  which  are  available  for 
Operator  selection.  The  selection  of  a  process  from  the 
menu  is  input  from  the  keyboard  into  system  variable  P0.  An 
example  of  the  menu  is  given  in  Appendix  ?. 

Lines  550-570  direct  the  program  flow  to  the 
appropriate  portion  of  DRIVER.  If  the  operator  has  selected 
one  of  the  first  four  subroutines,  and  this  is  an  initial 
run  of  the  OS,  then  the  time  for  a  full  screen  sweep  of  the 
DPO  must  be  inputted  to  the  program  from  the  keyboard.  The 
sweep  rate  is  put  in  in  seconds  in  ?1  in  line  640  and  must 
be  greater  than  or  equal  to  the  time  required  for  the  the 
DPO  beam  to  sweep  completely  across  the  DPO  screen. 

Program  control  is  now  passed  to  line  350-1070 
where  the  following  run  descriptive  parameters  are  input: 

1.  Date 

2.  The  distance  the  target  is  from  the  transmitting 
antenna 

3.  The  antenna  number  selected 

4.  The  target  number  selected 
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5.  Explanatory  remarks 

These  are  all  keyboard  entries.  An  example  of  the  messages 
that  query  entry  of  the  above  parameters  is  given  in 
Appendix  F. 

If  the  run  is  not  the  first  for  the  Operating 
System,  a  different  series  of  inputs  are  queried.  In  this 
case,  program  control  is  passed  to  lines  670-840  which 
generate  messages  on  the  3S  asking  the  operator  if  there  is 
a  change  in  the  original  location  of  the  target,  if  a  new 
antenna  is  to  be  used,  or  if  a  defined  target  is  to  be 
changed.  The  operator  inputs  the  alpha  character  "Y"  or  "N" 
for  "yes"  or  "no"  from  the  keyboard  to  reflect  the  current 
situation.  Appendix  F  displays  these  messages. 

Lines  1000-1040  generate  the  input  of  the 
antenna  parameters.  These  parameters  are  drawn  from  files 
stored  on  Data  Tape  2.  Each  antenna  file,  located  in  one  of 
the  files  numbered  5-24,  is  a  63-element  array,  returned  in 
A5.  In  elements  1  through  60  are  located: 

1.  the  run  number  for  this  particular  antenna-target 
combination 

2.  the  antenna  length  compensating  parameter  for  the 
antenna-target  combination 
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3.  the  antenna  amplitude  compensator  for  the 
antenna-target  combination 

4s  can  be  seen,  antenna  parameters  are  actually  three 
elements  within  the  larger  array,  each  three  elements 
identifying  a  unique  antenna-target  combination.  The 
parameters  in  elements  61-63  are  system  antenna  parameters 
used  for  unknown  targets.  They  are  simply  the  sum  of  all 
run,  antenna  length  compensator  updates  and  antenna 
amplitude  compensator  updates.  The  average  of  the  last  two 
elements  is  applied  to  signals  acquired  if  the  target  is  not 
defined  from  the  target  library.  This  is  shown  in  lines 
102 4-1025. 

The  antenna  amplitude  compensator  is  a  very 
sensitive  function  of  the  distance  the  target  is  from  the 
transmitter  due  to  EH  propagation  attenuation.  In  line 
1025,  for  unknown  targets,  and  in  line  4152  for  defined 
targets,  these  losses  are  taken  in  consideration  by  finding 
the  square  of  the  ratio  of  the  targets  true  location  on  the 
plana  with  reference  to  the  transmitting  antenna  and  the 
control  distance  (1.27m). 

lfter  the  definition  of  variables  is  complete, 
DBI7EB  passes  control  to  the  appropriate  function  in  line 
1090.  Line  1100  returns  the  display  to  the  Henu  for  further 
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function  selection  or  to  line  1170  where  the  DRIVES  passes 
control  to  MATH  and  terminates  any  further  INPOT  processing. 

In  summary,  INPUT  DRIVER  requires  input  of 
selected  program  routines  and  initial  parameters  by  the 
operator  from  the  keyboard,  and  requires  the  OS  to  determine 
whether  the  current  run  is  the  initial  for  the  day  or  a 
subsequent  one.  The  parameters  are  saved  for  later  graph 
reguiremts.  Funtion  selection  drives  the  appropriate 
subroutines  and  routines.  A  flow  diagram  for  GRAPH  DRIVER 
is  provided  in  Pigure  4.4. 

b.  Subroutine  ANTL0C/TGTL3C 

ANTLOC  and  TGTLOC  are  similiar  subroutines  in 
that  they  locate  the  appropriate  antenna  and  target  storage 
arrays  filed  on  Tape  2.  Antenna  parameters  are  located  in 
Files  5-24.  Target  parameters  are  located  in  files  26-45. 
Lines  4000-4060  contain  ANTL0C.  TGTLOC  is  found  in  lines 
4100-4160.  TGTLOC  has  the  additional  function  of  defining 
the  system  antenna  parameters  for  the  current  run,  as  shown 
in  lines  4149-4153.  Thera  are  no  operator-procecessor 
'interactions  during  execution  of  these  subroutines. 

c.  Subroutine  DISAv/INVAV/AOGSAV 

DIWAV,  CN9AV  and  AUGHAV  are  similar  subroutines. 
They  drive  ACA7G  to  input  the  direct  waveform,  the  incident 
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Figure  4.4.  Flow  Diagram — INPUT  DRIVER 


waveform,  and  the  augmented  waveform.  They  allow  the 
selection  by  the  Operator  of  the  number  of  signals  to  be 
acquired  for  each  waveform  by  the  DPO  (total  composite 
waveforms  from  the  sampler)  and  averaged  by  the  OS.  2ach  of 
these  subroutines  is  accompanied  by  a  message  on  CRT 
describing  the  procedure  to  be  employed  in  acquiring  the 
particular  signal.  These  messages  are  reproduced  in 

Appendix  P.  DINA7  is  located  at  lines  3000-3170.  It  saves 
data  in  File  50  of  Tape  2.  INWAV  is  located  at  lines 
3200-3380.  It  saves  data  in  File  51.  AtJGWAV  is  located  at 
lines  3400-3590.  It  saves  data  in  File  52. 
d.  Subroutine  ACA7G 

ACA7G  acquires  and  a?erages  the  appropriate 
waveform  interactively  requested  by  the  operator  from  the 
processor  keyboard.  It  also  removes  any  dc  signal  that 
might  be  included  in  the  acquisition.  This  feature  allows 
the  operator  to  position  the  desired  signal  anywhere  on  the 
DPO  screen  and  still  maintain  a  true  zero  axis  in  later 
graphic  displays.  Note,  the  signal  waveform  must  be 
completely  within  the  DPO  screen  area.  ACA7G  is  74-steps 
long  (23%  of  the  program)  and  is  located  in  statements 
2000-2730.  It  is  driven  by  one  of  the  appropriate  input 
waveform  subroutines  previously  described.  It  returns 


results  to  that  subroutine  for  later  storage.  It  also 
displays  system  information  for  the  current  signal  being 
acquired  to  the  DPO  screen. 

Lines  2090-2130  initialize  ACA7G  by  requesting 
the  operator  to  input  the  number  of  waveform  signals  to  be 
averaged  in  M2.  This  value  is  than  returned  for 
verification.  In  lines  2140-2540,  the  OS  writes  the  current 
acquisiton  number  of  the  waveform  to  the  DPO  screen, 
acquires  and  stores  the  signal  into  DPO  memory  and  then 
transfers  the  signal  to  the  4052  for  further  processing. 

Lines  217  0-2370  perform  the  count  function. 
Depending  on  the  number  of  waveforms  desired,  the  count  will 
be  for  every  signal  (N2  less  than  21)  ,  for  every  five 
signals  (M2  less  than  101)  or  for  every  ten  signals  (M2 
greater  than  or  equal  to  101).  This  is  done  to  reduce  the 
time  required  for  acquisition  of  all  signals. 

Lines  2410-2430  instruct  the  DPO  to  store  and 
hold  the  waveform  in  DPO  memory.  Line  2420  stops  processing 
to  allow  for  the  DPO  to  fully  acquire  the  waveform  in  a 
single  sweep.  T1  is  the  length  of  the  wait. 

In  lines  2440-2450,  the  GS  instructs  the  DPO  to 
transfer  the  digitized  signal  from  DPO  memory  to  the 
microprocessor  temporary  storage  array  X0.  Lines  2470-2500 
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then  search  the  first  sight  values  of  the  array  for  any 
obviously  out  of  range  values  input  by  poor  DPO  blanking  on 
the  signal  front  end.  If  such  a  value  is  found,  the  array 
element  is  set  to  discrete  value  511,  corresponding  to  a 
scaled  value  of  0.  The  waveform  is  next  added  to  the  sun  of 
the  previous  waveforms  for  the  current  acquisition  in  line 
2530.  The  resultant  matrix  is  the  time  sampled  values  of 
the  signal  in  BO: 

N 

BO  -  Z  X0.  (4.1) 
i-1  1 

If  further  acguistions  are  necessary,  the  entire  process 
described  above  is  repeated.  If  no  further  acguisitons  are 
to  be  done  on  the  current  signal,  the  ensemble  average  tine 
sampling  matrix  X0  is  determined  in  line  2550: 


X0  -  B0/N2 


(4.2) 


Any  dc  offset  is  eliminated  in  line  2560. 

It  is  interesting  to  note  the  time  required  to 
acquire  and  average  signals  using  the  solid  state  OHP  pulse 
generator.  For  ten  averages  of  a  signal,  42  seconds  of 
elapsed  time  is  needed,  or  about  4.2  seconds  per  acquistion. 
Of  this,  about  1  second  is  needed  to  write  information  to 
the  DPO  screen  per  acquisition.  Por  21  averages,  about  65 
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seconds  of  acquisition  time  is  required,  or  about  3.1 
seconds  per  acquisition.  This  is  about  a  35$  reduction  in 
acquisition  time  overall  per  signal.  The  improved  time  per 
acquisition  is  due  to  the  reduced  frequency  with  which  data 
in  the  form  of  the  vavecount  is  written  to  the  DPO  screen. 
For  a  relatively  small  increased  investment  of  time,  a 
tremendous  advantage  is  realized  in  averaging  signals  and 
reducing  polluting  noise.  A  further  decrease  in  acquisition 
time  per  signal  is  realized  when  101  or  more  averages  are 
specified.  However,  as  will  be  discussed  in  Chapter  V,  lack 
of  total  signal  stability  may  produce  errorneous  results  for 
such  a  large  number  of  acquisitions  due  to  the  length  of 
time  required. 

In  lines  2570-2820,  the  waveform  scale  factors 
are  determined  and  transferred  to  the  processor.  The  values 
are  returned  in  (vertical  magnitude  scale  factor)  and  T ^ 
(horizontal  time  scale  factor)  .  The  discrete,  digitized 
verticle  values  are  appropriately  scaled  to  the  correct  real 
time  values  in  line  2700.  The  time  between  a  single 
discrete  sample  on  the  horizontal  axis  (sampling  rate)  is 
returned  in  Z1  at  line  2820. 

In  summary,  ACA7G  does  the  physical  acquisition, 
averaging  and  scaling  of  each  signal.  It  returns  the 
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current  waveform  count  to  the  DPO  CRT  for  display,  and  the 
appropriately  scaled  waveform  in  XO  to  the  driving 
subroutine  for  storage  in  the  appropriate  file.  It  further 
returns  the  appropriate  horizontal  and  vertical  axis  labling 
scale  factors  in  T  and  M  respectively.  The  flow  diagraa 
for  ACAVG  is  provided  in  Pigure  4.5 

The  basis  for  this  subroutine  was  derived  by 
Capt.  Haaaond.  [Ref.  19]  His  prograa  in  its  entirety  has 
been  extensively  aodified  to  fit  the  current  reguireaents  of 
the  TDRL  OS. 

e.  Subroutine  NEWPAR/ANTIN/TGTIN 

NEffPAB,  ANTIS  and  TGTIS  are  interactive 
subroutines  which  allow  an  operator  to  input  new  paraaeters 
for  defined  antennas  and  targets  into  the  aass  storage 
library  on  Data  Storage  Tape  2. 

HERP AS  is  the  driving  routine  for  ANTIS  and 
TGTIS.  It  receives  its  input  froa  the  DRIVES  and  selects 
either  ANTIS  or  TGTIN.  It  returns  the  new  paraaeters  to  the 
appropriate  file  on  the  storage  tape.  NEWPAR  then  detaraines 
if  further  inputs  are  required.  If  not,  prograa  control  is 
returned  to  DRIVER. 

ANTIN  defines  the  values  of  the  paraaeters  of 
the  new  antenna.  It  reguests  the  operator  to  input  the 
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antenna  number.  An  array  of  63  elements  is  then  generated, 
each  third  element  being  filled  with  an  initial  antenna  run 
number  equal  to  1.  All  other  elements  are  set  tc  zero.  If 
the  operator  specifies  an  antenna  value  greater  than  20,  the 
following  message  is  returned. 

"STORAGE  AREA  LI HI TED  TO  20  ANTE!! AS.  HAST  TO  COHTIHOB?” 
This  is  necessary  as  current  space  limitations  limit  the 
number  of  antennas  to  be  defined. 

TGTIN  defines  the  values  of  the  parameters  of 
new  targets.  It  requests  the  operator  to  input  the  new 
target  number  as  well  a3  the  length  and  radius  of  the 
target.  The  carraige  return  must  be  pressed  after  each 
entry. 

NEWPAR,  ANTIN  and  TGTIN  are  located  in  lines 
4200-4800.  They  comprise  11*  of  the  program, 
f.  Subroutine  SCATWAV 

SCAT  SAV  determines  the  target  backscattered 
waveform  based  on  the  difference  of  the  augmented  and 
incident  waveforms.  It  returns  the  backscattered  waveform 
to  storage  at  file  53  on  Tape  2.  It  also  sets  flag  P7  to  the 
appropriate  value  required  for  either  Prony  processing 
(P7*i)  of  Time  Domain  processing  (p7=0) .  Program  control 
is  then  returned  to  INP0T  DRIVER . 
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2.  ssaaaa  am 


MATH  performs  most  of  the  significant  mathmatical 
computations  for  the  OS.  Its  correct  performance  is  central 
in  the  production  of  accurately  processed  signals.  MATH  is 
stored  on  Tape  1,  File  3.  It  requires  14544-bytes  of  CPU 
memory  space.  With  all  arrays  and  constants  defined,  a 
maximum  total  of  49108-bytes  of  memory  are  required. 

This  sub-program  is  initiated  under  program  control 
only  from  INPUT  and  GRAPH.  There  is  no  interaction  between 
operator  and  system  during  MATH'S  execution. 

Inputs  to  MATH  are  from  INPUT,  GRAPH  or  Tape  2. 
Values  are  returned  to  mass  strorage  or  GRAPH. 

MATH  performs  the  following  functions  for  the  OS; 

1.  Past  Fourier  transforms  the  Direct  and  Baclcscattered 
waveforms. 

2.  Determines  the  frequency  domain  impulse  response. 

3.  Performs  inverse  Fourier  transform  to  find  the  time 
domain  inpulse  response. 

4.  Determines  the  step  and  ramp  response  for  the  target. 

5.  Determines  the  physical  optics  shape  of  the  target. 

6.  Optimally  filters  the  frequency  domain  impulse 

response,  yielding  a  smoothed  time  domain  impulse 
response,  ramp  response  and  physical  optics  shape  of 
the  target  under  noisy  conditions. 
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7.  Compensates  the  ramp  response  and  optics  shape  for  EH 
propagation  attenuation. 

8.  Updates  optimal  antenna  parameters  for  the  antenna 
being  used,  based  on  the  targets  ramp  response, 
automatically. 

MATH  is  composed  of  four  subroutines  and  DRIVER.  A 
description  of  each  follows.  A  general  flow  diagram  for  HATH 
is  given  in  Eigure  4.6. 

a.  HATH  DRIVER 

The  driver  routine  for  HATH  is  located  between 
lines  100  to  1210  and  comprises  59*  of  the  total  program. 
In  addition  to  initializing  the  program  and  driving  the 
subroutines,  DRIVER  performs  most  of  the  significant 
mathematical  computations  in  straight  line  data  flow 
fashion. 

HATH  DRIVER  is  initiated  at  line  100  from  INPOT 
or  line  130  from  GRAPH.  Parametric  data  is  input  directly 
from  INPOT.  It  is  recovered  from  mass  storage  in  line  160 
if  GRAPH  initiates  HATH.  Lines  100-221  serve  to  initialize 
the  local  program  parameters. 

Line  250  inputs  into  X0  from  mass  storage  file 
53  the  signal  (backscatter ed  waveform)  waich  represents  the 
convolved  transient  response  of  the  transmitting  antenna. 
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isiure  4.6.  Gentral  Flow  Diagron — MATH 


the  receiving  antenna,  and  the  target.  This  data  is  then 
used  to  drive  subroutine  FFT.  XO  is  returned  froa  FPT  in  M 
and  P  as  the  frequency  domain  aagnitude  and  phase  and  saved 
in  arrays  E  and  S3  respectively. 

A  similar  procedure  is  performed  in  lines 
290-340  for  the  convolved  transient  response  of  the 
transmitting  and  receiving  antennas  without  target  on  the 
ground  plane.  The  aagnitude  and  phase  are  returned  in 
arrays  P  and  F3  respectively. 

In  lines  360-420,  the  frequency  donain  impulse 
response  of  the  target  and  noise  optimization  are  performed. 
Magnitude  of  the  impulse  response  is  returned  in  E0,  with 
phase  returned  in  F0.  The  algorithm  for  determining  tae 
impulse  response  and  performing  optimization  were  developed 
in  chapter  III  and  will  not  be  repeated  here. 

Lines  440-540  perform  the  inverse  Fourier 
transform  of  the  frequency  domain  impulse  response  of  the 
target  to  find  its  time  domain  equivalent.  The  frequency 
domain  data  entering  at  line  440  is  in  polar  form.  This 
form  is  most  suitable  for  the  previously  described 
mathematical  manipulations.  It  is  not  a  suitable  fora  for 
the  inverse  Fourier  transform  method  employed  by  the  4052 
system.  In  order  to  perform  this  operation,  the  magnitude 
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and  phase  arrays  must  first  be  converted  into  two  source 
arrays  containing  real  and  imaginary  data.  Next,  the  first 
elenent  of  the  real  array  must  contain  the  signal  dc  tern, 
while  the  last  real  elenent  must  be  the  Nyguist  frequency 
value.  As  both  these  elenents  represent  pure  real  numbers, 
the  first  and  last  elenents  of  the  iaaginary  source  array 
muse  be  correspondingly  zero. 

Lines  440-470  establish  a  POR/NEX?  loop  that 
returns  real  terns  to  M  and  inaginary  terns  to  P.  Real 
terns  are  established  by  multiplying  the  polar  magnitude, 
E0 ,  by  the  cosine  of  the  corresponding  phase  P0: 

M(j )  »  E0 ( j )  x  (cos  (F0 ( j ) )  (4.3) 

Imaginary  terms  are  established  by  multiplying  the  polar 
magnitude  by  the  sine  of  the  corresponding  phase  tern: 

P(j)  -  E0 (j )  x  sin(F0(j})  C4-4) 

INPUT  had  previously  removed  all  dc  levels  prior 
to  NATH,  so  S(1)  and  N  (257) ,  the  dc  and  Nyguist  terms,  are 
set  equal  to  zero  in  lines  480-490.  Lines  500-510  set  the 
corresponding  imaginary  terms  to  zero.  Note  that  both  of 
these  arrays  aust  be  257  data  points  long,  as  the  output 
array  will  be  512  data  points. 
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Data  is  now  correctly  formatted  in  the  two 
source  arrays  so  that  the  command  "IULZAV"  can  interleave 
the  real  and  imaginary  data  into  a  single  destination  array, 
10,  whose  foraat  is  acceptable  to  the  GS  "IFT"  coamand. 

Line  540  calls  the  system  inverse  Fourier 
transform  command  "IFT"  which  transforms  the  spectral  data 
into  waveform  (time- doma ini  data.  The  spectral  data 
originally  in  10  is  overwritten  by  the  new  time  domain  data 
and  returned  in  10.  10  is  now  the  time  domain  impulse 
response. 

The  inverse  Fourier  transform  can  be  expressed 
mathematically  by  the  following  summation: 


X(n)  -  (1/N) 


N-l 

Z 

k-0 


Xd(k)  e 


j  2Trnk/N 


for  n-0,1, . . . ,N- 1  (4 . 5) 


In  the  above  equation,  N  refers  to  length  of  the  array 
argument,  and  n  is  an  index  used  in  generating  the  various 
elements  of  the  output  array.  X  (it)  is  the  k-th  complex 
Fourier  coefficient,  and  X (n)  is  the  n  ♦  1-th  element  of  the 
real  data  output.  [  Bef .  20]  The  4052  uses  an  extremely 
fast  method  for  IFT  computations,  the  Sande-Tukey 
decimalization- in- frequency  algorithm. 
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Line  S20-720  save  the  tine  domain  iapulse 

response  (file  54),  and  frequency  doaain  direct  waveform 

(file  55)  ,  the  backscattered  waveform  (file  56)  and  the 
iapulse  response  (file  57)  in  aass  storage. 

Sext,  the  tine  doaain  step  response  of  the 

impulse  is  returned  to  XQ  in  line  740  by  using  the  4052 

"INT"  routine.  This  is  an  integration  calculated  using  the 
trapezoidal  rule  for  a pproxinating  the  definite  integral  as 
follows: 

B  *  0 

B  (t)  »  3(t-1)  ♦  .5  *  (A  (t-1)  ♦  A  (t)  ) 

for  t  *  2,3, « «  .  ,  N 

where:  A  =  the  source  array 

3  »  the  destination  array  (result) 

N  =  the  number  of  elements  in  the  array 
Any  dc  value  in  the  result  is  removed  through  the  actions  of 
lines  750-76  0. 

Next,  the  time  doaain  raap  response  is 
determined  by  integrating  the  step  response  in  X0.  The  ramp 
is  returned  in  HO.  The  NATH  DRIVES  outputs  the  time  domain 
raap  and  its  ainiaua  value  in  R2  to  Subroutine  OPTICS  which 
returns  the  physical  optics  shape  (shape)  and  ramp  response 
compensated  for  electromagnetic  propagation  attenuation. 
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Sots  the  shape  does  not  as  yst  have  tha  correct  amplitude 
values  in  its  array  elements. 

If  the  cun  is  the  first  for  the  target  (prior  to 
optimization)  and  the  target  is  one  defined  in  mass  storage, 
MATH  DRIVER  inputs  the  last  computed  antenna  parameters  for 
the  compensator  for  shape  length  (SS)  and  radius  (S6) ,  as 
well  as  the  actual  defined  lenghth  (C6)  and  radius  (C5)  to 
subroutine  ANTDP.  S5  and  S6  are  then  automatically  updated. 


The  physical  optics 

array  is  then 

multiplied 

by 

the 

compensator  value  SS 

so  that  its 

ele 

ments  reflect 

the 

estimated  dimensions 

of  the  targets 

in 

inches. 

Thi3 

is 

based  on  the  experimentally  observed  fact  that  a  particular 
size  and  configuration  of  stationary  target  will  return  some 
backscattered  field  relative  to  its  size  and  shape.  That 
this  value  is  not  exactly  the  same  from  one  run  to  the  next, 
but  only  approximately  so,  is  due  to  certain  stastically 
determinable  effects  produced  by  noise,  small  differences  in 
orientation,  etc.  A  large  number  of  runs  should  produce  an 
averaging  value  which  takes  into  account  these  factors.  If 
the  target  is  unknown  (the  general  case) ,  or  the  current  run 
is  for  noise  reduction,  previously  computed  values  of  the 
antenna  parameters  are  recalled  from  mass  storage.  For  the 
unknown  target  these  compensators  are  the  averages  for  all 
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targets  for  all  runs  for  a  particular  antenna.  In  fact, 
OPTICS  is  not  part  of  the  data  processing  flow,  the 
compensators  having  been  previously  determined  in  INPOT  at 
lines  1024-1025. 

Following  a  known  target  run,  the  minimum  value 
of  BO  is  returned  to  B8.  If  the  run  is  not  for  optimization 
(L 1  =*  1 )  ,  then  36,  an  amplitude  compensation  optimization 
factor  is  set  equal  to  B8  and  the  actual  length  of  the  shape 
is  retained  in  Q9.  Por  graphing  purposes,  it  is  desirable 
that  B6  and  Q9  not  vary  from  optimization  run  to 
optimization  run.  HATH  DRIVER  will  input  BO,  Q9,  and  B6  to 
subroutine  OPTIH,  which  compensates  the  amplitude  and  length 
of  the  physical  optics  shape  for  the  low  pass  filtering  done 
to  reduce  noise.  An  arbitrary  target  is  treated  in  exactly 
the  same  manner.  (Rote,  this  effects  the  shape  ONLY  and  not 
the  time  domain  ramp  response  during  optimization.) 

The  above  system  has  worked  remarkably  well  in 
deducing  the  actual  physical  parameters  of  both  known  and 
unknown  targets.  For  known  targets,  accuracy  is  invariably 
10%  or  less.  For  unknown  targets,  this  accuracy 
deteriorates  to  about  15%  or  less.  If  the  data  is  treated 
in  the  raw  form,  with  no  compensation  applied,  accuracy  is 
reduced  to  25%  or  worse,  up  to  50%  for  the  unknown  target. 
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This  can  be  noted  on  the  graphic  output  of  the  shape  where 
the  compensated,  uncompensatad  and  known  (when  given)  values 
of  length  are  compared  run  to  run. 

Lines  1300-1  140  save  the  updated  antenna 
parameters  in  the  appropriate  files,  the  shape  array  in  file 
58,  the  time  domain  ramp  response  in  file  59,  and 
significant  system  parameters  in  file  60,  all  on  mass  data 
storage  Tape  2. 

Lines  1160-1190  is  MATH  DRIVER  ending  routine, 
deleting  variables  and  passing  control  of  the  program  to 
GRAPH  (file  4)  on  Tape  1. 

In  summary,  MATH  DRIVER  is  the  principal  vehicle 
for  OS  processing  for  inverse  scattering.  A  flow  diagram  is 
given  in  Pigure  4.7. 

b.  subroutine  FFT 

FFT  is  the  fast  Fourier  transform  subroutine 
used  tc  convert  the  time  domain  Direct  and  Backscattered 
waveforms  to  spectral  data.  This  subroutine  is  located  in 
lines  2000-2140  and  requires  7*  of  the  total  program  space. 

X0  contains  the  real-valued  signal  data  on  which 
the  Fourier  transform  is  to  be  performed.  It  is  input  from 
HATH  DRIVER.  Line  2060  calls  the  GS  command  "FFT"  and 
performs  a  fast  calculation  (about  4  seconds)  of  the  DFT 
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(discrete  Pourier  Transform) ,  expressed  lathematically  by 
the  summation:  [Ref.  21] 

N  - 1 

XdOO«‘l  X(n)ej2ffnk/N  ,  for  k-0,1 . N/2  (4.7) 

n*0 

where:  N  =  array  length 

ic  =  index  used  in  generating  various 
Fourier  coefficients 

Xj (k)  »  k-th  Pourier  coefficient 

X(n)  -  n  ♦  1-th  element  of  the  real  data 
input 

After  the  command  has  executed,  the  original 
real  data  will  be  overwritten  by  the  results  of  the  PFT 
computation.  The  output  data  is  in  complex  rectangular  fora 
with  the  first  element  of  the  array  XO  containing  a  real 
number  representing  the  dc  value  of  the  signal,  the  second 
element  containing  a  real  number  representing  the  value  of 
the  discrete  Pourier  transform  at  the  Nyquist  frequency,  and 
the  remaining  array  elements  containing  alternately  real  and 
imaginary  Fourier  coefficients. 

This  arrangement  is  not  suitable  for  the 
division  that  later  takes  place.  In  line  2070,  the 
rectangular  array  XO  is  converted  to  polar  form  by  calling 
the  GS  command  "POLAR".  Magnitude  components  are  returned 
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in  array  H  and  phase  components  are  returned  in  array  P.  At 
this  point,  subroutine  FFT  returns  control  to  MATH  D3I72R 
with  the  arrays  XO,  M  and  P  for  further  processing, 
c.  Subroutine  OPTIC 

Subroutine  OPTIC  manipulates  the  ramp  response 
of  the  target  to  recover  the  physical  optics  shape,  but  it 
does  so  without  effecting  in  anyway  the  raw  data.  It  later 
compensates  both  the  ramp  response  and  the  shape  for  EH 
propagation  attenuation  experienced  along  the  body  of  the 
target.  Subroutine  OPTIC  resides  in  lines  2200-2700  and 
comprises  25%  of  Program  MATH. 

The  ramp  response  RQ  and  the  location  of  its 
minimum  value  within  the  array  are  input  into  OPTIC  from 
3 ATH  DRIVER.  The  mean  value  of  the  response  is  calculated 
in  lines  2250-2260  and  returned  in  R3.  This  value  is  used 
later  in  determining  the  endpoints  of  the  shape  array. 

Lines  229  0-25  00  is  a  FOR/RE  XT  loop  that 

separates  out  the  targets  shape  from  the  ramp  array  data. 
The  search  for  valid  data  begins  at  the  point  in  the  array 
equal  to  the  location  cf  the  minimum  value.  It  is  assumed 
that  the  target  will  always  return  a  response  larger  than 

any  that  might  be  due  to  range  noise  or  other  excitation. 

For  the  size  of  the  targets  being  used,  and  the  the 
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configuration  of  the  present  source,  this  has  always  proven 
to  be  an  accurate  assumption.  The  FOR/MSX?  loop  searchs 
first  for  a  trailin g-ed ge  endpoint.  This  is  done  in 
ascending  time  steps  in  search  of  the  point  where  the  array 
elements  become  zero  or  positive,  representing  a  zero-axis 
crossing.  For  every  value  of  RO  less  than  zero,  a 
corresponding  value  is  input  in  BO,  the  shape  array.  If 
30  (k)  is  a  positive  value,  representing  a  zero-crossing,  the 
FOR/NEXT  loop  is  terminated.  If  BO  (k)  is  negative,  OPTIC 
then  checks  the  location  of  the  point.  It  is  important  to 
do  this  as  noise  and  instability  of  input  signals  from  an 
acguisition-to-acguisition  run  may  erroneously  shift-zero 
crossings  or  preclude  zero-crossings  altogether.  In  these 
cases,  determination  of  the  shape  must  be  somewhat  modified 
as  described  below. 

In  line  2350,  two  situations  are  checked. 
First,  is  the  preceding  array  element  of  RO  more  negative 
than  the  present  element  and  is  the  present  element  greater 
that  90%  of  the  ramp  response  average  value.  If 
affirmative,  the  next  array  element  is  fetched.  K  positive 
result  of  this  check  indicates  that  datum  is  on  the  skirt  of 
the  ramp  response,  and  that  the  data  magnitudes  are 
progressing  towards  a  zero-axis  crossing  normally.  If  the 
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the  present  array  value 


check  returns  as  negative,  that  is, 
is  aore  negative  than  the  preceding  element,  than  an  anomaly 
has  occurred  either  due  to  noise,  or  time-axis  shifting.  In 
this  case  a  special  endpoint  routine  must  be  employed. 

The  second  check  is  to  determine  if  the  array 
element  is  on  the  main  signal  portion  of  the  response  and 
not  the  skirts.  If  the  element  amplitude  is  greater  that 
90 %  of  the  mean  value,  than  it  is  assumed  that  the  element 
is  on  on  the  skirt  ir regardless  of  the  relationship  between 
the  amplitude  of  the  present  and  preceding  elements.  The 
next  datum  is  fetched.  It  can  be  seen  that  this  arrangement 
is  subject  to  some  error  if  a  noisy  element  is  encountered 
just  at  the  point  where  the  ramp  response  passes  the  90-th 
percentile  of  the  signals  mean  value.  This  error  is 

considered  relatively  minor. 

The  90-th  percentile  was  chosen  following 
numerous  observations  that  noise  or  time-shift  in  the  main 
ramp  response  normally  do  not  produce  fluctuations  that 
cross  this  boundary.  Zrror  in  the  shape  representation  will 
be  introduced  if  it  does,  in  fact,  cross  this  Level.  The 
result  will  be  a  fore- shortened  drawing  of  the  target. 

If  the  need  for  an  endpoint  routine  is 
indicated,  processing  goes  to  lines  2360-2470.  In  line 
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2360-2380,  the  distance  the  last  valid  data  element 
(HO  (<—  1 )  J  is  from  the  zero  axis  is  determined.  In  line 
2360,  the  distance  between  the  previous  two  valid  data 
elements  is  found  and  returned  in  w (3) : 

W  (3)  =  R0(k-2>  -  30  (1c—  1 )  (H  «) 

The  above  is  for  the  ascending  time  series  run  of  the 
FOR/NEXT  loop  (trailing  edge) .  For  the  descending  time 
series,  the  values  in  the  parenthesis  have  opposite  sign. 
This  difference  represents  an  approximation  of  the  slope  of 
the  skirt.  A  problem  with  this  approximation  is  that  it 
often  occurs  where  the  slope  is  not  as  great  as  the  average. 
The  result  is  a  somewhat  elongated  figure.  It  may  be 

possible  to  defeat  this  problem  by  talcing  more  than  two 
points  for  the  average. 

In  lines  2370-2380,  a  check  is  made  to  see  if 
the  last  valid  data  point  is  within  10-ti*es  or  less  of  the 
difference  determined  above.  If  it  is,  the  actual 
difference  is  returned  in  Q1  and  used  to  drive  the  values 
placed  in  B0  to  a  zero-crossing.  If  the  difference  is  less 
that  10%  of  the  amplitude  needed  to  effect  a  zero-crossing, 
then  one-tenth  of  the  actual  distance  is  returned  in  Q1.  In 
this  manner,  the  shape  rapidly  goes  to  zero  and  an  endpoint 


is  determined  close  to  the  estimated  real  endpoint  had  a 
natural  zero-crossing  occured.  The  endpoint  is  located  with 
a  FOR/NEXT  loop  between  lines  2420-2470.  This  FOR/NEXT  loop 
simply  iterates  the  values  of  BO  by  Q1  until  a  zero-crossing 
is  effected.  Whenever  BO  becomes  positive,  all  loop 
processing  is  terminated  and  the  endpoint  is  returned  in 
W ( 1 )  for  the  trailing  skirt  and  W(2)  for  the  leading  skirt. 

If  the  above  process  was  for  the  leading  edge  of 
the  ramp  response,  the  loop  is  reentered  and  the  trailing 
edge  check  of  the  data  is  performed.  The  method  employed  is 
identical  with  the  above  process  for  the  leading  edge.  Rote 
that  up  to  this  point,  the  ramp  response  in  RO  has  been 
totally  unaffected  by  any  processing  in  OPTICS .  It  has  been 
used  simply  as  source  of  data. 

Once  the  target  response  has  been  located  in  RO 
and  entered  in  BO,  then  both  arrays  can  be  compensated  for 
propagation  losses  experienced.  This  is  done  in  lines 
2570-2640.  The  val'*e  of  the  nominal  distance  the  target  is 
from  the  transmitting  antenna  is  returned  in  Z6.  h  FOR/NEXT 
loop  in  lines  2590-2640  corrects  the  arrays  for  the 
estimated  propagation  losses.  Line  2500  computes  the 
distance  along  the  body  the  wave  propagates  for  each  sample 
point  in  R0: 
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H  ■  (Q  -  H(1)  *  Z1  *  (10*8) /2 


(*»•<?> 

where:  (Q  -  11(1)  )  *  Z1  =  twice  the  tine  of  travel 

of  the  propagated  wave 

In  line  2610,  the  propagation  loss  compensation  factor  is 
returned  using: 

u  -  (4.10) 

22^ 

where:  z  =  range  from  transaitting  antenna  to 
target  face 

It  is  also  assumed  that  the  distance  between  transmitting 
and  receiving  antenna  is  Z.  Z  is  in  fact  the  measured  value 
of  distance  read  froa  the  ”IIHE/DIS7A»C2"  dial  on  the  DPO. 
Each  value  of  R0  and  BO  is  then  multiplied  by  this 
compensation  factor  in  turn.  Processing  continues  once  all 
values  for  the  computed  length  have  been  appropriately 
amended. 

Line  2260  is  a  FOB/NEXT  loop  that  transfers  the 
values  in  the  array  BO  to  the  beginning  of  the  array  to 
facilitate  graphing.  Once  this  is  accomplished,  program 
control  is  returned  to  HATH  DHI VEa  with  the  shape  unsealed 
for  actual  dimensions  of  the  target  (values  are  equal  in 
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amplitude  to  the  compensated  voltage  magnitudes  of  the  ramp 
response).  The  ramp  resonss  is  returned  uneffected  by  any 
of  the  above  processing  except  for  propagation  attenuation 
compensations.  Figure  4.8  provides  a  flow  diagram  for 
Subroutine  OPTICS. 

d.  Subroutine  kUTUP 

This  surboutine  updates  the  specific  antenna 
length  and  radius  compensation  values  if  the  target  is  one 
defined  from  storage  on  Tape  2.  If  the  target  is  arbitrary 
in  the  general  case,  no  updating  occurs.  ANTOP  is  located 
between  statements  3000-3080  and  comprises  4^  of  the 
program. 

The  unsealed  shape  array  B0  is  input  to  ANTtJP 
from  BATH  DR  ITER .  In  line  3060,  the  amplitude  compensating 
and  scaling  factor  is  computed  using; 

36=  (A  *■  S*  B)  /  (C*B) 

where:  S6  *  scaling/compensating  amplitude  factor 

C  =  radius  of  target 

3  *  minimum  magnitude  of  physical  optics 
array 

A  =  the  number  of  previous  times  this 
antenna/target  combination  was  used 

S  *  sum  of  previous  compensation  values 


137 


Figure  4.8.  Flow  Diagran — MATH  OPTICS 


The  result  is  simply  the  averaged  value  of 
compensation/scaling  for  all  runs  using  this  particular 
antenna  and  target  combination. 

Line  3070  performs  a  similar  computation  for  S5, 
updating  the  coapensation/scaling  factor  for  the  target 
length.  The  formula  is  identical  except  that  the  defined 
target  length  replaces  the  target  radius  and  the  observed  or 
measured  length  from  the  plane  replaces  the  minimum 
magnitude  of  the  shape  array.  The  concept  is  that  slight 
differences  in  target  returns  due  to  variations  in  target 
location  and  orientation,  noise,  etc.,  will  be  averaged  out 
over  the  long  term  for  any  particular  antenna.  »hen  the 
target  is  arbitrary  and  there  is  no  a  priori  information 
then  ANT0P  is  bypassed  and  the  average  of  all  the 
con pen sat ion /sea ling  factors  for  all  targets  for  the 
particular  antenna  is  used.  This  requires  that  any  new 
antenna  be  "seeded”  vith  coapensation/scaling  factors  from 
at  least  one  run  using  known  targets  prior  to  attempting  to 
run  an  arbitrary  target  on  the  antenna.  »hat  results  is  a 
"poor  nan's"  compensator  for  the  currents  on  the  body  of  the 
target  based  on  empiracal  observations  of  many  target 
responses.  It  is  not  as  accurate  as  the  full  solution  of 
the  time  domain  integral  equations  would  be,  nor  is  the 
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exact  target  shape  returned.  What  is  obtained  is  an 
excellent  estiaate  of  the  target,  whether  known  of  unknown, 
based  on  an  historical  compilation  of  data  for  the  TDRL 
system.  This  concept  can  be  generalized  to  any  system  using 
components  which  are  nominally  fixed,  such  as  large  scale 
bistatic  radars.  The  shape  itself,  at  least  for  simple 
axisymmetric  targets,  can  probably  be  returned  with  the 
development  of  the  proper  software  routines.  Perhaps  the 
rate  at  which  the  ramp  response  slopes  vary  may  be  analyzed 
for  this  purpose.  What  is  obtained  is  a  method  which 
employes  a  small  scale  processor  requiring  relatively  simple 
processing  to  obtain  accurate  time  domain  estimates  based  on 
the  ramp  response  of  the  target  almost  in  real  time,  and 
capable  of  being  generalized  to  almost  any  system  employing 
the  proper  source. 

The  compensation/scaling  factors  are  returned  to 
MATH  DRIVES  along  with  program  control.  BO  is  unaffected  as 
yet  by  S5  or  S6. 

e.  Subroutine  OPTIM 

OPTIM  compensates  the  shape  array  for  the 
effects  of  low-pass  filtering  of  the  target  responses  using 
Riad's  Optimization  technique.  Optimization  is  done  to 
remove  high-frequency  noise  terms  and,  for  the  shape,  to 
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improve  the  overall  shape  of  the  drawing.  In  reality, 
without  compensation,  the  drawing  is  less  representative  of 
the  true  targets  shape  when  optimized.  Low  pass  filtering 
will  lengthen  the  overall  shape  and  reduce  the  radius,  the 
amount  being  relative  to  the  value  of  lambda,  the 
optimization  variable.  OPTIH  allows  for  this  and  introduces 
factors  which  will  pin  the  optimized  shape  to  the  original 
form  while  allowing  the  benefits  (reduced  noise  ripple  on 
the  body)  of  low-pass  filter  noise  reduction  to  be  shown. 
OPTIM  is  located  in  lines  3100-3180  and  represents  4}  of  the 
program. 

BO,  B6,  (original  minimum  of  the  unoptimized 
shape) ,  B8  (the  present  minimum  of  the  shape)  and  Q9  (the 
original  unoptimized  target  length)  are  input  from  MATH 
DRIVER .  The  amplitude  OPTIM  compensation  is  found  in  line 
3140  by  dividing  B6  by  B8.  It  is  returned  in  B7.  The  new 
compensated  target  length  factor  is  returned  in  Q7  and  is 
found  by  dividing  the  present  value  of  the  length  by  the 
original  value.  Finally  the  whole  physical  optics  shape 
array  is  scaled  by  the  new  amplitude  factor.  The  new  values 
are  returned  to  MATH  DRIVER,  along  with  the  original  factors 
in  36  and  Q9.  In  fact,  after  the  initial  unoptimized  run, 
B6  and  Q9  are  constants  until  a  new  target  is  acquired. 
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GRAPH  provides  a  central  leans  to  visually  display 
the  results  of  processing  perforaed  in  HATH  .  The  correct 
evaluation  and  analysis  of  the  acquired  processed  data  is  a 
direct  function  of  the  effective  reproduction  of  the 
resultant  arrays  by  GRAPH  into  a  fora  aost  useful  and 
understandable  to  the  operator.  GRAPH  is  located  on  Tape  1, 
Pile  4 0  and  requires  4  4712-bytes  of  aenory  space. 

prograa  GRAPH  is  normally  initiated  under  prograa 
control  froa  Prograa  HATH,  in  line  100.  It  aay  be  initiated 
directly  by  the  operator  if  the  results  of  a  previous  run 
have  been  retained  in  aass  storage.  The  coaaands  issued  by 
the  operator  are: 

FIND  4 
CALL  BOLD 
BOH 

Outputs  are  to  the  CRT  or  a  return  of  control  to  one  of  the 
other  prograas  that  make-up  the  TDHL  OS.  Ho  results  are 
returned  to  aass  storage. 

GRAPH  subroutines  aay  be  divided  into  three 
categories: 

1.  DRI7BR--provides  for  initialization  of  the  prograa 

and  drives  subroutines.  It  also  passes  control  to 

the  next  prograa  in  the  OS. 
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2.  AXIS  and  Labling  subroutines--establishes  the  axis  of 
all  graphs  and  lables  as  appropriate.  Provides  for 
formatting  of  descriptive  parameters,  as  appropriate 
builds  grids  for  graphics.  Six  subroutines  comprise 
this  group.  They  are  AXIS,  TICKS,  A X LA B ,  TITLE, 
GRPAR,  GRID. 

3.  Inverse  scatter ing--Subroutines  in  this  section  graph 
the  results  of  processing  in  INPUT  and  HATH.  They 
reproduce  the  four  basic  input  waveforms,  the  Fourier 
transforms  of  the  Direct,  Bacicscatter  and  Impulse 
Response,  the  Physical  Optics  Shape  of  the  target, 
and  the  tine  domain  ramp  response.  Four  subroutines 
comprise  this  group.  They  are  GRIN,  GRFOUR,  GROP, 
and  GRAMP. 

A  general  flow  diagram  for  GRAPH  is  provided  in 
Figure  u.9.  A  disscussion  of  all  Subroutines  in  GRAPH 
follows. 

a.  Graph  DRIVER 

As  for  INPUT  and  NATH  DRIVER,  GRAPH  DRIVER  is 
the  principal  initializer  and  controling  routine  for  this 
program,  coordinating  the  flow  of  processing.  Unlike  HATH 
DRIVER,  GRAPH  DRIVER  does  not  of  itself  perform  any  data 
processing.  That  function  is  retained  for  the  subroutines 
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exclusively 


GRAPH  DRIVER  is  located  in  lines  100-1190  and 
comprises  19*  of  the  total  prograa.  GRAPH  is  coapletely 
independant  of  data  input  froa  the  other  prograas  in  the  OS. 
All  inputs  are  froa  aass  storage  Tape  2  or  the  GS  keyboard. 
All  results  are  returned  to  CRT  display  of  the  4052.  This 
routine  is  highly  operator- processor  interactive.  Appendix 
?  describes  a  typical  run  with  aessages  as  produced. 

GRAPH  DRIVER  is  always  initiated  in  line  100 
when  being  entered  external  to  the  GRAPH  Prograa.  It  is 
initiated  froa  line  290  by  returns  froa  internal 
subroutines.  lines  100-300  initialize  the  prograa  and  input 
parameter  data  froa  aass  storage.  lines  320-750  develops 
the  correct  Menu,  with  the  function/subroutine  selected 
input  to  P0  in  line  760. 

The  Menu  that  is  used  in  GRAPH  Iis  shown  in 
Apendix  ?.  Host  of  the  fiunctions  on  the  menu  are 
self-explanatory.  A  few  remarks  concerning  iteas  10,  11, 
12,  and  13  should  suffice  to  remove  any  ambiguity. 

Itea  10,  "CONTINUE  OPTIMIZATION",  initiates  the 
Riad  low-pass  filtering  optimization  of  the  impulse  response 
in  lines  950-1030.  Prograa  control  with  the  operator 
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Item  11,  "NEXT  TARGET",  passes  program  control 
to  INPUT  for  the  acquisition  of  an  entirely  new  target. 
Parameters  are  retained  and  passed  along  with  program 
control. 

Item  12,  "FINISHED  GRAPHING" ,  terminates 
processing.  The  program  can  be  reiniti'.'.ai  by  typing  "RON". 

Lines  820-900  select  the  appropriate  subroutine 
and  functions.  Line  940  returns  to  the  menu  for  further 
selection  for  all  items  except  10,11,  12,  13. 

The  next  series  of  subroutines  to  be  discussed 
perform  the  graphics  axis  formatting,  labling  and  titling. 
These  subroutines  set  the  environment  that  makes  the 
information  of  use  to  an  observer.  AXIS, TICKS,  AXLA3 ,  and 
TITLE  are  subroutines  based  on  programs  previously  developed 
by  Tektronix,  Inc.  [Ref.  22]  All  have  been  modified  to  a 
greater  or  lesser  degree  dependant  on  the  needs  of  the  TDRL. 
b.  Subroutines  AX IS/TICKS/AXLA8 

Subroutine  AXIS  drives  Subroutines  TICXS  and 
AXLA3.  They  perform  the  general  axis  definition  function 
for  the  output  graphics.  Subroutine  AXIS  serves  as  the 
imput  portal  for  data  from  the  function  subroutine  that  was 
chosen  by  the  operator.  The  results  returned  to  the  CRT  of 
the  5S  are  "neat"  ticks,  labled  and  scaled  to  the  proper 
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values  of  the  input  function.  Subroutine  AXIS  occupies 
stateaents  2000-2250,  TICKS  is  in  lines  2300-2610  and  AXLAB 
is  in  stateaents  2700-3190.  These  three  subroutines  occupy 
18“?  of  program  space. 

Data  input  from  the  function  to  be  graphed  are 
the  minimum  and  maximum  horizontal  values  of  the  viewport  in 
71  and  72,  and  the  minimum  and  maximum  vertical  values  of 
the  viewport  in  73  and  74.  These  constants  are  in  Graphic 
Display  Units  (GDO) .  GDO's  are  a  measure  of  the  resolution 
capability  of  the  GS  screen.  The  four  arguments  defined 
above  refer  to  the  actual  limits  for  the  graphic  information 
of  the  display.  The  GS  will  draw  no  lines  outside  these 
limits. 

Also  input  to  Subroutine  AXIS  from  the  driving 
function  subroutine  are  the  minimum  and  maximum  horizontal 
values  of  the  function  itself,  71  and  72,  and  the  minimum 
and  maximum  vertical  values,  73  and  74.  These  values  are  in 
User  Data  Units  (UDU's).  That  is,  in  units  that  the  user 
defines,  such  as  millivolts,  seconds,  meters,  etc. 

Finally,  AXIS  receives  the  value  of  the  flag  SO. 
This  specifies  whether  the  vertical  axis  (S0*0)  or 
horizontal  axis  (SO* 1 )  is  being  labled.  The  magnitude  of  the 
horizontal  scaling  function  for  the  input  wave,  SI  (received 
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from  INPOT)  ,  and  the  definitions  in  ODU's  for  the  vertical 
axis  and  horizontal  axis  are  required. 

«•/ 

\ 

The  perforaance  of  AXIS,  TICKS,  ar/>i  AXLAB  is 
adequately  described  in  the  Tektronix  source  referances. 
[Bef.  23]  This  report  will  describe  only  those  statements 
which  have  been  changed  or  modified  to  produce  the  desired 
output  for  TDRL  purposes. 

Principally,  the  changes  effect  the  manner  in 
which  the  axis  is  lable  with  UDtPs  in  Subroutine  AXL&3. 
These  changes  lie  in  statements  2860-3190.  Three  conditions 
aust  be  checked: 

1.  Horizontal  (S0»0)  or  vertical  axis  (SO* f )  lables? 

2.  Graphics  of  a  Fourier  transform?  (4<P0<8) 

3.  Check  amplitude  of  maximum  vertical  value  (P6)  in 

000* s. 

TABLE  17  defines  the  relationship  between 
statements  in  Subroutine  AXLAB  and  the  functions  to  be 
graphed. 

c.  Subroutine  TITLE 

This  subroutine  positions  the  title  of  the 
graphic  above  the  plotted  values.  It  i3  located  in  lines 
3200-3310  and  comprises  2%  of  the  program.  The  value  of  the 
maximum  horizontal  coordinate  32  and  the  graphic  computed 
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TABLE  IV 


Graphic  Labling  Relationships — GRAPH  AXLA3 


GRAPH  STATEMENTS 


GRIN 

GRAMP 


HORIZONTAL 

AXIS 


REMARKS 


HORIZONTAL 

AXIS 


2880-2890 


30  90 
3110 

3140-3160 

3180-3190 

VERTICAL 

AXIS 

2860 

2940-2970 
2980-2990  " 
3000-3010 
3130-3150 
3180-3190 


GRFOOR  HORIZONTAL 
AXIS 

l - 

2861-2862 

2863-2864 

2865-2866 


3180-3190 

VERTICES 

AXIS 

2  a ''0-3  170 

"3T73-3T9O 


S1*10>1000.  Lable  Axis 
S1*10<1000.  Lable  Axis 
Position  Diaension  Lable 
Diaension  Lable 
Print  Diaension  Lable 
Return 

VERTICAL  — 
AXIS 


Start 

P  (6)  <  1 .  Lable  axis 
1<P  (6) <10.  Lable  axis 
P  (6) >10.  Lable  axis 
Print  Diaension  Lable 
Return 


HORIZONTAL 

AXIS 


"Magnitude"  Axis  Label 

"Phase"  Axis  Lable 

Print  Diaension  Lable 

Return 

VERTICLE 

AXIS 

Saae  as  for  Horizontal  Axis 
Return 


NOTE:  Ail  graphic  routines  use  lines  3020-3040 
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vertical  maximum  P(6)  are  input  by  the  graphic  subroutine. 
The  graph  title  centered  over  the  graph  is  returned  to  the 
CRT . 

4.  Subroutine  GRPAS 

GRPA8  graphs  the  parameter  input  by  INPUT  and 
developed  by  MATH.  This  subroutine  is  located  in  lines 
3400-3870  and  comprises  8X  of  the  program.  All  graphs  have 
the  run  number,  date,  target  distance  from  transmitter, 
antenna  number,  target  number,  averages  taken,  optimization 
and  any  specific  remarks  printed.  The  specialized 
parameters  printed  for  particular  graphs  alone  with  the 
statement  number  applicable  are  given  in  TABLE  v. 

TABLE  V 

Specialized  Parameter  Prints — GRAPH  GRPAH 


( 

1* 


GRAPHIC 

PARAMETERS 

LINES 

GRIN 

G3AMP 

3510-3600 

MAXIMUM  PEAK  7 ALOE 

MINIMUM  PEAK  VALUE 

SMS  TALUE 

MEAN  VALUE 

G20P 

3710-3750 

DIAMETER  (INCHES) 

_ 

WIDTH  (INCHES) 
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Subroutine  GRID 


This  simple  subroutine  draws  referencing  grid 
lines  on  the  appropriate  graphic.  The  computed  graph 
maximum  and  minimum  are  input  from  the  graphic  subroutine 
along  with  the  axis  tick  intervals.  Grid  lines  are  output 
to  the  CRT  from  GRID.  Lines  3900-4010  contain  the 
subroutine.  GRID  requires  2%  of  program  space. 

This  concludes  the  discussion  of  general 
graphics  formating  subroutines.  The  driving  graphic 
subroutines  will  now  be  discussed. 

Subroutines  GRIN,  GRFOU5,  and  GRAMP  are  all 
similar  in  formatting  structure.  Kith  some  judicious  choice 
of  programming  flags  and  routine  definition,  it  should  not 
be  too  difficult  to  combine  these  three  subroutines  into  a 
single  driving  graphic  subroutine.  GROP  is  formatted 
differently,  but  still  retains  a  few  features  common  to  all 
the  graphic  driving  subroutines.  The  similarities  between 
the  subroutines  will  be  discussed  first.  The  particular 
differences  applicable  will  be  discussed  secondly. 

Each  graphic  subroutine  is  composed  of  four 
basic  blocks.  They  are: 

1.  Input  waveform  data  from  mass  btorage 

2.  Define  waveform  and  graphic  boundaries 
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3.  Draw  the  graph  of  the  function  using  special  defining 
functions  when  necessary 

4.  Labling  the  graph  with  the  correct  title  (s) 

GRAPH  DRIVER  inputs  the  function  value  PO  which 
is  used  to  identify  the  location  of  the  data  storage  file  in 
which  the  appropriate  data  rests.  This  single  row  array  is 
usually  returned  in  X.  The  waveform  window  and  viewport 
parameters,  bounding  the  output  graphic  display,  are  next 
determined.  The  maximum  horizontal  axis  values  are  found  by 
mutiplying  the  distance  between  a  single  data  point  (in 
UDO's)  by  the  total  number  of  data  points  in  the  array  along 
the  X-axis.  For  those  waveforms  listing  the  mean  graph 
values  as  part  of  the  descriptive  parameters,  the  sum  of  all 
the  array  element  values  is  obtained  and  then  divided  by  the 
discrete  number  of  array  points: 

N 

Mean  ■  W5  «  2  x(i)/N  (4.11) 

To  ensure  that  the  function  is  graphed  along  the  zero  axis 
of  the  formatted  graphic,  the  means  are  subtracted  from  the 
array  elements.  In  similar  fashion,  the  root  mean  square  is 
found: 
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(4.12) 


RMS  -  W6  -  {  Z  x(i)2/S12]1/2 
i-1 


The  aariaum  and  minimum  values  for  the  input 
data  are  found  by  implementing  the  4052  sail  "MAX'*  and  "MIN" 
routines.  These  values  are  returned  in  W4  and  43 
respectively.  The  values  are  then  output  to  the  axis 
defining  subroutines  for  establishment  of  appropriately 
scaled  and  labled  axis. 


Next, 

the 

GS  cursor  is 

positioned 

to 

the 

beginning 

of  the 

arr  ay 

with  a  "»0VE" 

coaaand. 

For 

the 

number  of 

discrete 

data 

in  the  array. 

the  function 

is 

graphed  on  the  GS  CRT.  Finally,  the  graphic  is  titled  and 
identified  by  the  appropriate  parameters.  Figure  4.10 
provides  a  general  flow  diagraa  of  the  graphic  driver 
routines. 

The  particular  aspects  of  each  graph  will  now  be 
briefly  discussed. 

f.  subroutine  GRIN 

This  subroutine  is  the  driving  graphic  for  the 
display  of  the  inputted  Direct,  Incident  and  Augmented 
Waveforms,  and  derived  Backscatter  wavefora.  It  resides  in 

statements  5000-5410  and  aakes-up  7%  of  the  program. 
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Figure  4.10.  General  Flow  Diagratt — Graphics  Subroutine 


The  general  description  of  graphic  routines 
provided  above  is  aost  typified  by  GRIN. 

Two  special  notes  need  to  be  nade.  First, 

graphing  begins  with  the  ninth  eleaent  of  the  array.  This 
is  to  prevent  the  display  of  an  erroneous  spike  input  by  the 
DPO.  (These  erroneous  values  have  a  ainiaal  effect  on  the 
processing).  Secondly,  the  function  is  graphed  in  discrete 
steps  equal  to  1/N  of  the  total  array  value  along  the  axis: 

Graphed  Point  =  (Q  *  Z1)  ,  x  (Q)  Cm. 

where:  Q  -  array  nuaber 

X (Q)  *  array  eleaent 

The  title  of  the  graph  is  selected  by  the 
routine  located  in  lines  5290-5380. 

The  horizontal  diaenion  for  the  graphic  is  tiae 
is  the  vertical  axis  is  scaled  in  volts, 
g.  Subroutine  GRFOUR 

This  subroutine  is  also  fairly  typical  of  the 
general  graphic  nethod.  It  displays  the  aagnitude  and  phase 
coaponents  of  the  spectral  data  in  polar  fora  for  the  Direct 
and  Backseat ter ed  vaveforas  and  for  the  inpulse  response  of 
the  target. 
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Only  the  first  60  elements  of  the  magnitude 
array  are  graphed.  Values  beyond  the  first  60  are  of  a 
random  nature  due  to  the  processing  performed  to  derive  the 
transforms.  spectral  content  much  beyond  10  GHZ  is, 

therefore,  negligible.  All  values  of  the  phase  array  are 
plotted,  although  the  note  concerning  the  randomness  of  the 

data  is  also  applicable  here. 

GRFOOR  is  actually  performed  in  two  distinct 
operations.  The  first  graphs  the  magnitude  of  the 

transform.  The  horizontal  axis  dimensions  are  converted  to 
frequency  from  time  by  finding  the  frequency  of  each  da ^ 
point; 

Graphed  Point  =  Q/(Z1*257),  X  (Q)  Cm.im) 

The  second  operation  graphs  the  transform  phase. 
The  horizontal  axis  dimensions  are  converted  to  radians; 

Graphed  Point  *  Q*27T)/257,  X  (Q) 

G8F0UR  is  programmed  in  lines  5500-6080.  It  is  10%  of  the 
total  program. 

h.  GBOP 

GROP  is  the  graphic  that  reproduces  the  shape 
representing  the  target.  That  is,  the  target  shape  based  on 
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physical  optics  principals  and  uncorrected  for  field  effects 
induced  by  residual  currents.  It  represents  a 
two-dimensional  cross-sectional  area  cut  of  the  target. 
Therefore,  only  target  length  and  height  are  shown. 

This  subroutine  is  constructed  generally  in  the 
same  format  as  the  other  graphic  subroutines.  One  important 
difference  is  that  a  true  one-to-one  relationship  must  be 
established  for  the  coordinate  axis.  However,  before  this 
is  done,  the  metric  dimensions  of  the  target  length  must  be 
bounded.  This  is  done  in  lines  6210-6240.  Note  that  the 
array  BO,  source  for  the  shapes  storage  in  NATH,  has  already 
been  scaled  for  the  metric  dimension  in  the  y-direction  or 
amplitude.  The  x-dimension,  or  length,  is  scaled  for  metric 
measures  in  line  6220.  This  is  accomplished  by  multiplying 
the  time  between  samples  by  the  total  number  of  samples  in 
the  array  representing  the  target.  This  value  is  then 
compensated  by  the  averaging  antenna  parameter  established 
in  HATH,  S5: 


W2 


2xl010) fS5) 


(4. 16) 


where:  Z  *  sampling  time  between  single  sasples 
L  =  measured  length  in  number  of  samples 
S5  *  antenna  length  compensator 
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Q  *  low-pass  filter  compensator 
<2  is  always  one  for  zero  optimization,  or  no  low-pass 

filtering.  It  has  effect  only  when  optimization  is 

performed,  and  is  designed  to  compensate  for  the  spreading 
of  the  target  dimensions  in  the  x-direction  during 
filtering.  Mote  also,  that  all  defined  targets  are  measured 
in  inches.  All  outputs  are  scaled  inches.  It  should  be  a 
relatively  simple  procedure  to  generalize  the  program  to 
handle  either  metric  or  english  units. 

The  next  major  difference  between  GROP  and  other 
graphic  routines  is  that  all  graphing  is  done  in  GROP  itself 
without  calling  AXIS  or  any  other  external  formatting 
subroutine.  As  mentioned  above,  the  important  concept  is  to 
properly  relate  the  co-ocdinant  axis  on  a  true  one-to-one 
plot.  This  is  done  is  statements  6290-6360  and  is  called 
"axis  scaling".  Axis  scaling  is  done  in  two  steps.  First 
the  general  x-axis  scaling  factor  AO  is  found.  This  is  done 
in  6290  by  dividing  the  maximum  horizontal  bounding  CJDtl  by 
the  viewport  established  by  the  display: 

AO  *  (M2  -  HI)  /  (V2-71-5)  (M.ll) 

This  yields  ODO's/GDO  in  the  x-direction.  The  y-axis 
scaling  factor  A1  is  then  fixed  to  this  value.  The 
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dimensions  in  the  y-direction  are  then  checked  by  the 
program  in  lines  6320  to  ensure  that  "clipping”  or  windowing 
out  of  the  data  does  not  occur.  ‘This  is  necessary  as  the 
vertical  viewing  axis  is  only  62S  that  of  the  x-axis.  If  the 
y-dimension  were  more  than  S2X  of  the  x-dimension,  data 
would  be  lost  during  graphing.  To  prevent  this,  the  x-axis 
values  are  fixed  so  that  the  maximum  x-dimension  will  always 
be  plotted.  This  is  an  effective  reduction  in  the  size  of 
the  total  displayed  shape.  To  retain  the  true  one-to-one 
axis  relationship,  the  y-direction  dimensions  must  be 
reduced  by  a  corresponding  percentage.  This  reduction 
factor  is  computed  in  line  6330  and  applied  to  the  array  in 
6440.  Note  that  in  6330  twice  the  maximum  y-dimension  is 
used.  This  is  done  because  the  graphic  output  produces  the 
mirror  image  of  a  axisymmetric  object  while  the  array  itself 
is  single-sided. 

Once  scaling  is  established  in  statement  6380, 
the  array  is  plotted.  As  data  is  graphed  it  must  be 
dimensioned  in  english  units,  compensated  for  optimization 
affects,  and  compensated  for  any  viewport  clipping.  This  is 
done  on  a  point-by-point  basis  in  line  6410.  Note  that  the 
array  has  already  been  totally  dimensioned  and  compensated. 
It  i3  in  fact  the  plotting  point  which  now  must  be  taken 
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into  consideration.  This  compensated  plotting  point  is 
returned  in  tJ  1  and  is  used  in  drawing  the  graph.  Finally, 
the  airror  image  of  each  point  is  graphed  tc  produce  a 
"true"  iaage.  A  flow  diagraa  for  GROP  is  given  in  Pigure 
4.  11. 

i.  Subroutine  GRASP 

GRASP  graphs  the  tiae  doaain  raap  response  in 
lines  6600-6910.  It  coaprises  5X  of  the  program. 

GRASP  is  identical  to  GRIN  in  practically  every 
detail.  Note  that  the  array  1  is  aultiplied  by  1000  in  line 
6670.  This  converts  array  amplitudes  to  millivolts. 
GRIN  also  employs  this  conversion. 
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2*  §XES£Ifl5fiIU  OBSERVATIONS 


In  this  chapter,  the  results  of  several  selected 
measurements  will  be  presented.  The  observations  made  were 
chosen  to  demonstrate  the  flexibility  of  the  Operating 
System  as  well  as  verify  experimental  results  with  the 
theoretical. 

A.  BEASUREBERT  1 — 7ABY  BOSESISHTED  TARGET  DISTANCE 

1.  Description  o£  setup 

The  purpose  of  this  measurement  is  to  observe  what 
effect  locating  the  target  a  distance  different  from  the 
transmitting  antenna  then  that  specified  for  the  receiving 
antenna,  but  maintaning  the  boresighted  relationship  of  the 
antennas  has  on  measurements  taken  by  the  OS.  The 
receiving  antenna  is  placed  2.54a  from  the  transmitting 
antenna.  The  target  is  defined  from  the  target  library  on 
Tape  2. 

2  •  He ^Its 

The  results  are  given  in  Figures  5.1  to  5.3  with  the 
backscattered  waveform,  ramp  response  and  physical  optics 
shapes  being  presented. 
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The  ramp  response  is  seen  to  be  typical  for  a 
cylinder  although  somewhat  reduced  in  size.  This  is  to  be 
expected  since  the  propagation  attenuation  losses  are 
greater  for  the  distance  at  which  the  target  is  located  in 
this  experiment  than  it  is  for  the  normal  control  case. 
With  the  distance  taken  info  consideration  by  the  Operating 
System  software,  it  can  be  seen  that  the  resultant  physical 
optics  shape  is  a  close  approximation.  The  target  width 
(length)  is  within  10*  of  the  true  length  and  the  diameter 
is  15.4*  of  the  true  diameter.  The  uncorrected  width  is  49* 
greater  than  the  true  length.  The  backscattered  waveform 
shows  the  effect  of  some  spreading  in  comparison  to  the 
control  case  and  also  displays  a  reduced  S/N. 

3 •  Conclusions 

The  distance  between  the  receiving  antenna  and 
target  can  be  arbitrary  and  the  Operating  System  will  be 
able  to  provide  a  reasonably  accurate  interpretation  of  the 
target.  is  the  distance  increases,  this  accuracy  can  be 
expected  to  decrease  directly  to  the  limit  of  the  Systems 
ability  to  seperate  the  target  and  noise. 
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target:  run  date  dist  ant 

2  8-7-81  2.34M  1 

NAXINUN  PEAK  UALUE . 

NINIHUN  PEAK  UALUE . 

RHS  UALUE . 

MEAN  UALUE . 

NUN8ER  OF  UAUEFORHS  AUERAGED  -  38 


TGT  REMARKS 

1  LOCATION  DIFFERENCE 
♦1.88  nU 
-1.48  nU 
♦8.82  nU 
8.88  nU 

OPTIMIZATION  UALUE  •  8 


BACK  SCATTER  HAUE  FORM 


Figura  5.1.  HE1S08EHEHT  1  —  Backscatterad  iavefora. 
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TARGET:  RUN  DATE  DtST  ANT  TGT  RENARKS 

2  8-7-91  2.54H  1  1  LOCATION  DIFFERENCE 

NAXINUH  PEAK  VALUE .  ♦6.86  nV-N 

XIHINUn  PEAK  VALUE .  -196.33  *U-H 

RNS  UALUE .  ♦112.76  nU-H 

NEAN  UALUE .  -98.84  nU-N 

NUMBER  OF  MAVEFORNS  AVERAGED  -  38  OPTIMIZATION  VALUE  -  9 

RAMP  RESPONSE  --  TINE 


Figure  5.2.  92ASUBEBEB?  1~Sasp  Response 


target:  run  date  dist  ant  tgt  remarks 

2  8-7-81  2.S4M  1  l  LOCATION  DIFFERENCE 

DIAMETER  < INCHES-CORRECTED)  6.9623  (  6.8008) 

WIDTH  < INCHES-CORRECTED)  13.2047  (12.0000) 

WIDTH  (IHCHES-UNCORRECTED)**  17.3711  (12.0000) 

NUMBER  OF  WAUEFORMS  AUERAGED  -  30  OPTIMIZATION  VALUE  -  0 

RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Mgacs  5.3,  HE1SOBSHEMT  1 —  Physical  Optics  Shape 
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B.  HEASOREHENT  2—  TARGET  AND  RECEIVING  ANTENNA  NOT 

BOBESIGHTED 

1 .  Description  of  Setup 

An  attempt  is  aade  to  observe  the  effect  of  placing 
a  target  at  an  arbitrary  location  on  the  TDRL  ground  plane. 
The  transmitted  field  is  not  axially  incident  on  the  target. 

For  this  experiment,  the  target  (cylinder  1)  was 
placed  2.11  meters  from  the  transmitting  antenna  and  at  a 
90°  angle  (approximately)  to  the  recaiving/transmitting 
antenna  axis.  The  receiving  antenna  is  located  in  the 
control  area  of  the  TDHL.  The  target  is  at  about  a  60® 
angle  to  the  receiving  antenna  and  located  2.44  meters  from 
it.  The  experimental  setup  is  diagramed  in  Figure  5.4. 

2.  Insults 

The  results  ara  given  in  Figures  5.5  to  5.7 
representing  the  ramp  response,  physical  optics  shape  and 
backscattered  waveform.  Note  that  the  distance  input  to  the 
system  is  incorrect.  The  correct  value  as  read  from  the  DPO 
should  be  1.57  meters.  The  diameter  should  be  multiplied  by 
the  factor  1.53.  Doing  so  yields  a  diameter  of  6.3341 
inches  as  the  true  corrected  value.  The  resulting 
compensated  values  are  thus  very  close  to  the  targets  actual 
dimensions.  These  figures  may  be  suspect  as  subsequent  runs 
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TftANswtTTiHfc 


Figure  5.4.  Measurement  2 — Physical  Setup 

never  yielded  the  same  degree  of  accuracy,  although  they 
approached  it.  The  important  item  is  that  the  physical 
optics  shape  is  a  close  approximation  of  that  expected  to  be 
representative  of  a  cylinder. 

3 •  Conclusions 

It  would  appear  that  the  targets  location  can  be 
generalized  if  the  shape  only  and  not  the  dimensions  is  the 
significant  item  of  interest.  However,  time  precluded 
further  verification  of  this  observation.  It  cannot  be  said 
with  absolute  certainty  that  this  observation  is  completely 
general.  Other  configurations  with  the  target  located  at 
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different  distances  from  the  boresighted  transmitting  and 
receiving  antennas  should  be  observed.  Note  that  for  proper 
compensation  for  propagation  attenuation,  the  distance  the 
target  is  from  the  receiving  antenna  is  determined  by 
reading  the  DPO  Distance  Scale.  This  should  be  the  OS 
distance  input  value  vice  the  straight  line  measured 
distance  from  transmitting  antenna  to  target. 
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TRRGET:  RUN  DATE  ojgr 

3  8-7-91  1,2?N 

ham i nun  peak  vaujb . 

MINIMUM  PEAK  VALUE . 

RHS  VALUE . . 

HE.)H  VALUE . 

whber  of  hamepqkks  averaged  ■ 


8Nt  TCT  remarks 
1  1  NOT  BORESIGHTEO 

•  •  *47. 39  mV-H 

••  “329.95  mU-h 

••  *181.31  mV-N 

1  •  “133.98  mV— N 

J _ OPTIMIZATION  value  -  0 


Figure  5.5.  Beasuresent  2 — Rasp  Response 
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TARGET:  RUN  DATE  DtST 

3  8-7-81  1.27N 

DIAMETER  < INCHES-CORRECTED) 
MIDTH  < INCHES-CORRECTED) 

WIDTH  <INCHES-UNCORRECTED)M 
NUMBER  OF  WAVEFORMS  AVERAGED  -  4 


ANT  TGT  REMARKS 
1  1  NOT  B0RE3IGHTED 

4.1447  (  8.8808) 

11.8218  <12.8808) 

13.8789  <12.0808) 
OPTIMIZATION  VALUE  •  0 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Pigor«  5.6.  Hsasarssent  2— Physical  Optics  Shaps 


target:  run  date  dist 

3  8-7-81  1.27N 

NAN I NUN  PEAK  VALUE. . 

NIHIHUN  PEAK  VALUE . 

RHS  UAUJE . 

NEAN  VALUE. . 


ANT  TGT  RENARKS 
1  1  HOT  BORESIGNTED 

..  +2.48  "U 

..  -2. 48  mV 

. .  *9. 98  nV 

..  8.00  nU 


NUMBER  OF  WAVEFORMS  AVERAGED  -  4  OPTIMIZATION  VALUE  -  0 


BACK  SCATTER  WAVE  FORM 


Figac*  5.7.  Seasure««Bt  2— BackscattarAd  Ravaf ora 


C.  3EAS0REHEHT  3— ANTENNA  EFFECTS  OR  TARGET  RESPOHSES 

1 .  Description  o£  setug 

The  purpose  of  this  measurement  is  to  verify  that 
the  effects  of  the  antenna  are  in  fact  deconvolved  out 
during  processing.  Host  previous  runs  have  been  performed 
using  antenna  1  (TEN  Horn)  .  This  measurement  is  performed 
using  three  different  length  aonopole  antennas.  The 
antennas  used  were  antenna  2  (see  TABLE  II)  (85  mm 
aonopole) ,  antenna  8  (205  mm  monopole)  and  antenna  12  (285 
mm  monopole)  .  The  target  used  vas  target  1  (see  TABLE  I). 
Both  target  and  antenna  were  placed  in  the  control 
locations,  1.27  meters  from  the  transmitting  antenna. 
21 -acquisitions  were  averaged  for  each  run.  The  results 
obtained  should  be  compared  with  those  obtained  previously 
for  the  TEH  horn. 

2 .  Results 

The  results  are  presented  in  Figures  5.7  to  5.24 
Bach  antenna  has  the  results  of  its  run  shown  in  six  graphs. 
The  first  graph  depicts  the  direct  waveform  as  seen  by  the 
individual  antenna,  the  second  is  the  ramp  response 
resulting,  the  third  figure  is  the  physical  optics 
representation  of  the  target  shape,  and  the  final  three 
figures  are  the  optimized  frequency  domain  impulse,  time 
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domain  impulse  and  ramp  response  for  the  target.  Figure  5.7 
to  5.12  are  for  uonopole  antenna  2.  Figures  5.13  to  5.18 
are  for  monopole  antenna  8  and  figures  5.19  to  5.24  are  for 
aonopole  antenna  12. 

The  first  obvious  result  is  that  in  comparison  of 
the  ramp  responses,  the  smaller  the  antenna  is  in  relation 
to  the  target,  the  less  typical  is  the  ramp  and, 
correspondingly,  the  physical  optics  shape.  However,  for 
antennas  considerably  taller  than  the  height  of  the  target, 
a  significant  amount  of  noise  pollutes  the  ramp  response. 
In  comparing  the  target  impulse  responses,  it  can  be  seen 
that  the  exact  opposite  appears  to  be  true.  The  most 
representative  impulse  is  for  the  the  shortest  monopole. 
The  frequency  domain  impulse  graphs  show  that  the  longer 
monopoles  accept  a  higer  frequency  content  than  the  shortest 
aonopole. 

3 .  Conclusions 

First,  the  closer  the  antenna  is  to  natural 
resonance  with  the  incident  pulse,  the  more  faithful  appears 
the  resulting  impulse  response  to  theory.  But  it  does  not 
naturally  follow  that  the  corresponding  ramp  response  will 
be  as  true  a  repsentation.  The  opposite  appears  to  be  true. 
It  may  be  that  the  short  antenna  is  shadowed  by  a  portion  of 
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the  target.  More  likely,  the  short  antenna  is  just  not 
large  enough  to  accept  the  full  radiated  backscatter  of  the 
target.  Conversely,  the  235-mm  monopole  is  sensitive  to 
environmental  noise.  As  can  be  seen  in  figures  5.20  and 
5.14,  the  rasp  responses  display  the  relatively  greater 
sensitivity  of  the  285-aa  monopole  to  environmental  noise 
than  the  205-mm  aonopole.  Although  the  observation  is 
limited  to  just  these  runs,  it  aay  be  possibly  concluded 
that  there  is  a  "best"  aonopole  antenna  length  for  a 
particular  target  size. 

In  conclusion,  it  can  be  said  that  the  effects  of 
the  antenna  are  in  general  deconvolved  out.  However,  other 
effects  indirectly  related  to  the  target/antenna 
configuration  may  introduce  factors  which  materially  effect 
the  reproduction  of  the  target  shape  and  its  overall 
responses. 
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TARGET:  RUN  DATE  DIST  ANT  TGT"  REMARKS 

13  8-11-81  1.27H  2  1  DIPOLE  ANTENNA  «2<85MH) 

MAXIMUM  PEAK  UALUE . ♦76.31  nU 

MINIMUM  PEAK  UALUE .  -81.33  nU 

RMS  UALUE . ♦31.38  hU 

MEAN  UALUE .  8.38  nU 

NUMBER  OF  UAUEFORMS  AUER AGED  -  21  OPTIMIZATION  UALUE  •  0 


DIRECT  UAUEFORN 


Figure  5.8.  Seasareaent  3— Direct  faveform  (85 
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TARGET:  RUN  DATE  OIST  ANT  TGT  REMARKS 

13  S-U-ai  1 . 27M  2  1  DIPOLE  ANTENNA  •2<83HH) 

MAXIMUM  PEAK  UALUE .  ♦732. 43  nU-M 

MINIMUM  PEAK  UALUE .  -539. S3  nU-M 

RMS  UALUE . . .  ♦374.80  nU-N 

MEAN  UALUE .  ♦163.31  nU-H 

NUMBER  OP  UAUEPORMS  AUERAGED  -  21  OPTIMIZATION  UALUE  -  3 


RAMP  RESPONSE  —  TIME 


Figure  5.9.  Beasareaent  3--Raap  Response  (85  aa) 


TARGET:  RUN  DATE  DIST  ANT 
13  8-11-91  1.27M  2 

DIAMETER  ONCHES-CORRECTED) 
MIDTH  ONCHES-CORRECTED) 

WIDTH  < INCHES -UNCORRECTED)** 
NUMBER  OF  WAVEFORMS  AVERAGED  -  21 


TGT  REMARKS 

1  DIPOLE  ANTENNA  «2<83MH> 
6.6879  (  6.8880) 

12.3233  02.8888) 

18.1367  02.8888) 

OPTIMIZATION  VALUE  ■  8 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 
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F igare  5.10.  Seasaresent  3— Physical  Optics  Shape  (85 
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TARGET:  RUN  DATE  0 1ST  ANT  TGT  REMARKS 

13  8-11-81  1.27M  2  1  DIPOLE  ANTENNA  *2<83MM> 

NUMBER  OF  UAUEFORMS  AUERACED  »  21  OPTIMIZATION  UALUE  -  0.3 


IMPULSE  FFT— PHASE 


Figure  5.12.  Baaaareaant  3 — Fraguency-Douain  Iapulse  Bea- 
ponse  (85  aa) 
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TARGET:  RUM  OATE 

u  a-u-ai 

MAXIMUM  PEAK  UALUE. 
MINIMUM  PEAK  UALUE. 

RMS  UALUE . 

MEAN  UALUE . 


DIST  ANT  TCT  REMARKS 

1.27M  8  1  DIPOLE  ANTENNA  *8  (28SNN> 

. +184.92  hU 

.  -116.61  nU 

.  +39.50  mU 

8.00  nU 

OPTIMIZATION  UALUE  ■  0 


NUMBER  OF  HAUEFORMS  AUERAGED  -  21 


TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

14  8-11-91  1.27(1  8  1  DIPOLE  ANTENNA  *8  <283HH> 

MAXIMUM  PEAK  UALUE .  +5.33  piU-M 

MINIMUM  PEAK  UALUE .  -889.89  nU-H 

RMS  UALUE .  >337.89  nU-M 

MEAN  UALUE .  -299.54  nU-M 

NUMBER  OF  UAUEFORMS  AUER AGED  -  21  OPTIMIZATION  UALUE  ■  8 

RAMP  RESPONSE  —  TIME 


Flyura  5.15 


Haasareaa&t  3— Ba»p  Response 
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target:  run  date  dist  ant  tgt  remarks 

14  3-1 1-81  1.27M  a  1  DIPOLE  ANTENNA  «8  (203HH) 

DIAMETER  < INCHES-CORRECTED)  3.3341  (  6.3000) 

WIDTH  (INCHES-CORRECTED)  12.0734  (12.0000) 

UIDTH  (INCHES-UNCORRECTED)**  19.3945  (12.0000) 

NUMBER  OF  WAUEFORNS  AUERAGED  -  21  OPTIMIZATION  UALUE  -  0 

RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.16.  Measurement  3— Physical  Optics  Shape  (205  ae) 


TARGET:  RUN  OATE  DIST  ANT  TGT  REMARKS 

14  8-11-91  1.27M  8  1  DIPOLE  ANTENNA  98  <203HH> 

MAX! HUN  PEAK  UALUE .  8.0898  nU 

NININUN  PEAK  UALUE .  -9.0898  nU 

RNS  UALUE .  9.9082  nU 

MEAN  UALUE .  8.0008  nU 

NUMBER  OF  UAUEFORHS  AUERAGED  »  21  OPTIMIZATION  UALUE  •  8.S 


nU 
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8.58  1.88  1.39  2.98  2.30  3.98  3.59  4.98  4.30  3.98  nS 
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Figaro  5.17.  Haasuraaant  3— riaa-Doaain  lapulsa  Basponsa 
(205  aa) 
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TARGET:  RUN  DATE  0 1ST  ANT 
14  8-11-81  1.27H  8 

NUM8ER  OF  UAUEFORMS  AUERACED  •  21 


TGT  REMARKS 

1  DIPOLE  ANTENNA  >8  <203HM) 
OPTIMISATION  UALUE  ■  8.3 


IMPULSE  FFT— MAGNITUDE 


8.88 

-3.88 


8.8  8.3  1.8  1.3  2.8  2.3  3.8  3.3  4.8  4.3  3.8  3.3  0.8  RAO 
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Figure  5.18.  Beasureient  3--Freauency-Doiain  Inpulse  Res¬ 
ponse  (2  05  an) 
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TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

14  8-11-81  1.27M  8  1  DIPOLE  ANTENNA  IS  <289X10 

MAXINUN  PEAK  UALUE .  >8. 19  nU-M 

MINIMUM  PEAK  UALUE .  -537.18  nV-M 

RMS  UALUE .  >318.88  hU-H 

MEAN  UALUE .  -295.31  nU-H 


NUMBER  OF  MAUEFORMS  AUERAGED  ■  21  OPTIMIZATION  UALUE  *  8.9 


8.98  1.88  1.98  2.88  2.98  3.88  3.98  4.08  4.38  3.88  nS 


Figure  5. 19.  Measurement  3 — Raap  Optimized  (205  mm) 


target:  run  date  dist  ant 

12  8-11-91  1.27M  12 

MAXIMUM  PEAK  UALUE . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . . . 

MEAN  UALUE . . . 
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1  DIPOLE  ANTENNA  «12<28SMH> 
♦101.31  nU 
-122.78  «U 
♦67.32  mU 
8.88  nU 

OPTIMIZATION  UALUE  -  8 


DIRECT  HAUEFORN 


Pigara  5.20.  Beasaraaent  3— Direct  latrefora  (285  aa) 
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target:  run  pate  dist  ant 

12  8-11-31  1.27N  12 

MAXIMUM  PEAK  VALUE . 

MINIMUM  PEAK  VALUE . 
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OPTIMIZATION  VALUE  »  8 


RAMP  RESPONSE  —  TIME 


?igur«  5.21.  Haasoraaant  3 — 8a«p  Haspoase  (285 
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target:  rum  date  dist  ant  tgt  REMARKS 

12  8-11-31  1.27M  12  1  DIPOLE  ANTENNA  «12<203HH> 

DIAMETER  aNCNES-CORRECTED)  6.8000  <  6.0000) 

WIDTH  < INCHES-CQPRECTED  ‘  12.8000  <12. 0000) 

WIDTH  < INCHES-UNCORRECTED)**  16.2303  <12.0000) 

NUMBER  OF  WAUEFORMS  AUERAGED  *  21  OPTIMIZATION  UALUE  »  0 _ 

RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.22.  Heasureeent  3— Physical  Optics  shape  (285 
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TARGET:  RUN  DATE  D1ST  AMT  TGT  REMARKS 

12  8- it -31  1 . 2711  i  3  1  DIPOLE  ANTENNA  il2<283NH> 

MAXIMUM  PEAK  UALUE.  . .  3.0006  nU 

MINIMUM  PEAK  UALUE .  -3.0333  nU 

RMS  UALUE .  0.3332  nU 

MEAN  UALUE.  •••  0.0000  nU 

NUMBER  OF  WAUEFORMs’ AVERAGED' ■  ' 21  *  OPTIMIZATION  UALUE  •  3.5 


nU 


3.0333 


8.3080 


-8.0033 


0.33  1.08  1.53  2.09  2.59  3.99  3.39  4.99  4.39  3.99  nS 


( 


Figure  5.23. 


Heasareaent  3 — Tiae-Doaain  Iapulse  Response 
(285  aa) 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -1963- A 


TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

12  3-11-81  1.27M  12  1  DIPOLE  ANTENNA  I12<283HH> 

NUMBER  OF  WAVEFORMS  AUER AGED  >21  OPTIMIZATION  VALUE  -  8.3 

IMPULSE  FFT— MAGNITUDE 


IMPULSE  FFT— PHASE 


C 


Figure  5.28. 


Heasureaent  3— Preqaency-Doa&in  Iapulse  Res¬ 
ponse  (285  a  a) 
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TARGET:  RUN  DATE  DIST  ANT 
12  8-11-91  1.27M  12 

MAXIMUM  PEAK  VALUE . 

MINIMUM  PEAK  VALUE . 

RMS  VALUE . 

MEAN  VALUE . . . 

NUMBER  OF  WAVEFORMS  AVERAGED  -  21 


TGT  REMARKS 

1  DIPOLE  ANTENNA  •12<28SHN> 
♦28.80  rV-H 
-396.67  rU-M 
♦368. 85  rV-M 
-314.14  rU-M 

OPTIMIZATION  VALUE  ■  8.3 


0.30  1.00  1.30  2.00  2.30  3.00  3.30  4.00  4.30  3.00  nS 


Figar*  5.25.  Saasareaent  3--Haap  Optimized  (285 
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0.  MEASURES  ENT  4— ARBITRARY  COMPOSITE  TARGET 

1 .  Description  of  Set^p 

This  measurement  was  made  to  observe  the  response  of 
the  OS  to  a  completely  arbitrary,  but  axisymmetric  target 
not  encountered  by  the  system  previously.  The  composite 
target  was  a  combination  of  target  11,  a  cone,  and  target  3, 
a  cylinder.  A  drawing  of  the  composite  is  given  in  Figure 
5.  25. 


Figure  5.26.  Heasurement  4— Composite  Targ* 

2.  assays 

The  results  are  given  in  Figures  5.26  and  5.27,  the 
physical  optics  shape  and  ramp  responses  respectively.  The 
corrected  diameter  is  11*  larger  than  actual,  and  the 
corrected  length  is  355  shorter  than  actual.  The 
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uncorrected  length  is  135  greater  than  actual.  The  shape  is 
a  close  approximation  of  the  targets  actual  shape,  although 
somewhat  shortened.  The  cone  portion  appears  to  be  about 
one-half  of  the  the  total  figure,  as  it  was  in  the  actual 
target. 

3 .  conclusions 

Based  on  this  observation  and  the  many  others  taken 
during  the  implementation  of  the  OS,  it  can  be  said  that  the 
OS  can  be  relyed  upon  to  produce  recognizable  representative 
shapes  for  known  or  unknown  axisymmetric  metallic  targets. 


195 


TARGET:  RUN  DATE  0 1  ST  ANT 
27  9-11-81  1.27N  1 

DIAMETER  < INCHES-CORRECTED) 
WIDTH  < INCHES-CORRECTED) 

WIDTH  vINCHES-UNCORRECTEDX* 
UMBER  OF  WAUEFORMS  AUERACED  ■  21 


RAMP  RESPONSE  — 


TGT  REMARKS 

???  COMPOSITE  ARBITRARY  TARGET 
6.SS69  (  0.9008) 

8.9084  (  0.8880) 

13.S3S2  <  0.0888) 

OPTIMIZATION  VALUE  -  8 


PHYSICAL  OPTICS 


?igaca  5.27.  Hsasureasat  4 — Physical  Optics  Shape 
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TARGET:  RUN  DATE  DIST 
27  8-11-91  1.27M 


ANT  TCT  REMARKS 

1  ???  COMPOSITE  ARBITRARY  TARGET 


MAXIMUM  PEAK  UALUE. 


♦38.86  nU-M 


MINIMUM  PEAK  UALUE .  -338.77  nU-M  ^ 

RMS  UALUE .  ♦  171. 95  nU-M 

MEAN  UALUE . . .  -115.28  nU-M 

NUMBER  OF  HAUEFORMS  AUER AGED  -  21  OPTIMIZATION  UALUE 


Figure  5.28.  Beasuraient  4 — Basp  Response 
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E.  SEASORBHENT  5— DIELECT BIC  CTLIMDEHS 


1.  Description  of  Satu£ 

Two  runs  were  made  using  different  dielectric 
cylinders  in  place  of  metallic  objects  as  targets.  The 
dielectric  cylinders  used  are  similar  in  configuration  to 
cylinder  targets  2  and  3  machined  from  aluainum  stock,  but 
are  aade  of  plastic.  ?or  this  aeasureaent,  the  dielectric 
cylinders  will  referei  to  as  target  2  and  3,  respectively. 

The  purpose  of  the  aeasureaent  was  to  observe  the 
response  of  the  dielectric  cylinders  to  axially  incident 
electroaagnetic  fields  and  to  also  see  how  the  OS  responds 
to  a  target  which  was  non-aetallic.  The  targets  were  placed 
in  the  standard  control  location,  1.27  meters  from  the 
transmitting  antenna.  The  runs  assumed  the  targets  were  of 
arbitrary,  unknown  3hape  for  3S  observation. 

2.  aasalia 

Target  2  was  run  first.  The  results  are  given  in 
Pigures  5.28  to  5.31.  The  physical  optics  shape  in  Figure 
5.28  is  a  very  poor  representation  of  the  true  targets 
shape.  The  reason  for  thi3  can  be  observed  in  the  raap 
response  of  Figure  5.29.  The  S/V  is  extreaely  low.  A  noise 
spike  at  about  1.6  nsecs  has  caused  the  entire  resultant 
shape  to  be  foreshortened.  The  aaplitude  discrepancy  can  be 
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attributed  to  the  fact  that  the  amplitude  compensator  was 
developed  for  the  relatively  large  valued  returns  of  the 
aetallic  cylinders,  nevertheless,  when  the  ramp  response  is 
optimized,  the  general  shape  of  the  target  can  be  seen  as  in 
Figure  5.30  and  5.31. 

The  results  using  Target  3  are  shown  in  Figures  5.32 
to  5.35.  Because  this  target  is  physically  larger  than 
target  2,  it  has  returned  a  substantially  larger  signal, 
large  enough  so  that  system  noise  is  less  of  a  significant 
factor  (although  still  seen  on  the  waveform) .  The  amplitude 
is  again  much  below  that  to  be  expected  for  a  metallic 
target  of  similiar  size. 

3 .  conclusions 

Although  the  information  provided  by  this 
measurement  is  limited,  it  appears  that  dielectric  targets 
might  provide  returns  which  can  be  identified  by  time  domain 
techniques.  This  could  be  an  important  area  for  more 
intensive  research  as  the  use  of  high  density  polymers  and 
resins  in  real  targets  of  interest  is  of  growing  importance 
to  surveillance  systems  attempting  to  locate  and  identify 
such  targets. 
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TARGET:  RUN  DATE  DIST  ANT  TGT  REMARKS 

24  9-11-91  1.27N  1  T??  DIELECTRIC  CYLINDER  <3X12  IN) 

DIAMETER  (INCHES-CORRECTED)  9.7962  (  9.9099) 

WIDTH  (INCHES-CORRECTED)  2.1396  (  9.9099) 

WIDTH  ( INCHES-UNCORRECTED''**  3.2913  (  9.9999) 

NUMBER  OF  WAUEFORMS  AUERAGED  ■  21  OPTIMIZATION  UALUE  -  9 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.29. 


Seasuceaeat  5 — Physical  Optics  (Target  2) 
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target:  run  date 

DIST 

ANT 

TGT  REMARKS 

24  8-11-31 

1.27* 

1 

???  DIELECTRIC  CYLINDER 

<3X12  IN) 

MAXIMUM  PEAK  VALUE. . 

►  •  • 

♦11.15  nU-H 

MINIMUM  PEAK  VALUE.. 

>  •  • 

-38.39  nV-M 

RMS  VALUE . 

•  e  • 

♦17.07  MU-M 

MEAN  VALUE . 

»  •  • 

-13.94  nU-H 

NUMBER  OF  UAVEFORMS  AVERAGED  • 

21 

OPTIMIZATION  UALUE  - 

0 

a. 30  1.00  1.90  2. 00  2.90  3.00  3.50  4.00  4.50  9.00  n8 


Figaro  5.30.  He&sareoeat  5  —  Haop  Response  (Target  2) 


TARGET:  RUN  PATE  PI  ST  ANT  TGT  REMARKS 

24  9-11-31  1.27M  1  ???  DIELECTRIC  CYLINDER  <3X12  IM> 

DIAMETER  (IMCHES-CORRECTED)  0. 7662  <  0.0000) 

WIDTH  < IMCHES-CORRECTED)  2.139©  (  3.0000) 

WIDTH  < INCHES-UNCORRECTED ) tS  9.3020  (  0.0000) 

NUMBER  OF  WAVEFORMS  AVERAGED  -  21  OPTIMIZATION  VALUE  -  1 _ 

RAMP  RESPONSE  --  PHYSICAL  OPTICS 


Figure  5.31.  Measurement  5 — Dptisized  Physical  Optics  (Tar 
get  2) 


TARGET:  RUN  DATE  DIST  ANT 
24  3-11-81  1.27H  1 

MAXIMUM  PEAK  UALUE . 

MINIMUM  PEAK  UALUE . 

RMS  UALUE . 

MEAN  UALUE . 

NUMBER  OF  MAUEFORMS  AUER AGED  -  21 


TGT  REMARKS 

???  DIELECTRIC  CYLINDER  <3X12  IH> 
+8.29  nU-N 


-27.71  nU-M 
♦13.28  kU-M 


-18.39  mU-M 

OPTIMIZATION  UALUE  -  I 


RAMP  RESPONSE  --  TIME 


C 


Figure  5.32.  Seasureient  5— Optimized  Baap  (Target  2) 
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TARGET:  RUN  DATE  DtST  ANT  TGT 

23  3-11-91  1.27N  1 

DIAMETER  < INCHES-C0R9ECTE0) 

WIDTH  \ INCNES-CGRPECTED ' 

UIDTH  < INCHES-UNCORRECTED) tt 
NUMBER  OF  WAUEFORNS  AVERAGED  -  2t 


REMARKS 

TT?  DIELECTRIC  CYLINDER  (EXE  IH> 
1.4134  (  0.0000) 

3.7073  (  0.0000) 

0. S719  (  0.0000) 

OPTIMIZATION  VALUE  ■  0 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.33.  Seasareaent  5~-Physical  Optics  (Target  3) 


204 


TARGET:  RUN  PATE  DIST  ANT  TGT  RENARKS 

25  8-11-81  1.27N  1  tTf  DIELECTRIC  CYLINDER  <6X6  IN> 

MAX I HUM  PEAK  VALUE. . .  >94.83  nU-H 

HI MI NUN  PEAK  VALUE. . . -71.29  «U-M 

RNS  VALUE..... . >35.88  nV-N 

NEAN  VALUE .  -3.58  nU-N 

NUMBER  OF  WAVEFORMS  AVERAGED  -21  OPTIMIZATION  VALUE  -  8 

RAMP  RESPONSE  —  TIME 


Figure  5. 3* 


He&sareeent  5  —  Raip  Response  (Target  3) 


rU 

0.0002 

e.eoee 

-0. 0002 


IMPULSE  RESPONSE 


?igora  5.35.  Haasureaant  5— Optinized  Physical  Optics  (Tar- 
gat  3) 
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Figure  5.36,  Measurement  5— Optiaized  Saep  (Target  3) 


VI.  SOHMARI 


A  software  Operating  System  has  been  developed  fcr  the 
Time  Domain  Radar  Laboratory.  The  System  nas  been  designed 
to  be  relatively  easy  for  inexperienced  operators  to  use 
after  an  initial  introduction  to  the  OS  code  and  the 
constituent  hardware  controlled.  The  OS  extends  the 
capability  of  the  TDRL  to  the  point  where  accurately 
processed  transient  time-domain  data  is  presented  in  a  form 
suitable  for  operator  evaluation.  The  OS  should  serve  as  a 
primay  research  tool  for  the  study  of  electromagnetic 
transients  at  the  Naval  Postgraduate  School. 

The  OS  is  highly  operator-processor  interactive. 
Although  the  code  that  aakes-up  the  OS  is  relatively 
complicated,  its  actual  implementation  and  use  is  not 
difficult.  There  are  numerous  directives  issued  by  the 
System  and  displayed  on  the  Graphic  Systems  Crt  to  the 
operator  requesting  interactively  the  input  of  information 
through  the  GS  keyboard  at  vital  points  in  the  programs. 

The  OS  is  composed  of  three  subprograms.  They  are: 

1.  INPOT  for  the  inputting  of  signal  information  from 
the  TDRL  ground  plane 
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2.  MATH  for  performing  the  significant  mathematical 
computations  for  the  OS 

3.  GRAPH  for  the  presentation  of  the  processed 
information  in  a  form  suitable  for  user  evaluation 
and  analysis 

INPUT  receives  the  desired  signal  in  discrete  digitized 
fora  from  the  systems  Digital  Processing  Oscilloscope,  the 
principal  analog  to  digital  converter  used  in  the  TDRL.  The 
source  signal  is  a  rapid  rise  time  gaussian  pulse  simulating 
an  impulse.  The  acguired  information  is  averaged  a  number 
of  times  as  interactively  reguested  by  the  operator.  This 
is  done  to  reduce  the  white  gaussian  noise  content  of  the 
acquired  analog  signal. 

NATH  performs  the  Fourier  conversion  and  signal 
decon volution  to  determine  the  impulse  response  of  the 
target.  It  also  implements  a  noise  reduction  algorithm  to 
optimize  the  processed  data.  The  amount  of  optimization  is 
under  the  control  of  the  operator  and  is  dependant  on  the 
degree  to  which  loss  of  small  surface  features  can  be 
tolerated  and  yet  obtain  a  good  representation  of  the 
targets  true  shape  and  dimensions.  The  step  and  ramp 
responses  are  also  obtained  and  stored  for  later  evaluation. 
MATH  determines  the  Physical  Dptics  Shape  of  the  target 
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based  on  the  rasp  response  uncorrected  for  the  effects  of 
induced  currents  on  the  body  of  the  target.  It  further 
computes  and  applies  compensation  factors  for  length  and 
width  based  on  historically  observed  results  from 
pre-defined  targets.  The  compensators  are  applied  directly 
to  unknown  targets  resulting  in  an  accuracy  of  10%  or  better 
in  determination  of  the  true  target  dimensions. 

GRAPH  outputs  the  processed  data  to  the  GS  CRT  of  the 
4052  Graphic  System.  Information  graphed  includes  the  real 
time  direct,  incident,  augmented  and  backscattered 
waveforms,  the  time-dcmain  ramp  response,  physical  optics 
shape  of  the  target  and  the  frequency  domain  spectral  data 
for  the  direct  and  backscattered  waveforms  as  well  as  the 
frequency  domain  impulse  response.  A  means  for  inputting 
the  optimization  variable  and  for  the  return  of  program 
control  to  INPUT  to  acquire  a  new  signal  is  also  provided. 

The  total  result  is  a  highly  versatile  and  powerful 
transient  analysis  tool.  The  OS  has  flexiblity  in  its 
internal  organization  to  allow  for  future  expansion  and/or 
modification  so  as  to  allow  the  control  of  constituent 
hardware  and  software  in  applications  involving  other 
aspects  of  time-domain  target  identification.  This  is 
demonstrated  in  the  complete  integration  of  a  Prony's  Method 
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Hybrid  Solution  for  the  location  of  a  targets  poles  and 
residues  with  the  requirements  for  pure  transient 
time-domain  analysis  procedures. 

The  code  itself  is  complete  as  it  stands.  However,  there 
are  areas  in  which  it  could  oe  mads  mere  efficient. 
Particularly,  a  number  of  graphic  subroutines  could  be 
combined  into  one  routine  so  as  to  reduce  memory 
requirements.  also,  a  method  for  the  display  and  selection 
of  menu  items  using  placement  of  the  screen  cursor  to  select 
the  item  might  be  developed.  This  would  allow  for  the 
display  of  more  information  graphing  processes  than  is  now 
allowed. 

One  of  the  specialized  algorithms  developed  for  the  OS 
is  a  method  for  expanding  certain  displayed  graphs  while 
retaining  a  true  plot  of  relevant  information.  The  result 
is  that  the  GS  has  essentially  infinite  resolution  and  would 
allow  for  a  close  look  at  particular  areas  cf  the  display 
without  disturbing  the  true  graph  relationships.  Although 
presently  applied  to  only  one  subroutine  in  GHAPK,  it  could 
easily  be  adapted  to  any  other  graphic  subroutine. 

A  second  specialized  routine  is  one  whereby  a  true 
one-to-one  relationship  between  coordinate  axis  is 
maintained.  The  unique  feature  of  this  algorithm  is  its 


211 


generality  in  being  applicable  to  any  size  data  field  on  any 
size  viewport  without  ever  losing  inforaation  due  to 
"viewport  clipping". 

The  results  of  field  observations  using  ■‘■he  OS  have 
proven  remarkably  faithful  to  theoretical  results.  However, 
the  present  use  of  the  OS  has  been  principally  directed 
towards  verifying  that  the  system  was  up  and  working.  As 
such,  data  was  collected  for  a xi symmetric  targets  almost 
exclusively  as  there  was  a  large  amount  of  literary  support 
in  this  research  area.  There  still  is  a  need  to  explore 
system  operation  with  this  configuration  of  targets.  In 
particular,  the  OS  needs  to  be  extended  so  as  tc  be  aole  to 
solve  the  Hagnetic  Field  Integral  Equation  for  EM  transients 
on  the  surface  of  a  target.  This  calculation  would  then  be 
compared  with  field  results  to  obtain  a  corrected  and 
accurate  representation  of  the  target  both  in  shape  and 
dimension. 

Hardware  has  proven  very  well  suited  to  the  needs  of  the 
TflHL.  However,  a  monostatic  antenna  more  directional  chan 
the  present  omnidirectional  radiator  presently  used,  may 
make  the  system  more  generally  applicable.  Improvements  in 
reception  by  better  match  of  the  systems  antenna  components 
should  allow  for  the  determination  of  finer  feature  on 
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smaller  targets.  A  more  rapid  rise  time  generator  having 
greater  power  outputs  than  presently  supplied  by  the  UHF 
impulse  generator  would  be  useful.  Finally,  a  fully 
computerized  analog  digital  processor  with  a  microprocessor 
controlled  sweep  rate  and  distance  windowing  function  would 
add  much  to  the  automation  and,  consequentially,  consistency 
of  results  in  the  TDRL. 

This  project  is  dynamic  and  holds  promise  for 
potentially  practical  applications  of  applied  research. 
There  remain  many  areas  to  explore  by  using  the  OS 
imaginatively.  The  system  is  now  ready.  The  true  need  is 
for  the  refinement  of  procedures  in  using  the  OS  and 
applying  it  to  research  areas  having  merit  both  instructive 
and  practical. 
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DESCRIPTION  DP  £2R  IPH3R  A  I,  GRAPHIC  3Y3IM  ROM 

As  stated  in  Chapter  17,  three  peripheral  ROM  packs  have 
been  added  to  the  Graphic  System  to  tailor  it  to  the 
specialized  needs  of  the  TDRL.  Th?y  are  the  Signal 

Processing  ROM  Pack  Mo.  1,  the  Signal  Processing  RCM  Pack 
Mo.  2  (FFT)  and  the  editor. 

The  Signal  Processing  ROM  Packs  provide  fifteen 
mathmatical  operations  that  are  essential  to  advanced  signal 
processing  and  data  analysis.  Signal  processing  is 
performed  on  waveforms  or  data  stored  in  one-diaensionsal 
arrays.  These  capabilities  do  no  alter  the  operation  of  the 
GS  or  occupy  any  of  its  available  RAM  space.  A  summary  of 
each  ROM  command  follows.  [Ref.  24  25] 

A.  signal  Processing  ROM  Pack  Mo.  1 

1.  MAX  --  finds  the  maximum  of  an  array. 

2.  MIN  --  finds  the  minimum  of  an  array. 


3.  CROSS  --  finis  the  location  of  a  specified  crossing 
level  within  an  array 


4.  DIP2  —  performs  a  two-point  dif f ecentiatior.  of  an 
array . 

5.  DIF3  —  performs  a  three-point  differentiation  of  an 

array . 

6.  INT  --  integrates  an  array. 

7.  DISP  —  displays  a  graph  of  an  array  in  raw  form 
(without  axis)  . 

B.  Signal  Processing  SOM  Pack  So.  2 

1.  FFT  —  computes  the  fast  Fourier  transform  of  an 
array,  placing  the  results  in  the  same  array. 

2.  IFT  --  computes  the  inverse  Fourier  transform  of  an 
array,  placing  the  results  in  the  same  array. 

3.  CONVL  —  convolves  two  input  arrays,  placing  the 

results  in  a  third  array. 

4.  COHH  —  correlates  two  input  arrays,  placing  the 

results  in  a  third  array. 

5.  POLAR  --  converts  an  array  of  fft  data  from 

rectangular  form  (reals  and  imaginaries)  to  polar 

form  (magnitude  and  phase). 

6.  TAPER  —  Multiplies  an  array  by  a  cosine  window  of 
program-selectable  tapering  weights.  When  tapering 
is  selected  to  be  50*,  it  provides  a  Hanning  window. 
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7.  ONLEAV  —  sorts  an  array  of  interleaved  FFT  data  into 
two  arrays,  one  containing  real  and  one  containing 
imaginary  components. 

8.  INLEA7  —  interleaves  the  real  and  imaginary  data 
from  two  input  arrays  into  a  third  array  whose  formal, 
is  acceptable  to  the  IFI  command. 

The  EDITOR  is  a  ROM  device  that  contains  firmware 
routines  that  allow  ASCII  magnetic  files  no  be  altered  or 
edited  off-line.  [Ref.  26] 
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APP2SDIX  3 

OSNERAL  INFORMATION  ZZ  22S.L  2EMATIN.2  SYSTEM 

A.  GRAPHING  VARIABLES 

The  following  named  variables,  numeric  and  constant.,  are 
used  in  more  than  one  subprogram.  A  brief  description  of 
their  purpose  and  the  file  in  which  they  are  initiated  is 
given  below.  Note  specific  program/subroutine  for  local 
numeric  variables,  constants  and  arrays. 


VARIABLE 

NAME 

PURPOSE 

- - — 

FILE 

INITIATED 

STRING  CONSTANTS 

E  $ 

'  1 

DATE  3F  EON 

2 

?$ 

TARGET  DISTANCE 

2 

3$ 

ANTENNA  NUMBER 

2 

HS 

TARGET  NUMBER 

2 

MS 

VERTICAL  SCALE  FACTOR  -  mVOLTS 

2 

PS 

RON  NUMBER 

2 
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VARIABLE 

NAME 

PURPOSE 

HIT"! 

*  ■+.  ±j  Zi 

INITIATED 

R$ 

"RUN'1 

.  . .  ..  .  .  . 

2 

L  - 

r$ 

_j-.m_.nni_-.-..  '  ' 

HORIZONTAL  SCALE  FACTOR  -  SECONDS 

2 

xi 

SPECIAL  REMARKS 

2 

NUMERIC  CONSTANTS  AND  VARIABLES 

36 

OPTIMIZATION  COMPENSATION 

3 

37 

OPTIMIZATION  COMPENSATION 

3 

B8 

OPTIMIZATION  COMPENSATION 

3 

C5 

TARGET  LENGTH 

2 

C6 

TARGET  RADIUS 

2 

Cl 

ANTENNA  PARAMETER  AVERAGING  VAL 

2 

L 

OPTIMIZATION  FACTOR 

3,4 

L  1 

FLAG 

3 

32 

_  _ 

SIGNALS  AVESAGED/WAV3F0RM 

2 

FLAG 

2 

P8 

RUN  NUMBER 

_ 

2 
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VARIABLE 

NAME 

PURPOSE 

FILS 

INITIATED 

37 

OPTIMIZATION  COMPENSATION 

3 

39 

ENDPOINT  DISTANCES  OPTICS 

3 

3  1 

1  1 

MAG.  TIME  SCALE  FACTOR 

2 

S4 

FLAG 

2 

35 

ANTENNA  PARAMETER-LENGTH 

3 

36 

ANTENNA  PARAMETER-WIDTH 

3 

n 

DPO  SWEEP 

2 

i  W  i  MISC  CONSTANTS  ARRAY 

L  _L 

21 

— 

SAMPLING  RATE  DPO 

2  i 

APPENDIX  C 

LIST  OP  £NPtJI  PROGRAM 

A.  GENERAL  DESCRIPTION 

1.  Location:  Pile  2 

2.  Memory  space  required  (less  mass  storage): 
24552-bytes 

3.  Maximum  no.  of  steps:  272  (212  less  remarks) 

4.  Description:  INPOT  provides  the  means  to  initialize 
and  drive  the  TD3L  Operating  System.  It  allows  input 
of  all  waveforms,  averages  acquisitions  and  stores 
data.  It  allows  for  input  of  antenna  target 
parameters  and  initializes  most  graphing  parameters. 
Inputs  are  from  the  keyboard  of  the  GS,  mass  storage, 
and  File  4.  Output  is  to  mass  storage  and  File  3. 
This  program  is  highly  operator-processor 
interactive.  Numerous  messages  and  directives  guide 
the  user. 
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3.  LOCAL  PROGRAM  TASIABLES  AND  CONSTANTS 


1  •  Significant 


—  ■ 

TYPE 

r - 1 

NAME 

OSE 

REMARKS 

ARRAY 

3 

DIHA7 

(512) — TEMPORARY  ARRAY  FOR 

DIRECT  WAVEFORM 

ARRAY 

C 

AUG  WA V 

(512) —TEMPORARY  ARRAY  FOR 

,  . 

AUGMENTED  WAVEFORM 

ARRAY 

xo 

ACAV3 

(512) — TEMPORARY  ARRAY  FOR 

INN  AV 

SIGNAL  ACQUISITION 

DIWAV 

AND  AVERAGING 

AU3WA7 

SCATNA7 

CONSTANTS 

PO 

ALL 

MENU  ITEM 

2 •  Miscellaneous 

a.  ARRAY:  30  (512) 

b.  STRING  CONSTANTS:  AS,  3S,  CS,  D$,  J$,  LS,  K$ 

S3,  Y$ 

c.  NUMERIC  VARIABLES  AND  CONSTANTS:  Cl,  C2,  C3,  C4, 
C5,  I,  MO,  N ,  N1 ,  P9,  SO,  T,  X 
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C.  MEMORY  REQUIREMENTS 


1.  4052  Internal  Data  Storage 


—  —  - 

DESCRIPTION 

3YTE3 

EACH  ACQUIRED  WAVEFORM  TO  BE  AVERAGED 

4114 

AUGMENTED  WAVEFORM  ARRAY 

4114 

DIRECT  WAVEFORM  ARRAY 

41  14 

VARIOUS  NUM  ERIC /STRING  VARIABLES  AND  CONSTANTS 

5575 

TOTAL 

17917 

MAXIMUM  MEMORY  REQUIREMENTS  (INTERNAL)  FOR  INPUT 

42469 

2.  External  Mass  Storage  Requirements  (Data  Tape) 


DESCRIPTION 

— 

BYTES 

FILE  # 

DIRECT  WAVEFORM  ARRAY 

4114 

50 

INCIDENT  WAVEFORM  ARRAY 

4114 

51 

AUGMENTED  WAVEFORM  ARRAY 

4  114 

52 

BACKSCATTERED  WAVEFORM  ARRAY 

4114 

53 

TOTAL 

_ 

16456 
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D.  SUBROUTINES 


r -  -i - 1 - 1 — - 1 - 1 - J 


CL  ASS 

NAN  E 

LINS 

STEPS 

* 

DESCRIPTION 

■  ■ 

DRIVER 

DRIVER 

100- 

1380 

129 

40 

INITIALIZE  AND  DRIVE 

PROGRAM.  PASS  CONTROL 

TO  NEXT  FILE 

INPUT 

ACAVG 

2030- 

2730 

74 

23 

ACQUIRES,  STORES  AND 

AVERAGES  DIWAV,  I  SWA V 

AUGWAV  INPUTS 

INPUT 

INWAV 

_ 

3000- 

3170 

18 

19 

6 

6 

INITIALIZES  INPUTS. 

STORES  INCIDENT  WAVE 

INPUT 

DIWAV 

3200- 

3380 

19 

6 

INITIALIZES,  INPUTS, 

STORES  DIRECT  WAVE 

INPUT 

AUG  WAV 

3400- 

3590 

20 

6 

_ 

INITIALIZES,  INPUTS, 

STORES  AUGMENTED  WAVS 

INPUT 

SCATWAV 

3600 

3630 

9 

1 

3 

INPUTS  TARGET  SCATTER 

PROM  DIFF.  OF  INWAV  - 

AUGWAV.  STORES  DATA 

— 

CLASS 

— 

NAME 

— 

LINE 

— 

STEPS 

— -  —  —  — 

— 

% 

DESCRIPTION 

SPECIAL 

ANTLOC 

4000- 

4060 

7 

2 

INPUTS  ANTENNA  PARA. 

FROM  MASS  STORAGE 

SPECIAL 

TGTLOC 

- 

4100- 

4160 

7 

2 

TARGET  DIMENSIONS 

(I?  KNOWN)  FROM  MASS 

STORAGE 

SPECIAL 

NEW  PAR 

4200- 

4330 

4 

1 

INITIALIZES  ANT IN  AND 

TGTIN  FOR  STORAGE  OF 

NEW  ANTENNA  AND 

TARGET  PARAMETERS 

SPECIAL 

ANT  IN 

4400- 

4600 

21 

7 

INPUTS  NEW  ANTENNA 

PARA.  TO  MASS  STORAGE 

SPECIAL 

TGTIN 

4700- 

4800 

1 1 

3 

INPUTS  NEW  TARGET 

PARA.  TO  MASS  STORAGE 
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APPHNDi^  D 


till  OP  MAJH  PROGRAM 


GENERAL  DESCRIPTION 

1.  Location:  Pile  3 

2.  Memory  space  required  (less  aass  storage): 
14544-bytes 

3.  Maximum  no.  of  steps:  206  (169  less  remarks) 

4.  Description:  MATH  perforins  :he  significant 

mathematical  computations  for  the  operating  system. 
These  are  fast  Fourier  and  inverse  Fourier 
transformations,  determination  of  impulse,  step,  and 
ramp  responses,  noise  optimization  and  compensation, 
and  antenna  parameter  determination.  Inputs  are  from 
Pile  2,  File  4  and  mass  storage.  Output  is  to  File  4 
and  mass  storage.  This  program  is  competely 
automated  with  no  operator-processor  interactions. 
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B.  LOCAL  PROGBAH  VARIABLES  AMD  COMSTAHTS 


1  •  Sign  if  leant 


— —  1 

TYPE 

— — — i 

NAME 

- - 

OSE 

1  . . .  . . .  "i 

REMARKS 

ARRAY 

30 

OPTIC 

DRIVER 

(512) --METRIC  PHYSICAL  OPTICS 

FROM  RAMP  RESPONSE 

ARRAY 

E 

DRIVER 

(257) — FREQUENCY  DOMAIN 

BACKSCATTE5  MAG. 

ARRAY 

EO 

DRIVER 

(257) --FREQUENCY  DOMAIN 

IMPULSE  RESPONSE  PHASE 

ARRAY 

E3 

DRIVER 

(257) --FREQUENCY  DOMAIN 

BACKSCATTER  PHASE 

AR3AY 

F 

DRIVER 

(257) --FREQUENCY  DOMAIN 

DIRECT  NAVE  MAG. 

ARRAY 

FO 

DRIVER 

(257) --FREQUENCY  DOMAIN 

IMPULSE  RESPONSE  MAG. 

ARRAY 

F3 

DRIVER 

(257) — FREQUENCY  DOMAIN 

DIRECT  WAVE  PHASE 

1  ARRAY 

10 

DRIVER 

(512) --TIME  DOMAIN  TARGET 

IMPULSE  RESPONSE 

. 
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r- 

TYPS 

— — — i 

NAME 

— 

USE 

■  . . .  . . . — ~  ■■ 

REMARKS 

ARRAY 

RO 

DRIVER 

(512) --TIME  DOMAIN  RAMP 

RESPONSE 

CONSTANTS 

MISC 

J 

ALL 

MISC  CONSTANTS  AND 

VARIABLES 

J 

2 •  Miscellaneous 

a.  ARRAY:  XO  (512)  ,  M  (257)  ,  p  (257) 

b.  NUMERIC  VARIABLES  AND  CONSTANTS:  I,  13,  X,  Ml,  M2, 
N,  N 1 ,  PI,  P2,  0,  Ql,  22,  Q4,  32,  S2,  0 

C.  M SHORT  REQUIREMENTS 

1 «  4052  Internal  Memory  Requirements 


. . .  i  -■  - —  . . . — - 

DESCRIPTION 

r———i 

BYTES 

EACH  ACQUIRED  WAVEFORM  TO  BE  AVERAGED 

4114 

FREQUENCY  DOMAIN  BACK  SCATTER  PHASE 

2074 

FREQUENCY  DOMAIN  BACKSCArTER  MAG 

2074 

FREQUENCY  DOMAIN  DIRECT  WAVE  PHASE 

2074 

FREQUENCY  DOMAIN  DIRECT  WAVE  MAG 

2074 

FREQUENCY  DOMAIN  IMPULSE  PHASE 

2074 

- - 

DESCRIPTION 


FREQUENCY  DOMAIN  IMPULSE  MAS 


EACH  ACQUIRED  WAVEFORM  TO  BE  AVERAGED 


FREQUENCE  DOMAIN  TEMP  STORAGE  PHASE 


FREQUENCY  DOMAIN  TEMP  STORAGE  MAG 


TIME  DOMAIN  TARGET  IMPULSE  RESPONSE 


TIME  DOMAIN  RAMP  RESPONSE 


PHYSICAL  OPTICS  FROM  RAMP  RESPONSE 


SCATCH  PAD  TEMPORARY  STORAGE 


VARIOUS  NUMERIC/STRING  VARIABLES  AND  CONSTANTS 


TOTAL 


BYTES 


2074 


41  14 


2074 


2074 


4114 


41  14 


4114 


4114 


1516 


34564 


MAXIMUM  MEMORY  REQUIREMENTS  (INTERNAL)  FOR  INPUT 


49108 


D.  EXTERNAL  HASS  STORAGE  REQUIREMENTS  (DATA  TAPE) 


DESCRIPTION 

1 

BYTES 

"  1  1 

FILS  # 

FREQUENCY  DOHAIN  IMP0LSE  RESPONSE 

4  las 

56 

FREQUENCY  DOHAIN  DIRECT  WAVEFORM 

4148 

55 

FREQUENCY  DOHAIN  BACK  SCATTER 

4148 

57 

TIME  DOHAIN  IMPULSE  RESPONSE 

4  114 

54 

PHYSICAL  OPTICS  FROH  RAHP  RESPONSE 

4  114 

58 

— j 

TIME  DOHAIN  RAMP  RESPONSE 

4114 

-  - 

59 

GRAPHING  PARAMETERS 

8364 

60 

TOTAL 

29036 
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iPPSgDIX  3 

LI S2  OF  GH^PH  ggOGRAM 

A.  GENERAL  DESCRIPTION 


1. 

Location:  Pile 

4 

2. 

Memory  space 

required 

(less 

mass  storage) : 

44712-bytes 

3. 

Maximum  no.  of 

steps:  584 

(4  7  3  less 

remarks) 

4. 

Description: 

GRAPH  prov 

ides  the  means  to  visually 

display  the 

results  of 

signal 

acquisition  and 

process ing. 

Subroutines 

are  div 

ided  into  three 

classes: 

a)  GENERAL — those  subroutines  th'.t  are  applicable  to 
most  of  the  other  routines  in  GRAPH.  They 
initialize  the  program,  format  graphics,  title, 
lable,  and  overlay  a  grid. 

b)  INVERSE  SCATTERING — those  routines  applicable  to 

the  inverse  scattering  solution  of  time  domain 
electromagnetic  transients.  They  are  the 

secondary  drivers  for  the  development  of  graphics. 
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c)  PHONY'S  TEC  HNIQUE-- routines  applicable  to  the 
Prony  method  of  solution. 

Inputs  to  GRAPH  are  from  File  3,  File  7,  the  GS  keyboard 
and  mass  storage.  Outputs  are  to  the  GS  CRT,  mass  storage. 
File  2,  File  3,  File  6,  and  File  8.  This  program  is  highly 
operator-processor  interactive.  Numerous  messages  and 
directives  guide  the  user.  Signicant  graphs  are  provided 
for  operator  evaluation  of  processing  results. 

B.  LOCAL  PROG HI 8  VARIABLES  AND  CONSTANTS 


1 •  £A££i£icanf 


r~  - ~i 

TYP3 

HAME 

OSE 

REMARKS 

ARHAY 

X 

VARIOUS 

(512) --TEMPORARY  STORAGE 

AND  PROCESSING 

ARRAY 

r 

VARIOUS 

(512) --TEMPORARY  STORAGE 

AND  PROCESSING 

ARRAY 

p 

AXIS 

TICKS 

AXLAB 

(8)  —AXIS  LABLING 

. . . .  , 

250 


TYPE 

NAME 

USE 

• - 

REMARKS 

CONSTANTS 

SO 

'  1  J  1 

ALL 

FLAG 

CONSTANTS 

LS 

ALL 

HORIZONTAL  AXIS  LAB LE 

CONSTANTS 

NS 

ALL 

VERTICAL  AXIS  LABIS 

CONSTANTS 

PO 

ALL 

MENU  ITEM 

CONSTANTS 

L2 

DRIVES 

OPTIMIZATION  FACTOR 

CONSTANTS 

AO 

GSOP 

HORIZONTAL  SCALING 

FACTOR 

CONSTANTS 

A  1 

G80P 

VERTICAL  SCALING 

FACTOR 

CONSTANTS 

n 

GROP 

SCALING  VIEMPORT 

COMPENSATOR 

.  _ 

2 •  Miscellaneous 

a.  AS  8  M :  none 

b.  STRING  CONSTANTS:  AS,  3$,  S$ 

C.  NUMERIC  VARIABLES  AND  CONSTANTS:  D5,  Gl,  G2,  G3, 
G«,  C5 ,  I,  J,  Q,  PI,  P2,  P3 ,  P4 ,  VI,  V2,  V3,  V4, 
Ml,  M2,  M3,  M4 ,  M5,  M6 
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C.  MEMORY  REQUIREMENTS 


1  •  405  2  I  lit arna  1  Data  Storage 


1  — 1 

DESCRIPTION 

r— - 1 

BYTES 

VARIOUS  TEMPORARY  STORAGE /PROCESSING 

3228 

MATRIX  LA3LING  ARRAY 

82 

VARIOUS  NUN  ERIC/STRING  VARIABLES  AND  CONSTANTS 

1566 

TOTAL 

9866 

MAXIMUM  MEMORY  REQUIREMENTS  (INTERNAL)  FOR  INPUT 

54573 

2.  fiitiiaal  mss  stosaai.  R§aui £SJ»eats  (Data  Ta£e)  :  none 


D.  SOBROOTI BBS 


CLASS 

r~ - 

NAME 

LINE 

STEPS 

% 

—  —  "  ~~i 

DESCRIPTION 

DRIVER 

DRIVER 

100- 

1190 

110 

19 

INITIALIZE  AND  DRIVE 

PROGRAM.  PASS  CONTROL 

TO  NEXT  FILE 

GRAPH 

AXIS 

2000- 

2250 

26 

4 

DRIVES  AXIS  GRAPHING, 

GRAPH 

TICKS 

2300- 

_ 

22 

5 

AXIS  TICKS  INTERVALS 

25  2 


CLASS 

— 

NAME 

LINS 

-  -  .  J 

STEPS 

r— — - - -  —  — 

DESCRIPTION 

GRAPH 

AXLAB 

n  - 1  ^ 

2700- 

3190 

50 

9 

LA3LE  VERTICAL  AND 

HORIZONTAL  AXIS 

GRAPH 

TITLE 

3200- 

3310 

12 

2 

POSITION  GRAPH  TITLE 

GRAPH 

GRPAR 

_ 

3400- 

3870 

4  8 

8 

PRINTS  PARAMETER 

VALDES  AS  REQUIRED 

GRAPH 

GRID 

3900- 

4010 

12 

2 

DRAWS  GRIDS 

INVERSE 

SCATTER 

GRIN 

.  _ 

5000- 

5410 

42 

7 

GRAPHICS  FOR  DIRECT, 

INCIDENT,  AUGMENTED 

AND  3ACXSCATTER 

TIME  DOMAIN  WAVES 

INVERSE 

SCATTER 

GRFOOR 

5500- 

6080 

59 

10 

GRAPHICS  FOR  FREQ 

DOMAIN  PHASE 

AND  MAGNITUDE  OF 

DIRECT,  BACKSCATTER 

AND  IMPULSE  RESPONSE 

INVERSE 

SCATTER 

GROP 

6100- 

5410 

43 

8 

GRAPHICS  FOR  PHYSICAL 

OPTICS  OF  TARGET 

INVERSE 

SCATTER 

.  ■  ■  . 

GRAHP 

- - 

6600- 

6910 

32 

5 

GRAPHICS  FOR  RAMP 

RESPONSE 
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— 

CLASS 

r— - — 

NAME 

r— - — — — 

LINS 

p  i 

STEPS 

r— - 

% 

DESCRIPTION 

PROMT 

PRIMP 

7000- 

7340 

35 

6 

GRAPHICS  FCR  IMPULSE 

RESPONSE 

PHONY 

PROTEAN 

7400- 

7710 

32 

5 

GRAPHICS  FOR  PHONY 

TRANSIENT  RESPONSE 

PHONY 

PROLES 

7800- 

8300 

51 

_ 

9 

_ 

GRAPHICS  FOR  FRONT 

POLES  OF  IMPULSE 
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APPENDIX  F 

TYPICAL  0PE3A2INS  SYSTEK  PROCEDURE 

In  this  Appendix  is  presented  a  typical  run  in  which  a 
sphere  and  a  cone  are  illuminated  by  the  incident 
elect omagnet ic  field.  The  messages  presented  and  the 
results  displayed  are  those  as  would  be  encountered  during  a 
normal  run  of  the  system. 

The  OS  is  first  used  to  acquire  an i  display  a  normal 
signal  (the  sphere)  .  This  signal  is  than  optimized  to 
remove  noise.  The  results  of  the  effected  waveforms  are 
then  represented.  Finally,  the  system  is  returned  to  INPUT 
and  a  new  target,  a  cone  is  acquired.  The  results  for  the 
cone  are  presented.  Note  that  the  impulse  response  is 
optimized  by  a  factor  of  .5  while  all  other  graphs  for  the 
cone  are  shown  unoptimized. 
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